Hello, Could you print the value of info (you can use gdb to do that). That would help us understand what's the problem.
Best, Bruno Le mer. 2 févr. 2022 à 07:21, Joss G. <[email protected]> a écrit : > > 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 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 a topic in the Google > Groups "deal.II User Group" group. > To unsubscribe from this topic, visit > https://groups.google.com/d/topic/dealii/uq1J51MCS80/unsubscribe. > To unsubscribe from this group and all its topics, 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. -- 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/CAGVt9eNBQ6ye-su59XGb_Oyyy0mNLonZv6ehAhN%3Dk7K6kZgCpQ%40mail.gmail.com.
