Hi,
In all linear eqtn solvers there is an initial stage where memory allocations take place given the "dimensions" of the problem. This is very beneficial for using the same memory (as in a PDE setting) for reusing the same solver without paying the penalty of memory allocations.
With direct sparse solvers though, I see that the very rhs and sparse matrix are used for initialization.
Since I have overwritten the mkl allocator with another one (tbb), I tried to see, if any memory allocations happen at all, by putting a breakpoint, but got nothing. Since I never had an opportunity to check my overwritting of the MKL allocators-still looking for a function totest it-, though, I cannot be sure.
Can you please confirm, if this is the case, that is, that no memory allocations happen internally in the direct sparse solver ?
TIA for your help,
Petros
Update: it's been a little while I asked this. anyone noticed ? ;- ))