> I think DenseMatrix::condense() will work to zero out values, but if you > want to put some nonzero value in, e.g., the diagonal of the assembled > matrix, it seems like you need to do that after an initial assembly, > don't you?
Not sure if this is the same concern as the previous - about adding the diagonal multiple times, but if it is that is not a problem because while you do add the diagonal multiple times you also add the RHS multiple times. So for a DOF shared by say 5 elements rather than [0 ... 0 1 0 ....0] dof = val you get instead [0 ... 0 5 0 ....0] dof = 5*val -Ben ------------------------------------------------------------------------------ Achieve unprecedented app performance and reliability What every C/C++ and Fortran developer should know. Learn how Intel has extended the reach of its next-generation tools to help boost performance applications - inlcuding clusters. http://p.sf.net/sfu/intel-dev2devmay _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
