Re: [DuMuX] Eigenvalues of the jacobian matrix
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: dumux-boun...@listserv.uni-stuttgart.de [mailto:dumux-boun...@listserv.uni-stuttgart.de] Im Auftrag von Christoph Grüninger Gesendet: Donnerstag, 3. April 2014 07:17 An: dumux@listserv.uni-stuttgart.de 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 Dumux@listserv.uni-stuttgart.de https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
[DuMuX] Eigenvalues of the jacobian matrix
Hello Everyone, To investigate some instabilities in my model I would like to have a look at the eigenvalues of the jacobian matrix. Is there an example how this can be done? Thanks Georg -- Deutsches Zentrum für Luft- und Raumfahrt e.V. (DLR) Institut für Technische Thermodynamik | Computergestützte Elektrochemie | Pfaffenwaldring 38-40 | 70569 Stuttgart Dipl.-Ing. Georg Futter | Doktorand Telefon 0711/6862-8135 | georg.fut...@dlr.demailto:georg.fut...@dlr.de www.DLR.dehttp://www.dlr.de/ ___ Dumux mailing list Dumux@listserv.uni-stuttgart.de https://listserv.uni-stuttgart.de/mailman/listinfo/dumux
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 * signature.asc Description: OpenPGP digital signature ___ Dumux mailing list Dumux@listserv.uni-stuttgart.de https://listserv.uni-stuttgart.de/mailman/listinfo/dumux