Re: Using external subroutines in OpenMP?

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