Franco Le ven. 19 juil. 2019 à 10:04, Franco Milicchio <[email protected]> a écrit : > Is there a small documentation that I can read? For instance, I see there are > several namespaces like LinearAlgebra, PETSc/Epetra/Trilinos/CUDAWrappers, > LinearAlgebraDealii and so on, that might be restructured without > compromising old codes, adding some uniformity to classes. No unfortunately we don't have such things. However, most of the new effort on linear algebra is in the LinearAlgebra namespace.
> -------------------------------------------------------- > An error occurred in line <3781> of file > </Users/sensei/Downloads/src/dealii-9.1.0/include/deal.II/lac/affine_constraints.templates.h> > in function > void dealii::AffineConstraints<double>::distribute_local_to_global(const > FullMatrix<number> &, const Vector<number> &, const std::vector<size_type> &, > MatrixType &, VectorType &, bool, std::integral_constant<bool, true>) const > [number = double, MatrixType = dealii::BlockSparseMatrix<double>, VectorType > = dealii::BlockVector<double>] > The violated condition was: > static_cast<typename ::dealii::internal::argument_type<void( typename > std::common_type<decltype(global_matrix.m()), > decltype(global_vector.size())>::type)>::type>(global_matrix.m()) == > static_cast<typename ::dealii::internal::argument_type<void( typename > std::common_type<decltype(global_matrix.m()), > decltype(global_vector.size())>::type)>::type>(global_vector.size()) > Additional information: I don't know why this error is so unreadable but if you look at line 3781 in include/deal.II/lac/affine_constraints.templates.h you will see that the assert checks that the global matrix has as many rows has the global vector. To me it looks like the problem is that the global vector as a size of zero. If you forget about the static_cast in the error message just above you can actually see what happens: The violated condition was: global_matrix.m() == global_vector.size() Best, Bruno -- 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/CAGVt9eM%2Bg5HhuJOyQ3%3D-D__EbUyf7tZS5ZqyTu32rxVN1S-vXw%40mail.gmail.com.
