Hi all, Just refreshing this one, because I am sort of stuck in inserting a "rigid" BCs. It seems simple, but I cannot make it work.
Some help would be handy... Maybe adding an extra 'virtual' element would be the best idea? How can I do that? Thanks upfront. Renato On Mon, Feb 18, 2019 at 5:52 PM Renato Poli <rebp...@gmail.com> wrote: > 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
