Good afternoon, Thank you for your answer. When using the function invert, I get sometimes the following running error:
*An error occurred in line <1849> of file </zhome/32/9/115503/dealii-candi/tmp/unpack/deal.II-v9.3.1/include/deal.II/lac/full_matrix.templates.h> in function* * void dealii::FullMatrix<number>::gauss_jordan() [with number = double]* *The violated condition was: * * info == 0* *Additional information: * * (none)* I am not sure what is happening, as the matrix is invertible. Sometimes after compiling few times and run, it works. Could you please help me understand what the problem could be? Thank you Regards, El miércoles, 26 de enero de 2022 a las 14:44:53 UTC+1, [email protected] escribió: > Hello, > > You can use the invert > <https://dealii.org/current/doxygen/deal.II/classFullMatrix.html#a940ea10266174ad366d10bd806177ace> > > function to get the inverse of the matrix and then you can use vmult with > the new matrix. > > Best, > > Bruno > > On Wednesday, January 26, 2022 at 8:19:46 AM UTC-5 [email protected] > wrote: > >> Good afternoon, >> >> I have a basic question. I can see that for sparse matrix the function >> SparseDirectUMFPACK.vmult() >> solves Ax=b. I am trying now to solve a system with full matrix. I can >> see that the function FullMatrix.vmult() does not apply the inverse as >> in the sparse analogy. What is the best way in Deal.ii to solve Ax=b with >> fullmatrix using a direct solver approach? >> >> Thank you >> >> -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum options, see https://groups.google.com/d/forum/dealii?hl=en --- You received this message because you are subscribed to the Google Groups "deal.II User Group" group. To unsubscribe from this group and stop receiving emails from it, send an email to [email protected]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/a31b89c2-b1fa-49e8-be3f-b2eca2b54ec0n%40googlegroups.com.
