Derek Gaston wrote:
> Thanks for the patch... I'll take a look at all of this shortly and
> commit it.
>
> As for the Jacobian in parallel... let's call not_implemented(). I
> would rather get an error than get completely different behavior and not
> know why.
Here's another small patch for epetra_vector.h. I had to change this in
order to print out the vector using print() - it segfaults if one uses
_map->... directly.
Norbert
Index: trilinos_epetra_vector.h
===================================================================
--- trilinos_epetra_vector.h (revision 3093)
+++ trilinos_epetra_vector.h (working copy)
@@ -743,7 +743,7 @@
{
libmesh_assert (this->initialized());
- return _map->MinMyGID();
+ return _vec->Map().MinMyGID();
}
@@ -754,7 +754,7 @@
{
libmesh_assert (this->initialized());
- return _map->MaxMyGID()+1;
+ return _vec->Map().MaxMyGID()+1;
}
-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Libmesh-devel mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-devel
