open TELEMAC-MASCARET The mathematically superior suite of solvers

Hi
I would say yes
The better think is to check this by printing the information for the first element...

Regards

Thanks a lot!

I defined a BIEF_OBJ in a subroutine as

TYPE(BIEF_OBJ) :: EVOL

It passed the compilation but failed to use OS('X=0 ',X=EVOL) on it or assign values to EVOL%R. Do I have to initialize it in this subroutine before using it? I know for other global BIEF_OBJ, They are initialized in BIEF_ALLVEC(...).

Kind regards,
Qilong

Hi Qilong

Yes of course you should allocate it.
It's like all other vector types you declare it and you allocate it.
If you add a call to bief_allvec with the right parameters this should works.

Regards

Thanks for your suggestions.

I also printed out the IKLE vector. The order is not like what I thought before.

Actually it's like

[ELEM_1(1) ELEM_2(1) ... ELEM_N(1) ELEM_1(2) ELEM_2(2) ... ELEM_N(2) ...]

I also found a piece of code that confirms this.

Kind regards,
Qilong

Hello,

In Fortran in an array with two dimensions, in the memory the first dimension varies first, before changing the second, so e.g. IKLE(1,1), IKLE(2,1) up to IKLE(NELEM,1) and then IKLE(1,2), etc.


In C language it is the contrary, too bad!

With best regards,

Jean-Michel Hervouet

Thanks for your explanation!
It's much clear to me now.

With best regards,
Qilong