[aspect-devel] Direct Solver in Aspect

[Jonathan Perry-Houts]

-----BEGIN PGP SIGNED MESSAGE-----
Hash: SHA1

So, I guess the direct solver isn't working because of something about
enumerating block vectors and index sets when pressure and velocity
are in the same block. That would be easier to debug if the vectors
are always organized that way, rather than making the direct solver a
special case.

Is there a technical reason why pressure is in a separate block from
velocity or can they be grouped in both the direct solver case and the
default, iterative solver, case? It would eliminate some special cases
and tests when setting up the introspection object, and make it easier
to verify that the correct indices are being accessed.

Cheers,
Jonathan

On 05/13/2014 10:47 AM, Wolfgang Bangerth wrote:
> 
>>> I have no idea how difficult it would be to support any of this
>>> (I don't even recall who wrote the current interfaces...)
>> 
>> 1. configure Trilinos with it 2. change one line in SolverDirect.
>> :-)
> 
> Interesting. Maybe we should change the installation instructions
> for Trilinos to make this possible in the future.
> 
> 
>>> Yes, that occurred to me as well after I sent the last email.
>>> We currently drive the MPI_InitFinalize class with only a
>>> single thread, but maybe we should be using 
>>> n_cores/n_mpi_jobs_on_this_host instead? This would give us
>>> exactly the behavior you describe, while reverting back to the
>>> current state on machines that we fully load with MPI 
>>> processes...
>> 
>> Yes. I was going to try things out first though. Maybe this could
>> be a runtime parameter...
> 
> That's too late -- we need to know at the time we call the 
> MPI_InitFinalize constructor. At least I haven't found a way yet
> to change the number of threads TBB creates at a later time...
> 
> Best W.
> 
> 
-----BEGIN PGP SIGNATURE-----
Version: GnuPG v1
Comment: Using GnuPG with Thunderbird - http://www.enigmail.net/

iQEcBAEBAgAGBQJTct57AAoJEGe6xJ1FYRpRqX0H/iy16jCiydYpT7XQAhGzrHCc
TBRJCBMDAorNS9CF2sl27tT/DOQ5y1Y8GBEdZSuThXS1m0LjM3MUTP4EK5xiYd5G
tpITSp1M/AfADmHIhQ/p7LqydmwF52gYG0LME/E8OLi1zvA88S08fmMA0GaFrZc+
DXyHunFqFrinn1KfYDShOFBrQWr01sDX/Zu5RBedQ8So5CwlQJCmT5ZHsLvR6jE8
hNrX9JwEXTqVQ/3lzByNDsXRPPVSytOEHKX3+tWVP1ou8ODfnHOyPVWEXiVuHacx
UU2DFcFWVZ31doxRl0dKqojSXpgFuG1AGavzX9vm2vBHOKyb36Kzx1VPz3lKZUQ=
=mxpF
-----END PGP SIGNATURE-----