Hi, I need to assemble L2 projection manually, byhand, not using "project" methods. As the resulting mass matrix is symmetric and positive definite (should be), I use CG+SSOR as a solver. However I get "loss of precision error" - I perform fully distributed computation using TrilinosWrappers for vectors, matrices, solvers. I guess, that my matrix is not symmetric and(or) positive definite. Could the problem originate from the improperly configurated constraint matrix, which would result into bad matrix? I have no boundary condition on the system. How to configure constraint matrix so that the properties of matrix to be solved by iterative solver would be retained?
Thank You MareK C -- 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]. For more options, visit https://groups.google.com/d/optout.
