Poll

Maximum number of variable "valp" one can use in user element

Maximum number of variable "valp" one can use in user element
0 (0%)
no
1 (100%)

Total Members Voted: 1

Author Topic: Maximum number of variable "valp" one can use in user element  (Read 3563 times)

Dadadong

  • New Member
  • *
  • Posts: 2
Currently I have 18 variables to output, and I did like this in the user element:
          valp(1)=rho
          valp(2)=rho_ele
          valp(3)=rho_hol
           .........
          valp(17)=sigma(2)
          valp(18)=sigma(3) 
Now the output file works just fine, however when feap proceed to get vtk file for paraview, such an error occurs: *** Error in `feap': corrupted double-linked list: 0x0000000001113ce0 ***
when I change back to only 9 valps, it works again.
I did change the "stress" options in umacr04, like this:
          if(foundnp) then
            write(17,'(a,2i10,a)') 'stress',12,numnp,' float'
            do n=1,numnp
              write(17,'(12f30.11)') (hr(begnp-1+i*numnp+n),i=1,12)
            enddo
          endif
and all these efforts didn't work.

So, is there a way to give more than 9 valps?

Prof. R.L. Taylor

  • Administrator
  • FEAP Guru
  • *****
  • Posts: 2649
Re: Maximum number of variable "valp" one can use in user element
« Reply #1 on: April 01, 2017, 10:39:47 AM »
In your user element you need to set two parameters to increase the number of projected items.  To get the parameters being set include 'strnum.h':

1.  In the 'isw = 1' part:  set   istv = number you want to project (12?)

2. In the projection routine set 'iste = number being projected for this element'

This will allocate the memory feap needs for the projections.

Dadadong

  • New Member
  • *
  • Posts: 2
Re: Maximum number of variable "valp" one can use in user element
« Reply #2 on: April 02, 2017, 08:25:26 AM »
Dear Prof. Taylor,

Thank you very much for the reply, it did easily solve my problem!

Best regard!