Hi Paul, Following up on our discussion on FD Jacobian through PetscNonlinearSolver.
It is not the coloring you need to generate, just the nonzero structure. >> (From that PETSc computes the coloring) so do as I suggest in my other >> email. >> > FYI, this information is contained in the libMesh::CouplingMatrix which > you can get from the libMesh::DofMap. > I tried to use CouplingMatrix in DofMap for building a "wrong" Jacobian just for the non-zero structure. Apparently the DofMap class contains a NULL pointer to the Coupling matrix. Do I need to call any function to initialize that? Thanks, Harshad On Fri, Jun 5, 2015 at 2:49 PM, Paul T. Bauman <[email protected]> wrote: > > > On Fri, Jun 5, 2015 at 2:42 PM, Harshad Sahasrabudhe <[email protected]> > wrote: > >> It is not the coloring you need to generate, just the nonzero >>> structure. (From that PETSc computes the coloring) so do as I suggest in my >>> other email. >> >> >> Thanks, I'll try that out. >> > > FYI, this information is contained in the libMesh::CouplingMatrix which > you can get from the libMesh::DofMap. > ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
