Welcome, Guest
Username: Password: Remember me

TOPIC: How to run MASKOB in parallel mode

How to run MASKOB in parallel mode 11 years 8 months ago #5126

  • fcebron
  • fcebron's Avatar
Hello

I want to run a model in parallel mode, using the MASKOB routine :

MASKEL(I)=0.D0, where I the the index of the element taht I want to mask.


How can I transform the global index of my elements to the local index in each sub-mesh ?

For the point, I konw the MESH%KNOLG%I(I) instruction, but what is the good one for element ?

Regards
The administrator has disabled public write access.

Re: How to run MASKOB in parallel mode 11 years 8 months ago #5132

  • riadh
  • riadh's Avatar
Hello,

You can use the following integer arrays:
MESH%KNOLG%I : gives the original global number of nodes
MESH%KNOGL%I : gives the local number of node as a function of their original number
This mens that MESH%KNOGL%I(MESH%KNOLG%I(J))=J

You can see the description of these arrays in the commented lines of routine paragl.f (bief)

With my best regards.

Riadh
The administrator has disabled public write access.

Re: How to run MASKOB in parallel mode 11 years 8 months ago #5133

  • fcebron
  • fcebron's Avatar
Thanks....

MESH%KNOGL%I gives the local number of node as a function of their original number

But I need a function that gives thne local number of ELEMENT as a function of their original number, since i want to use the MASKEL .....
The administrator has disabled public write access.

Re: How to run MASKOB in parallel mode 11 years 8 months ago #5137

  • jmhervouet
  • jmhervouet's Avatar
Hello,

There is no such information for elements, I'm afraid you will have to do a loop on all elements, to find the one that has the correct 3 global numbers of points, and there you will have to use KNOLG. Depending on where you do this, you can use the connectivity IKLE as IKLE(NELEM,3), or if not available as such the 3 global numbers of points of the element will be MESH%IKLE%I(IELEM),MESH%IKLE%I(IELEM+NELEM),MESH%IKLE%I(IELEM+2*NELEM), this for the local numbering, and then apply KNOLG to find the original numbering.

I hope this helps, it is rather cumbersome but actually the KNOLG for elements is never used except for your case.

R>Egards,

Jean-Michel Hervouet
The administrator has disabled public write access.
Moderators: pham

The open TELEMAC-MASCARET template for Joomla!2.5, the HTML 4 version.