> > How far non-local and how sparse are you talking? If i and j are on > neighboring elements we have some options for automatically expanding > the Jacobian that far. If i and j might be halfway across the mesh > then support is more limited and potentially much more expensive.
I basically have a dense Jacobian in some parts of the mesh. Though I can first try using a Jacobian which spans the element and neighbors. I already tried a diagonal Jacobian though, and it doesn't work in my problem. What are the options to expand the Jacobian? ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users