Dear all, The method NumericVector::add(const T) seems to be completely unused; at least the doxygen documentation states no reference to that method. Actually, I'm quite sure the the implementation PetscVector::add(const T) would not work for parallel vectors, because VecGetArray() supplies an array that uses local indices, whereas VecSetValues() requires a global index. Confusing enough of course. (There is also a VecSetValueLocal() function.) In my patch that I'm preparing, I will also fix that.
Best Regards, Tim -- Dr. Tim Kroeger [email protected] Phone +49-421-218-7710 [email protected] Fax +49-421-218-4236 Fraunhofer MEVIS, Institute for Medical Image Computing Universitaetsallee 29, 28359 Bremen, Germany ------------------------------------------------------------------------------ This SF.net email is sponsored by: SourcForge Community SourceForge wants to tell your story. http://p.sf.net/sfu/sf-spreadtheword _______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
