Changeset [d022dba2e105ab10b58767087497466aa2c2f80b] by Jed Brown

November 21st, 2008 @ 08:02 PM

Solve projection problems with two elements. [#1 status:resolved]

Two main weaknesses:

  • Uses -snes_mf for projection, we should also assemble a matrix and be able to apply the Jacobian matrix-free.

  • Discretization error is not yet checked, but evaluation at quadrature points demonstrates spectral convergence so the projection should as well.

Signed-off-by: Jed Brown jed@59A2.org http://github.com/jedbrown/dohp/...

Committed by Jed Brown

  • M CMakeLists.txt
  • M include/dohpfs.h
  • M include/dohpjacobi.h
  • M include/dohpmesh.h
  • M include/dohptype.h
  • M include/private/fsimpl.h
  • M include/private/jacimpl.h
  • M src/fs/impls/cont/cont.c
  • M src/fs/impls/cont/cont.h
  • M src/fs/interface/fs.c
  • M src/fs/mesh/interface/mesh.c
  • M src/fs/tests/CMakeLists.txt
  • M src/fs/tests/ex1.c
  • M src/jacobi/impls/tensor/ruletopo.c
  • M src/jacobi/interface/jacobi.c
  • M src/jacobi/tests/CMakeLists.txt
New-ticket Create new ticket

Create your profile

Help contribute to this project by taking a few moments to create your personal profile. Create your profile »

An implementation of the ``dual order hp'' version of the finite element method. This project targets parallel domain-decomposition methods for strongly coupled nonlinear problems with PDE constraints.