> No unfortunately we don't have such things. However, most of the new
> effort on linear algebra is in the LinearAlgebra namespace.
>
>
Thanks, Bruno, I've started looking into that namespace. There are several
options to make it work in a coherent way, I'll ask for comments on these
solutions as I figured all classes out (and how to use MPI seamlessly, at
first I am thinking about tag dispatching to avoid breaking legacy codes).
>
> 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()
>
>
I am lost in this, I have made some stupid error here but I don't see how.
I am posting the code here:
https://gist.github.com/fmilicchio/2dac430d7c071f98b4f66d79d2101b04
As far as I can tell, I've initialized every sparsity pattern, matrix, and
vector as they were with Trilinos, but I am surely seeing the code I
intended to write, not the code I have actually written.
Do you see any stupid error at first sight?
Thank you!
Franco
--
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/ce549f8c-9a44-462f-a66f-fcfab2794b3e%40googlegroups.com.
