Re: Using external subroutines in OpenMP?

On Apr 28, 5:21 pm, jomarbue...@[EMAIL PROTECTED]
 wrote:
> On Apr 28, 11:57 am, Tobias Burnus <bur...@[EMAIL PROTECTED]
> wrote:
>
>
>
> > On Apr 28, 6:26 pm, jomarbue...@[EMAIL PROTECTED]
 wrote:
>
> > > > > In the parallel ****tion, I tried making private copies of the
array
> > > > > that the subroutines change. I also tried making private copies
of all
> > > > > arrays the subroutines use. Neither approach worked.
>
> > Sometimes compiler decide to use static memory for local variables in
> > subroutines. You could try to mark those subroutines as RECURSIVE.
>
> > In gfortran -fopenmp implies -frecursive, without the latter option,
> > larger local arrays might end up in static memory (rather than on the
> > stack), which can cause the problems [cf. -fmax-stack-var-size].
>
> > I think most if not all compilers don't use static local variables
> > when OpenMP is enabled. However, I believe you compiled the netlib
> > libraries without -fopenmp (ifort: -openmp). Thus you should recompile
> > those libraries with -frecursive (ifort: -recursive) or you mark the
> > subroutines of interest as RECURSIVE and recompile.
>
> > If this does not help or you already compiled the libraries with the
> > needed options: Good luck at finding the true cause of the problem.
>
> > Tobias
>
> Dear Tobias,
>
> Thank you for your response. It seems that you hit the nail in the
> head. I re-made the libraries using the -recursive flag and the
> program no longer crashes or hangs in the subroutine were I found the
> problem. Thus far it seems that the combination or making all the
> arrays shared but giving each threads a column of the array and using
> the "-recursive" flag works. I have another subroutine that still uses
> a private array and it still crashes. I'll modify that second
> subroutine in the way I did the first.
>
> Thanks again,
>
> Jomar

Hi Tobias,

I spoke too soon. I had a couple of runs where my first subroutine was
succesfull. Since then it's been very erratic. Even the serial
****tions of the subroutine are not working.

Thanks again for your response

Jomar