Hi, I wish to have a rigid boundary condition. I am trying to do so with add_constrain_row. The idea would be to tie the Dofs to each other in the constrain funcion, befor calling the solver. Not sure it is the best idea, though. Each processor adds a constrain row to the system regarding the elements on its domain.
I see this error: [2]PETSC ERROR: --------------------- Error Message -------------------------------------------------------------- [2]PETSC ERROR: Argument out of range [2]PETSC ERROR: New nonzero at (72,1051) caused a malloc Use MatSetOption(A, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE) to turn off this check [2]PETSC ERROR: See http://www.mcs.anl.gov/petsc/documentation/faq.html for trouble shooting. ... Should I consider adding MatSetOption? Any advice? Thanks, Renato _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
