Hi Christoph,

Thank you for your quick answer. I was wondering whether I could use the 
functions in eigenvalues.hh. Is this not a good idea?

Regards
Georg

-----Ursprüngliche Nachricht-----
Von: [email protected] 
[mailto:[email protected]] Im Auftrag von Christoph 
Grüninger
Gesendet: Donnerstag, 3. April 2014 07:17
An: [email protected]
Betreff: Re: [DuMuX] Eigenvalues of the jacobian matrix

Hi Georg,
as far as I know, you can not get the eigenvalues out of the box.
Dune-ISTL has some eigenvalue computation support.
I'd advice you to hack the DuMuX assembler and print everything into a 
Matlab/Octave file. Using their capabilities it's easy to obtain the 
eigenvalues.
- you can get the matrix for example in
dumux/implicit/common/implicitassembler.hh at the end of the assemble() method
- use dune-istl's Dune::writeMatrixToMatlab

Bye
Christoph

--
... as Sir Cyril Hinshelwood has observed ...fluid dynamicists were divided 
into hydraulic engineers who observed things that could not be explained and 
mathematicians who explained things
that could not be observed.                 -- James Lighthill
*********************************************
CMWR 2014: 10th - 13th June 2014 in Stuttgart
         Please visit www.cmwr14.de
*********************************************

_______________________________________________
Dumux mailing list
[email protected]
https://listserv.uni-stuttgart.de/mailman/listinfo/dumux

Reply via email to