[petsc-users] DMPlex with spring elements

[Miguel Angel Salazar de Troya]

Hi all

I was wondering if it could be possible to build a model similar to the
example snes/ex12.c, but with spring elements (for elasticity) instead of
simplicial elements. Spring elements in a grid, therefore each element
would have two nodes and each node two components. There would be more
differences, because instead of calling the functions f0,f1,g0,g1,g2 and g3
to build the residual and the jacobian, I would call a routine that would
build the residual vector and the jacobian matrix directly. I would not
have shape functions whatsoever. My problem is discrete, I don't have a PDE
and my equations are algebraic. What is the best way in petsc to solve this
problem? Is there any example that I can follow? Thanks in advance

Miguel



-- 
*Miguel Angel Salazar de Troya*
Graduate Research Assistant
Department of Mechanical Science and Engineering
University of Illinois at Urbana-Champaign
(217) 550-2360
salaza11 at illinois.edu
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/petsc-users/attachments/20140923/270b8fc5/attachment.html>