[Nek5000-users] Subroutine opdssum

[nek5000-users at lists.mcs.anl.gov]

Dear Andreas,

gsh_fld(ifldt) refers to the gather-scatter handle for field "ifldt"

ifldt would be a flag refering to the field (1=velocity, 2=temperature, 
etc.) for which the gather-scatter (gs) operation is being performed.

At present, the situation is not entirely satisfactory because ifield
is a global variable that can be in an ambiguous state.  In the past,
opdssum (like all other "op" commands) would be applied only on a vector
field - for which there is only one mesh (or topology), namely, velocity.
In MHD, one also has the B-field but this field generally matches the
topology of the velocity.   This state of affairs will change in the near
future when we support a B-field external to the velocity field.

I hope this clarifies a bit...

Paul



On Sat, 28 May 2011, nek5000-users at lists.mcs.anl.gov wrote:

> Dear all,
>
> I was looking into the subroutine "vec_dssum", which is called in "opdssum".
> At a certain point there is a call to "gs_op_many" (a gather scatter operation I guess, isn't it?) and one of the arguments which is passed to this subroutine is "gsh_fld(ifldt)". I printed to screen the content of  "gsh_fld" but I don't understand what is the meaning of it.
>
> Thanks in advance for your help.
>
> Regards,
> Andreas.
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>