Hi, thank You for quick answer. Could the nonsymmetry be introduced by the fact, that I use the fully distributed setting (in the way of Step-40)? Could the method create_mass_matrix be used for the fully distributed case?
Thank You Marek C Dne středa 9. ledna 2019 0:36:11 UTC+1 Daniel Arndt napsal(a): > > Marek, > > if I were you, I would first check that the assembled matrix (without > constraints) is correct. > For a suitably small mesh you can compute the expected matrix by hand if > necessary. > Alternatively, there is MatrixCreator::create_mass_matrix ( > https://www.dealii.org/9.0.0/doxygen/deal.II/namespaceMatrixCreator.html#a298a5ae54ecf4709eca32d9b81405274 > ). > > If that works correctly, I would try to incorporate the constraints (check > correctness again if possible) > and use a direct (serial) solver with serial linear algebra classes. > > Best, > Daniel > -- 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.
