Re: Who uses clapack?

From: Salvatore (sfilippone_at_uniroma2.it)
Date: 12/11/04 

Date: 11 Dec 2004 08:35:24 -0800

Victor Eijkhout wrote:
> Ok, it's clear where this use of temporaries comes from. However,
it's
> not clear that dynamic allocation is the best solution. If your
routine
> is called once, with a large dataset, then it could do its own
> alloation. However, if the routine is called a large number of times
> with the same small problem size, then you want to do the allocation
> outside it.
>
> We probably have to take a step back to the application level, and
see
> which routines can be part of a computationally intensive inner loop.
>

If you do write a shell in Fortran 95 it would be easy to have an
optional WORK parameter that allows the user to do his/her own
allocation when needed, and the shell will allocate its own otherwise
before calling the inner code (be it legacy F77 or newly rewritten F95
or whatever).

Cheers
Salvatore

Relevant Pages

  • Re: Who uses clapack?
    ... > which do not support dynamic memory allocation? ... not clear that dynamic allocation is the best solution. ... if the routine is called a large number of times ... www cs utk edu tilde lastname ...
    (comp.lang.fortran)
  • Re: Who uses clapack?
    ... > which do not support dynamic memory allocation? ... not clear that dynamic allocation is the best solution. ... if the routine is called a large number of times ... www cs utk edu tilde lastname ...
    (sci.math.num-analysis)
  • Re: Who uses clapack?
    ... > not clear that dynamic allocation is the best solution. ... if the routine is called a large number of times ... in form of an allocatable array*, and the called routined checked whether it ... it would allocate at that proper size. ...
    (sci.math.num-analysis)
  • Re: Crash of dlincg
    ... is not allocation for the return inverse but the working space inside ... the routine itself so allocating that memory a priori isn't going to ... memory for the work arrays outside the routine and pass it by using the ... application to be able to allocate this work space, ...
    (comp.lang.fortran)
  • Re: Who uses clapack?
    ... > not clear that dynamic allocation is the best solution. ... if the routine is called a large number of times ... in form of an allocatable array*, and the called routined checked whether it ... it would allocate at that proper size. ...
    (comp.lang.fortran)
Quantcast