[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
[DuMuX] Flyspray Activity for Task 184 (Naming conventions)
THIS IS AN AUTOMATED MESSAGE, DO NOT REPLY. The following task has a new comment added: FS#184 - Naming conventions User who did this - Thomas Fetzer (Thomas Fetzer) -- Here is one possibility to get rid of undesired variable names (expemlarily for the ElementIterator type and eIt or eEndIt as target variable names). I think these are the basic commands Zheng has used. # get all lines with the type of the variable, which is to be renamed grep -nri ElementIterator dumux/* test/* tutorial/* file1.txt # remove all lines containing desired variable name (see excel sheet), -v inverses the grep command grep -iv ElementIterator eIt file1.txt file2.txt grep -iv ElementIterator eEndIt file2.txt file3.txt # remove all typedefs grep -iv typedef typename GridView::template Codim0::Iterator ElementIterator file3.txt file4.txt -- More information can be found at the following URL: http://www.dumux.org/flyspray/index.php?do=detailstask_id=184#comment445 You are receiving this message because you have requested it from the Flyspray bugtracking system. If you did not expect this message or don't want to receive mails in future, you can change your notification settings at the URL shown above. ___ 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