← Back to branch summary

~njansson/dolfin/hpc

~njansson/dolfin/hpc

« back to all changes in this revision

Viewing changes to src/demo/fem/convergence/Poisson2D_3.h

Committer: Anders Logg
Date: 2007-01-10 09:04:44 UTC
mfrom: (1689.1.221 trunk)
Revision ID: logg@simula.no-20070110090444-ecyux3n1qnei4i8h

Remove old head

files added:
doc/manual/chapters/design.tex

doc/manual/dia

doc/manual/dia/class-diagram-la.dia

doc/manual/eps/class-diagram-la.eps

dolfin.pc.in

m4/ax_boost_base.m4

m4/ax_compiler_vendor.m4

m4/libtool.m4

m4/libxml.m4

numpy.conf

petsc.conf

slepc.conf

src/bench

src/bench/Makefile.am

src/bench/Makefile.in

src/bench/README

src/bench/fem

src/bench/fem/Makefile.am

src/bench/fem/Makefile.in

src/bench/fem/assembly

src/bench/fem/assembly/AdvectionOperator_3D_3_FER.form

src/bench/fem/assembly/AdvectionOperator_3D_3_FER.h

src/bench/fem/assembly/AdvectionOperator_3D_3_FFC.form

src/bench/fem/assembly/AdvectionOperator_3D_3_FFC.h

src/bench/fem/assembly/Makefile

src/bench/fem/assembly/README

src/bench/fem/assembly/bench.log

src/bench/fem/assembly/bench.py

src/bench/fem/assembly/main.cpp

src/bench/fem/convergence

src/bench/fem/convergence/Makefile

src/bench/fem/convergence/Poisson2D_1.form

src/bench/fem/convergence/Poisson2D_1.h

src/bench/fem/convergence/Poisson2D_2.form

src/bench/fem/convergence/Poisson2D_2.h

src/bench/fem/convergence/Poisson2D_3.form

src/bench/fem/convergence/Poisson2D_3.h

src/bench/fem/convergence/Poisson2D_4.form

src/bench/fem/convergence/Poisson2D_4.h

src/bench/fem/convergence/Poisson2D_5.form

src/bench/fem/convergence/Poisson2D_5.h

src/bench/fem/convergence/Poisson3D_1.form

src/bench/fem/convergence/Poisson3D_1.h

src/bench/fem/convergence/Poisson3D_2.form

src/bench/fem/convergence/Poisson3D_2.h

src/bench/fem/convergence/Poisson3D_3.form

src/bench/fem/convergence/Poisson3D_3.h

src/bench/fem/convergence/Poisson3D_4.form

src/bench/fem/convergence/Poisson3D_4.h

src/bench/fem/convergence/Poisson3D_5.form

src/bench/fem/convergence/Poisson3D_5.h

src/bench/fem/convergence/README

src/bench/fem/convergence/bench.log

src/bench/fem/convergence/bench.py

src/bench/fem/convergence/main.cpp

src/bench/la

src/bench/la/Makefile.am

src/bench/la/Makefile.in

src/bench/la/bench.log

src/bench/la/bench.py

src/bench/la/sparse-matrix

src/bench/la/sparse-matrix/Makefile

src/bench/la/sparse-matrix/Poisson.form

src/bench/la/sparse-matrix/Poisson.h

src/bench/la/sparse-matrix/VectorPoisson.form

src/bench/la/sparse-matrix/VectorPoisson.h

src/bench/la/sparse-matrix/main.cpp

src/bench/la/vector

src/bench/la/vector/Makefile

src/bench/la/vector/main.cpp

src/bench/mesh

src/bench/mesh/Makefile

src/bench/mesh/bench.log

src/bench/mesh/bench.py

src/bench/mesh/main.cpp

src/bench/ode

src/bench/ode/Makefile.am

src/bench/ode/Makefile.in

src/bench/ode/reaction

src/bench/ode/reaction/Makefile

src/bench/ode/reaction/README

src/bench/ode/reaction/bench-domain.log

src/bench/ode/reaction/bench-domain.py

src/bench/ode/reaction/bench-overhead.log

src/bench/ode/reaction/bench-overhead.py

src/bench/ode/reaction/bench-plot.py

src/bench/ode/reaction/bench-tol.log

src/bench/ode/reaction/bench-tol.py

src/bench/ode/reaction/bench.log

src/bench/ode/reaction/bench.py

src/bench/ode/reaction/benchutil.py

src/bench/ode/reaction/checkerror.py

src/bench/ode/reaction/main.cpp

src/bench/ode/reaction/parameters-bench.xml

src/bench/ode/reaction/parameters-plot.xml

src/bench/ode/reaction/parameters-reference.xml

src/bench/ode/reaction/reference_1000.data

src/bench/ode/reaction/reference_16000.data

src/bench/ode/reaction/reference_2000.data

src/bench/ode/reaction/reference_4000.data

src/bench/ode/reaction/reference_8000.data

src/bench/ode/wave

src/bench/ode/wave/Makefile

src/bench/ode/wave/bench.log

src/bench/ode/wave/bench.py

src/bench/ode/wave/main.cpp

src/bench/ode/wave/parameters-bench.xml

src/bench/ode/wave/slit.xml.gz

src/demo/mesh/plot.py

src/demo/meshfunction

src/demo/meshfunction/Makefile

src/demo/meshfunction/main.cpp

src/demo/meshfunction/mesh2D.xml.gz

src/demo/meshfunction/meshfunction.xml

src/demo/nls/cahn-hilliard

src/demo/nls/cahn-hilliard/CahnHilliard2D.form

src/demo/nls/cahn-hilliard/CahnHilliard2D.h

src/demo/nls/cahn-hilliard/CahnHilliard3D.form

src/demo/nls/cahn-hilliard/CahnHilliard3D.h

src/demo/nls/cahn-hilliard/Makefile

src/demo/nls/cahn-hilliard/main.cpp

src/demo/nls/cahn-hilliard/plot.py

src/demo/nls/nonlinearpoisson/plot.py

src/demo/ode/aliev-panfilov

src/demo/ode/aliev-panfilov/Makefile

src/demo/ode/aliev-panfilov/main.cpp

src/demo/ode/aliev-panfilov/plot.py

src/demo/ode/complex/plot.py

src/demo/ode/courtemanche

src/demo/ode/courtemanche/Makefile

src/demo/ode/courtemanche/main.cpp

src/demo/ode/courtemanche/plot.py

src/demo/ode/harmonic

src/demo/ode/harmonic/Makefile

src/demo/ode/harmonic/main.cpp

src/demo/ode/harmonic/plot.py

src/demo/ode/homotopy/simple/plot.py

src/demo/ode/lorenz/plot.py

src/demo/ode/stiff/plot.py

src/demo/pde/bidomain

src/demo/pde/bidomain/Bidomain.form

src/demo/pde/bidomain/Bidomain.h

src/demo/pde/bidomain/Makefile

src/demo/pde/bidomain/main.cpp

src/demo/pde/convection-diffusion/plot.py

src/demo/pde/dg

src/demo/pde/dg/Makefile

src/demo/pde/dg/Poisson.form

src/demo/pde/dg/Poisson.h

src/demo/pde/dg/Projection.form

src/demo/pde/dg/Projection.h

src/demo/pde/dg/main.cpp

src/demo/pde/dg/plot.py

src/demo/pde/elasticity/plot.py

src/demo/pde/functional

src/demo/pde/functional/EnergyNorm.form

src/demo/pde/functional/EnergyNorm.h

src/demo/pde/functional/Makefile

src/demo/pde/functional/main.cpp

src/demo/pde/nonlinear/nonlinear-poisson/plot.py

src/demo/pde/poisson/plot.py

src/demo/pde/stokes/stabilized/plot.py

src/demo/pde/stokes/taylor-hood/plot.py

src/demo/scripting/pydolfin/output.py

src/demo/scripting/pydolfin/projection

src/demo/scripting/pydolfin/projection/projectiondemo.py

src/demo/scripting/pydolfin/solvers/elasticity-tpde

src/demo/scripting/pydolfin/solvers/elasticity-tpde/composite

src/demo/scripting/pydolfin/solvers/elasticity-tpde/composite/Coefficients.h

src/demo/scripting/pydolfin/solvers/elasticity-tpde/composite/Coefficients.i

src/demo/scripting/pydolfin/solvers/elasticity-tpde/composite/Elasticity.form

src/demo/scripting/pydolfin/solvers/elasticity-tpde/composite/elasticityform.py

src/demo/scripting/pydolfin/solvers/elasticity-tpde/composite/elasticitypde.py

src/demo/scripting/pydolfin/solvers/elasticity-tpde/composite/elasticitysolver.py

src/demo/scripting/pydolfin/solvers/elasticity-tpde/composite/plot.py

src/demo/scripting/pydolfin/solvers/elasticity-tpde/composite/tetmesh-4.xml.gz

src/demo/scripting/pydolfin/solvers/elasticity-tpde/manual

src/demo/scripting/pydolfin/solvers/elasticity-tpde/manual/Elasticity.form

src/demo/scripting/pydolfin/solvers/elasticity-tpde/manual/elasticitypde.py

src/demo/scripting/pydolfin/solvers/elasticity-tpde/manual/elasticitysolver.py

src/demo/scripting/pydolfin/solvers/elasticity-tpde/manual/plot.py

src/demo/scripting/pydolfin/solvers/elasticity-tpde/manual/tetmesh-4.xml.gz

src/demo/scripting/pydolfin/solvers/elasticity-updated

src/demo/scripting/pydolfin/solvers/elasticity-updated/diff

src/demo/scripting/pydolfin/solvers/elasticity-updated/diff/ElasticityLoad.form

src/demo/scripting/pydolfin/solvers/elasticity-updated/diff/ElasticityUpdatedEq.form

src/demo/scripting/pydolfin/solvers/elasticity-updated/diff/ElasticityUpdatedStress.form

src/demo/scripting/pydolfin/solvers/elasticity-updated/diff/diffpde.py

src/demo/scripting/pydolfin/solvers/elasticity-updated/diff/geometry.py

src/demo/scripting/pydolfin/solvers/elasticity-updated/diff/plot.py

src/demo/scripting/pydolfin/solvers/elasticity-updated/diff/problem.py

src/demo/scripting/pydolfin/solvers/elasticity-updated/diff/tetmesh-1c.xml.gz

src/demo/scripting/pydolfin/solvers/elasticity-updated/diff/transform.py

src/demo/scripting/pydolfin/solvers/elasticity-updated/direct

src/demo/scripting/pydolfin/solvers/elasticity-updated/direct/ElasticityDirect.form

src/demo/scripting/pydolfin/solvers/elasticity-updated/direct/ElasticityLoad.form

src/demo/scripting/pydolfin/solvers/elasticity-updated/direct/directpde.py

src/demo/scripting/pydolfin/solvers/elasticity-updated/direct/geometry.py

src/demo/scripting/pydolfin/solvers/elasticity-updated/direct/plot.py

src/demo/scripting/pydolfin/solvers/elasticity-updated/direct/problem.py

src/demo/scripting/pydolfin/solvers/elasticity-updated/direct/tetmesh-1c.xml.gz

src/demo/scripting/pydolfin/solvers/elasticity-updated/direct/transform.py

src/demo/scripting/pydolfin/solvers/elasticity-updated/plast

src/demo/scripting/pydolfin/solvers/elasticity-updated/plast/Coefficients.h

src/demo/scripting/pydolfin/solvers/elasticity-updated/plast/Coefficients.i

src/demo/scripting/pydolfin/solvers/elasticity-updated/plast/ElasticityLoad.form

src/demo/scripting/pydolfin/solvers/elasticity-updated/plast/ElasticityUpdatedEq.form

src/demo/scripting/pydolfin/solvers/elasticity-updated/plast/ElasticityUpdatedStress.form

src/demo/scripting/pydolfin/solvers/elasticity-updated/plast/diffpde.py

src/demo/scripting/pydolfin/solvers/elasticity-updated/plast/geometry.py

src/demo/scripting/pydolfin/solvers/elasticity-updated/plast/plot.py

src/demo/scripting/pydolfin/solvers/elasticity-updated/plast/problem.py

src/demo/scripting/pydolfin/solvers/elasticity-updated/plast/tetmesh-1c.xml.gz

src/demo/scripting/pydolfin/solvers/elasticity-updated/plast/transform.py

src/demo/scripting/pydolfin/solvers/poisson/Poisson.form

src/demo/solvers/convdiff/plot.py

src/demo/solvers/elasticity-updated/plot.py

src/demo/solvers/elasticity-updated/tetmesh-1c.xml.gz

src/demo/solvers/elasticity/plot.py

src/demo/solvers/heat/plot.py

src/demo/solvers/plasticity

src/demo/solvers/plasticity/Makefile

src/demo/solvers/plasticity/main.cpp

src/demo/solvers/plasticity/plot.py

src/demo/solvers/poisson/plot.py

src/demo/solvers/stokes/plot.py

src/kernel/fem/DofMap.cpp

src/kernel/fem/dolfin/DofMap.h

src/kernel/io/NewXMLMesh.cpp

src/kernel/io/XMLMeshFunction.cpp

src/kernel/io/dolfin/NewXMLMesh.h

src/kernel/io/dolfin/XMLMeshFunction.h

src/kernel/la/PETScEigenvalueSolver.cpp

src/kernel/la/PETScKrylovMatrix.cpp

src/kernel/la/PETScKrylovSolver.cpp

src/kernel/la/PETScLUSolver.cpp

src/kernel/la/PETScMatrix.cpp

src/kernel/la/PETScPreconditioner.cpp

src/kernel/la/PETScVector.cpp

src/kernel/la/dolfin/DenseMatrix.h

src/kernel/la/dolfin/DenseVector.h

src/kernel/la/dolfin/GenericMatrix.h

src/kernel/la/dolfin/GenericVector.h

src/kernel/la/dolfin/KrylovMethod.h

src/kernel/la/dolfin/LUSolver.h

src/kernel/la/dolfin/PETScEigenvalueSolver.h

src/kernel/la/dolfin/PETScKrylovMatrix.h

src/kernel/la/dolfin/PETScKrylovSolver.h

src/kernel/la/dolfin/PETScLUSolver.h

src/kernel/la/dolfin/PETScLinearSolver.h

src/kernel/la/dolfin/PETScMatrix.h

src/kernel/la/dolfin/PETScPreconditioner.h

src/kernel/la/dolfin/PETScVector.h

src/kernel/la/dolfin/SparseMatrix.h

src/kernel/la/dolfin/SparseVector.h

src/kernel/la/dolfin/uBlasDenseMatrix.h

src/kernel/la/dolfin/uBlasDummyPreconditioner.h

src/kernel/la/dolfin/uBlasILUPreconditioner.h

src/kernel/la/dolfin/uBlasKrylovMatrix.h

src/kernel/la/dolfin/uBlasKrylovSolver.h

src/kernel/la/dolfin/uBlasLUSolver.h

src/kernel/la/dolfin/uBlasLinearSolver.h

src/kernel/la/dolfin/uBlasMatrix.h

src/kernel/la/dolfin/uBlasPreconditioner.h

src/kernel/la/dolfin/uBlasSparseMatrix.h

src/kernel/la/dolfin/uBlasVector.h

src/kernel/la/dolfin/ublas.h

src/kernel/la/uBlasDummyPreconditioner.cpp

src/kernel/la/uBlasILUPreconditioner.cpp

src/kernel/la/uBlasKrylovMatrix.cpp

src/kernel/la/uBlasKrylovSolver.cpp

src/kernel/la/uBlasLUSolver.cpp

src/kernel/la/uBlasVector.cpp

src/kernel/mesh/BoundaryComputation.cpp

src/kernel/mesh/BoundaryMesh.cpp

src/kernel/mesh/CellType.cpp

src/kernel/mesh/Facet.cpp

src/kernel/mesh/GeometricPredicates.cpp

src/kernel/mesh/Interval.cpp

src/kernel/mesh/LocalMeshCoarsening.cpp

src/kernel/mesh/LocalMeshRefinement.cpp

src/kernel/mesh/MeshConnectivity.cpp

src/kernel/mesh/MeshEditor.cpp

src/kernel/mesh/MeshEntity.cpp

src/kernel/mesh/MeshEntityIterator.cpp

src/kernel/mesh/MeshGeometry.cpp

src/kernel/mesh/MeshTopology.cpp

src/kernel/mesh/TopologyComputation.cpp

src/kernel/mesh/UniformMeshRefinement.cpp

src/kernel/mesh/dolfin/BoundaryComputation.h

src/kernel/mesh/dolfin/BoundaryMesh.h

src/kernel/mesh/dolfin/CellType.h

src/kernel/mesh/dolfin/Facet.h

src/kernel/mesh/dolfin/GeometricPredicates.h

src/kernel/mesh/dolfin/Interval.h

src/kernel/mesh/dolfin/LocalMeshCoarsening.h

src/kernel/mesh/dolfin/LocalMeshRefinement.h

src/kernel/mesh/dolfin/MeshConnectivity.h

src/kernel/mesh/dolfin/MeshEditor.h

src/kernel/mesh/dolfin/MeshEntity.h

src/kernel/mesh/dolfin/MeshEntityIterator.h

src/kernel/mesh/dolfin/MeshFunction.h

src/kernel/mesh/dolfin/MeshGeometry.h

src/kernel/mesh/dolfin/MeshTopology.h

src/kernel/mesh/dolfin/TopologyComputation.h

src/kernel/mesh/dolfin/UniformMeshRefinement.h

src/kernel/mesh/predicates_init.h

src/kernel/mesh/rounding.h

src/kernel/mf

src/kernel/mf/Makefile.am

src/kernel/mf/Makefile.in

src/kernel/mf/MatrixFactory.cpp

src/kernel/mf/dolfin

src/kernel/mf/dolfin/ConvectionMatrix.h

src/kernel/mf/dolfin/LoadVector.h

src/kernel/mf/dolfin/Makefile.am

src/kernel/mf/dolfin/Makefile.in

src/kernel/mf/dolfin/MassMatrix.h

src/kernel/mf/dolfin/MatrixFactory.h

src/kernel/mf/dolfin/StiffnessMatrix.h

src/kernel/mf/dolfin/dolfin_mf.h

src/kernel/mf/ffc-forms

src/kernel/mf/ffc-forms/ConvectionMatrix2D.form

src/kernel/mf/ffc-forms/ConvectionMatrix2D.h

src/kernel/mf/ffc-forms/ConvectionMatrix3D.form

src/kernel/mf/ffc-forms/ConvectionMatrix3D.h

src/kernel/mf/ffc-forms/LoadVector2D.form

src/kernel/mf/ffc-forms/LoadVector2D.h

src/kernel/mf/ffc-forms/LoadVector3D.form

src/kernel/mf/ffc-forms/LoadVector3D.h

src/kernel/mf/ffc-forms/MassMatrix2D.form

src/kernel/mf/ffc-forms/MassMatrix2D.h

src/kernel/mf/ffc-forms/MassMatrix3D.form

src/kernel/mf/ffc-forms/MassMatrix3D.h

src/kernel/mf/ffc-forms/StiffnessMatrix2D.form

src/kernel/mf/ffc-forms/StiffnessMatrix2D.h

src/kernel/mf/ffc-forms/StiffnessMatrix3D.form

src/kernel/mf/ffc-forms/StiffnessMatrix3D.h

src/kernel/pde/TimeDependentPDE.cpp

src/kernel/pde/dolfin/TimeDependentPDE.h

src/modules/navierstokes/dolfin/NSEDualContinuity2D.form

src/modules/navierstokes/dolfin/NSEDualContinuity2D.h

src/modules/navierstokes/dolfin/NSEDualContinuity3D.form

src/modules/navierstokes/dolfin/NSEDualContinuity3D.h

src/modules/navierstokes/dolfin/NSEDualMomentum2D.form

src/modules/navierstokes/dolfin/NSEDualMomentum2D.h

src/modules/navierstokes/dolfin/NSEDualMomentum3D.form

src/modules/plasticity

src/modules/plasticity/DruckerPrager.cpp

src/modules/plasticity/Makefile.am

src/modules/plasticity/Makefile.in

src/modules/plasticity/PlasticityModel.cpp

src/modules/plasticity/PlasticityProblem.cpp

src/modules/plasticity/PlasticitySolver.cpp

src/modules/plasticity/ReturnMapping.cpp

src/modules/plasticity/VonMises.cpp

src/modules/plasticity/dolfin

src/modules/plasticity/dolfin/DruckerPrager.h

src/modules/plasticity/dolfin/Makefile.am

src/modules/plasticity/dolfin/Makefile.in

src/modules/plasticity/dolfin/Output2D.form

src/modules/plasticity/dolfin/Output2D.h

src/modules/plasticity/dolfin/Output3D.form

src/modules/plasticity/dolfin/Output3D.h

src/modules/plasticity/dolfin/Plas2D.form

src/modules/plasticity/dolfin/Plas2D.h

src/modules/plasticity/dolfin/Plas3D.form

src/modules/plasticity/dolfin/Plas3D.h

src/modules/plasticity/dolfin/PlasticityModel.h

src/modules/plasticity/dolfin/PlasticityProblem.h

src/modules/plasticity/dolfin/PlasticitySolver.h

src/modules/plasticity/dolfin/ReturnMapping.h

src/modules/plasticity/dolfin/Strain2D.form

src/modules/plasticity/dolfin/Strain2D.h

src/modules/plasticity/dolfin/Strain3D.form

src/modules/plasticity/dolfin/Strain3D.h

src/modules/plasticity/dolfin/Tangent2D.form

src/modules/plasticity/dolfin/Tangent2D.h

src/modules/plasticity/dolfin/Tangent3D.form

src/modules/plasticity/dolfin/Tangent3D.h

src/modules/plasticity/dolfin/VonMises.h

src/modules/plasticity/dolfin/dolfin_plasticity.h

src/modules/plasticity/dolfin/ep_strain2D.form

src/modules/plasticity/dolfin/ep_strain2D.h

src/modules/plasticity/dolfin/ep_strain3D.form

src/modules/plasticity/dolfin/ep_strain3D.h

src/modules/plasticity/dolfin/p_strain2D.form

src/modules/plasticity/dolfin/p_strain2D.h

src/modules/plasticity/dolfin/p_strain3D.form

src/modules/plasticity/dolfin/p_strain3D.h

src/pydolfin/assemble.py

src/pydolfin/dolfin_fem.i

src/pydolfin/dolfin_fem_post.i

src/pydolfin/dolfin_glue.cpp

src/pydolfin/dolfin_glue.h

src/pydolfin/dolfin_la.i

src/pydolfin/dolfin_mesh.i

src/pydolfin/dolfin_mesh_post.i

src/pydolfin/dolfin_mesh_pre.i

src/pydolfin/dolfin_ode.i

src/pydolfin/projection.py

src/test/la

src/test/la/__init__.py

src/test/la/algebraUBlas.py

src/test/la/createUBlas.py

src/test/la/solverUBlas.py

src/test/la/test.py

src/test/mesh

src/test/mesh/__init__.py

src/test/mesh/test.py

src/test/passembly

src/test/passembly/Makefile

src/test/passembly/main.cpp

src/test/test.py

src/test/tmp

src/test/tmp/L2Norm.form

src/test/tmp/L2Norm.h

src/test/tmp/Makefile

src/test/tmp/main.cpp

src/test/tmp/newmesh.xml.gz

src/test/tmp/oldmesh.xml.gz

src/utils/convert/test_gmsh.msh

src/utils/convert/test_old_xml.xml

src/utils/convert/test_tetgen.mesh

src/utils/emacs

src/utils/emacs/macros

src/utils/swig

src/utils/swig/Makefile.am

src/utils/swig/Makefile.in

src/utils/swig/dolfin-swig

src/utils/vim

src/utils/vim/macros

src/utils/xml/convertall

files removed:
doc/misc/cvs

doc/misc/releasing

scripts/preparedist

src/config

src/config/Makefile.am

src/config/Makefile.in

src/config/config.sh

src/config/dolfin-config.template

src/config/dolfin-swig

src/config/dolfin-swig.template

src/config/swig.sh

src/demo/fem/convergence

src/demo/fem/convergence/Makefile

src/demo/fem/convergence/Poisson2D_1.form

src/demo/fem/convergence/Poisson2D_1.h

src/demo/fem/convergence/Poisson2D_2.form

src/demo/fem/convergence/Poisson2D_2.h

src/demo/fem/convergence/Poisson2D_3.form

src/demo/fem/convergence/Poisson2D_3.h

src/demo/fem/convergence/Poisson2D_4.form

src/demo/fem/convergence/Poisson2D_4.h

src/demo/fem/convergence/Poisson2D_5.form

src/demo/fem/convergence/Poisson2D_5.h

src/demo/fem/convergence/Poisson3D_1.form

src/demo/fem/convergence/Poisson3D_1.h

src/demo/fem/convergence/Poisson3D_2.form

src/demo/fem/convergence/Poisson3D_2.h

src/demo/fem/convergence/Poisson3D_3.form

src/demo/fem/convergence/Poisson3D_3.h

src/demo/fem/convergence/Poisson3D_4.form

src/demo/fem/convergence/Poisson3D_4.h

src/demo/fem/convergence/Poisson3D_5.form

src/demo/fem/convergence/Poisson3D_5.h

src/demo/fem/convergence/README

src/demo/fem/convergence/main.cpp

src/demo/ode/bistable

src/demo/ode/bistable/Makefile

src/demo/ode/bistable/README

src/demo/ode/bistable/main.cpp

src/demo/ode/bistable/mesh.xml.gz

src/demo/ode/homotopy/economy/plotsolution.m

src/demo/ode/homotopy/simple/plotsolution.m

src/demo/ode/lorenz/plotlaps.m

src/demo/ode/mechanical

src/demo/ode/mechanical/Makefile

src/demo/ode/mechanical/main.cpp

src/demo/ode/modeling

src/demo/ode/modeling/Makefile.am

src/demo/ode/modeling/Makefile.in

src/demo/ode/modeling/lattice

src/demo/ode/modeling/lattice/Makefile

src/demo/ode/modeling/lattice/main.cpp

src/demo/ode/modeling/lattice/plotsolution.m

src/demo/ode/modeling/simple

src/demo/ode/modeling/simple/Makefile

src/demo/ode/modeling/simple/main.cpp

src/demo/ode/modeling/simple/plotsolution.m

src/demo/ode/reaction

src/demo/ode/reaction/Makefile

src/demo/ode/reaction/README

src/demo/ode/reaction/bench-parameters.xml

src/demo/ode/reaction/bench.log

src/demo/ode/reaction/checkerror-debug.py

src/demo/ode/reaction/checkerror.py

src/demo/ode/reaction/debug-parameters.xml

src/demo/ode/reaction/domain.log

src/demo/ode/reaction/main.cpp

src/demo/ode/reaction/plot_reaction.py

src/demo/ode/reaction/reference-gamma-100.data

src/demo/ode/reaction/reference-gamma-1000-N-20.data

src/demo/ode/reaction/reference-gamma-1000.data

src/demo/ode/reaction/run-bench.py

src/demo/ode/reaction/run-debug.py

src/demo/ode/reaction/run-domain.py

src/demo/ode/residual

src/demo/ode/residual/Makefile

src/demo/ode/residual/main.cpp

src/demo/ode/residual/plot.py

src/demo/ode/residual/run-debug.py

src/demo/ode/stiff/plotsolution.m

src/demo/ode/test

src/demo/ode/test/Makefile

src/demo/ode/test/main.cpp

src/demo/ode/wave

src/demo/ode/wave/Makefile

src/demo/ode/wave/main.cpp

src/demo/ode/wave/plotsolution.m

src/demo/ode/wave/slit.xml.gz

src/demo/scripting/pydolfin/solvers/elasticity

src/demo/scripting/pydolfin/solvers/elasticity/Elasticity.form

src/demo/scripting/pydolfin/solvers/elasticity/ElasticityMass.form

src/demo/scripting/pydolfin/solvers/elasticity/elasticitysolver.py

src/demo/scripting/pydolfin/solvers/elasticity/tetmesh-4.xml.gz

src/demo/scripting/pydolfin/solvers/poisson/Poisson2D.form

src/demo/solvers/convdiff/visualize_convdiff.py

src/demo/solvers/elasticity-updated/plotsolution.sh

src/demo/solvers/elasticity-updated/template.pov

src/demo/solvers/elasticity/pde2ko.m

src/demo/solvers/elasticity/plotsolution.sh

src/demo/solvers/elasticity/template.pov

src/demo/solvers/elasticity/writeko.m

src/demo/solvers/elasticity/writeko.sh

src/demo/solvers/stokes/plotsolution.m

src/kernel/common/dolfin/PArray.h

src/kernel/common/dolfin/PList.h

src/kernel/common/dolfin/Table.h

src/kernel/common/dolfin/Tensor.h

src/kernel/fem/LoadVector.cpp

src/kernel/fem/MassMatrix.cpp

src/kernel/fem/StiffnessMatrix.cpp

src/kernel/fem/dolfin/LoadVector.h

src/kernel/fem/dolfin/MassMatrix.h

src/kernel/fem/dolfin/MassMatrix2D.form

src/kernel/fem/dolfin/MassMatrix2D.h

src/kernel/fem/dolfin/MassMatrix3D.form

src/kernel/fem/dolfin/MassMatrix3D.h

src/kernel/fem/dolfin/StiffnessMatrix.h

src/kernel/fem/dolfin/StiffnessMatrix2D.form

src/kernel/fem/dolfin/StiffnessMatrix2D.h

src/kernel/fem/dolfin/StiffnessMatrix3D.form

src/kernel/fem/dolfin/StiffnessMatrix3D.h

src/kernel/la/EigenvalueSolver.cpp

src/kernel/la/KrylovSolver.cpp

src/kernel/la/LinearSolver.cpp

src/kernel/la/Matrix.cpp

src/kernel/la/Preconditioner.cpp

src/kernel/la/Vector.cpp

src/kernel/la/VirtualMatrix.cpp

src/kernel/la/dolfin/EigenvalueSolver.h

src/kernel/la/dolfin/VirtualMatrix.h

src/kernel/mesh/Boundary.cpp

src/kernel/mesh/BoundaryData.cpp

src/kernel/mesh/BoundaryInit.cpp

src/kernel/mesh/CellIterator.cpp

src/kernel/mesh/EdgeIterator.cpp

src/kernel/mesh/EdgeRefData.cpp

src/kernel/mesh/FaceIterator.cpp

src/kernel/mesh/GenericCell.cpp

src/kernel/mesh/MeshInit.cpp

src/kernel/mesh/MeshIterator.cpp

src/kernel/mesh/MeshRefinement.cpp

src/kernel/mesh/TetMeshRefinement.cpp

src/kernel/mesh/TriMeshRefinement.cpp

src/kernel/mesh/VertexIterator.cpp

src/kernel/mesh/dolfin/Boundary.h

src/kernel/mesh/dolfin/BoundaryData.h

src/kernel/mesh/dolfin/BoundaryInit.h

src/kernel/mesh/dolfin/CellIterator.h

src/kernel/mesh/dolfin/CellRefData.h

src/kernel/mesh/dolfin/EdgeIterator.h

src/kernel/mesh/dolfin/EdgeRefData.h

src/kernel/mesh/dolfin/FaceIterator.h

src/kernel/mesh/dolfin/GenericCell.h

src/kernel/mesh/dolfin/MeshInit.h

src/kernel/mesh/dolfin/MeshIterator.h

src/kernel/mesh/dolfin/MeshRefinement.h

src/kernel/mesh/dolfin/TetMeshRefinement.h

src/kernel/mesh/dolfin/TriMeshRefinement.h

src/kernel/mesh/dolfin/VertexIterator.h

src/kernel/ode/ParticleSystem.cpp

src/kernel/ode/dolfin/ParticleSystem.h

src/post

src/post/Makefile.am

src/post/Makefile.in

src/post/links.sh

src/post/post.sh

src/pre

src/pre/Makefile.am

src/pre/Makefile.in

src/pydolfin/SettingsGlue.cpp

src/pydolfin/SettingsGlue.h

src/pydolfin/dolfin.py

src/pydolfin/dolfin_wrap.cpp

src/pydolfin/dolfin_wrap.h

src/test/L2Norm.form

src/test/L2Norm.h

src/test/Makefile

src/test/dolfin-1.xml.gz

src/test/function.xml

src/test/main.cpp

src/test/reftri.xml.gz

files modified:
.hgignore

AUTHORS

ChangeLog

Makefile.am

Makefile.in

README

TODO

aclocal.m4

configure

configure.ac

data/meshes/cylinder.xml.gz

data/meshes/dolfin-1.xml.gz

data/meshes/dolfin-2.xml.gz

data/meshes/gear.xml.gz

data/meshes/slit.xml.gz

data/meshes/tetmesh-1.xml.gz

data/meshes/tetmesh-16.xml.gz

data/meshes/tetmesh-4.xml.gz

data/meshes/tetmesh-8.xml.gz

data/meshes/tetmesh_24_8_8.xml.gz

data/meshes/tetmesh_LES_10_2_7.xml.gz

data/meshes/tetmesh_backward_facing_step_32_8_8.xml.gz

data/meshes/tetmesh_bluff_body_32_8_8.xml.gz

data/meshes/tetrahedron.xml.gz

data/meshes/triangle.xml.gz

data/meshes/trimesh-1.xml.gz

data/meshes/trimesh-128.xml.gz

data/meshes/trimesh-16.xml.gz

data/meshes/trimesh-2.xml.gz

data/meshes/trimesh-32.xml.gz

data/meshes/trimesh-4.xml.gz

data/meshes/trimesh-64.xml.gz

data/meshes/trimesh-8.xml.gz

data/meshes/two_particles.xml.gz

doc/manual/bibliography.bib

doc/manual/chapters/about.tex

doc/manual/chapters/contributing-dolfin.tex

doc/manual/chapters/contributing.tex

doc/manual/chapters/functions.tex

doc/manual/chapters/installation-downloading.tex

doc/manual/chapters/installation.tex

doc/manual/chapters/io.tex

doc/manual/chapters/linearalgebra.tex

doc/manual/chapters/logsystem.tex

doc/manual/chapters/mesh.tex

doc/manual/chapters/nonlinearsolver.tex

doc/manual/chapters/ode.tex

doc/manual/chapters/parameters.tex

doc/manual/chapters/pde.tex

doc/manual/chapters/quickstart.tex

doc/manual/chapters/referenceelements.tex

doc/manual/chapters/solvers.tex

doc/manual/dolfin-user-manual.pdf

doc/manual/dolfin-user-manual.tex

doc/manual/eps/poisson.eps

doc/manual/fenicsmanual.cls

doc/progr_reference/html/Doxyfile

doc/progr_reference/html/myfooter.inc

doc/progr_reference/html/myheader.inc

m4/ac_python_devel.m4

scripts/compileforms

scripts/makedist

src/Makefile.am

src/Makefile.in

src/demo/Makefile.am

src/demo/Makefile.in

src/demo/fem/Makefile.am

src/demo/fem/Makefile.in

src/demo/fem/simple/Makefile

src/demo/fem/simple/main.cpp

src/demo/fem/simple/tetrahedron.xml.gz

src/demo/mesh/Makefile

src/demo/mesh/main.cpp

src/demo/mesh/mesh2D.xml.gz

src/demo/mesh/mesh3D.xml.gz

src/demo/nls/Makefile.am

src/demo/nls/Makefile.in

src/demo/nls/nonlinearpoisson/Makefile

src/demo/nls/nonlinearpoisson/NonlinearPoisson.form

src/demo/nls/nonlinearpoisson/NonlinearPoisson.h

src/demo/nls/nonlinearpoisson/main.cpp

src/demo/ode/Makefile.am

src/demo/ode/Makefile.in

src/demo/ode/complex/Makefile

src/demo/ode/complex/main.cpp

src/demo/ode/homotopy/Makefile.in

src/demo/ode/homotopy/economy/Makefile

src/demo/ode/homotopy/economy/ces.cpp

src/demo/ode/homotopy/economy/test.cpp

src/demo/ode/homotopy/simple/Makefile

src/demo/ode/homotopy/simple/main.cpp

src/demo/ode/lorenz/Makefile

src/demo/ode/lorenz/main.cpp

src/demo/ode/method-weights/Makefile

src/demo/ode/stiff/Makefile

src/demo/ode/stiff/TestProblem1.h

src/demo/ode/stiff/TestProblem2.h

src/demo/ode/stiff/TestProblem3.h

src/demo/ode/stiff/TestProblem4.h

src/demo/ode/stiff/TestProblem5.h

src/demo/ode/stiff/TestProblem6.h

src/demo/ode/stiff/TestProblem7.h

src/demo/ode/stiff/TestProblem8.h

src/demo/ode/stiff/TestProblem9.h

src/demo/ode/stiff/main.cpp

src/demo/ode/stiff/run-tests.sh

src/demo/pde/Makefile.am

src/demo/pde/Makefile.in

src/demo/pde/convection-diffusion/ConvectionDiffusion.form

src/demo/pde/convection-diffusion/ConvectionDiffusion.h

src/demo/pde/convection-diffusion/Makefile

src/demo/pde/convection-diffusion/Stokes.form

src/demo/pde/convection-diffusion/Stokes.h

src/demo/pde/convection-diffusion/main.cpp

src/demo/pde/elasticity/Elasticity.form

src/demo/pde/elasticity/Elasticity.h

src/demo/pde/elasticity/ElasticityStrain.form

src/demo/pde/elasticity/ElasticityStrain.h

src/demo/pde/elasticity/Makefile

src/demo/pde/elasticity/main.cpp

src/demo/pde/nonlinear/Makefile.in

src/demo/pde/nonlinear/nonlinear-poisson/Makefile

src/demo/pde/nonlinear/nonlinear-poisson/NonlinearPoisson.form

src/demo/pde/nonlinear/nonlinear-poisson/NonlinearPoisson.h

src/demo/pde/nonlinear/nonlinear-poisson/main.cpp

src/demo/pde/poisson/Makefile

src/demo/pde/poisson/Poisson.form

src/demo/pde/poisson/Poisson.h

src/demo/pde/poisson/README

src/demo/pde/poisson/main.cpp

src/demo/pde/stokes/Makefile.in

src/demo/pde/stokes/stabilized/Makefile

src/demo/pde/stokes/stabilized/Stokes.form

src/demo/pde/stokes/stabilized/Stokes.h

src/demo/pde/stokes/stabilized/main.cpp

src/demo/pde/stokes/taylor-hood/Makefile

src/demo/pde/stokes/taylor-hood/Stokes.form

src/demo/pde/stokes/taylor-hood/Stokes.h

src/demo/pde/stokes/taylor-hood/main.cpp

src/demo/quadrature/Makefile

src/demo/scripting/Makefile.in

src/demo/scripting/pydolfin/cow05b.xml.gz

src/demo/scripting/pydolfin/cow05c.xml.gz

src/demo/scripting/pydolfin/lademo.py

src/demo/scripting/pydolfin/meshdemo.py

src/demo/scripting/pydolfin/odedemo-harmonic.py

src/demo/scripting/pydolfin/odedemo.py

src/demo/scripting/pydolfin/settingsdemo.py

src/demo/scripting/pydolfin/solvers/poisson/poissonsolver.py

src/demo/scripting/pydolfin/solvers/stokes/L2Error.form

src/demo/scripting/pydolfin/solvers/stokes/Stokes2D.form

src/demo/scripting/pydolfin/solvers/stokes/stokessolver.py

src/demo/solvers/Makefile.am

src/demo/solvers/Makefile.in

src/demo/solvers/convdiff/Makefile

src/demo/solvers/convdiff/dolfin.xml.gz

src/demo/solvers/convdiff/main.cpp

src/demo/solvers/elasticity-updated/Makefile

src/demo/solvers/elasticity-updated/cow05.xml.gz

src/demo/solvers/elasticity-updated/cow06.xml.gz

src/demo/solvers/elasticity-updated/cow07.xml.gz

src/demo/solvers/elasticity-updated/main.cpp

src/demo/solvers/elasticity-updated/roterror01.xml.gz

src/demo/solvers/elasticity/Makefile

src/demo/solvers/elasticity/main.cpp

src/demo/solvers/elasticity/tetmesh-4.xml.gz

src/demo/solvers/elasticity/tetmesh-8.xml.gz

src/demo/solvers/heat/Makefile

src/demo/solvers/heat/main.cpp

src/demo/solvers/navierstokes/Makefile.in

src/demo/solvers/navierstokes/benchmark/Makefile

src/demo/solvers/navierstokes/benchmark/README

src/demo/solvers/navierstokes/benchmark/cylinder_2d_bmk.xml.gz

src/demo/solvers/navierstokes/benchmark/cylinder_3d_bmk.xml.gz

src/demo/solvers/navierstokes/benchmark/main.cpp

src/demo/solvers/poisson/Makefile

src/demo/solvers/poisson/main.cpp

src/demo/solvers/stokes/Makefile

src/demo/solvers/stokes/main.cpp

src/greeting/Makefile.in

src/kernel/Makefile.am

src/kernel/Makefile.in

src/kernel/common/Makefile.am

src/kernel/common/Makefile.in

src/kernel/common/Variable.cpp

src/kernel/common/dolfin/Array.h

src/kernel/common/dolfin/Makefile.am

src/kernel/common/dolfin/Makefile.in

src/kernel/common/dolfin/Variable.h

src/kernel/common/dolfin/dolfin_common.h

src/kernel/fem/AffineMap.cpp

src/kernel/fem/BoundaryValue.cpp

src/kernel/fem/FEM.cpp

src/kernel/fem/FiniteElement.cpp

src/kernel/fem/FiniteElementSpec.cpp

src/kernel/fem/Makefile.am

src/kernel/fem/Makefile.in

src/kernel/fem/dolfin/AffineMap.h

src/kernel/fem/dolfin/BoundaryValue.h

src/kernel/fem/dolfin/FEM.h

src/kernel/fem/dolfin/FiniteElement.h

src/kernel/fem/dolfin/Makefile.am

src/kernel/fem/dolfin/Makefile.in

src/kernel/fem/dolfin/P1tet.h

src/kernel/fem/dolfin/P1tri.h

src/kernel/fem/dolfin/P2tet.h

src/kernel/fem/dolfin/P2tri.h

src/kernel/fem/dolfin/P3tet.h

src/kernel/fem/dolfin/P3tri.h

src/kernel/fem/dolfin/P4tet.h

src/kernel/fem/dolfin/P4tri.h

src/kernel/fem/dolfin/P5tet.h

src/kernel/fem/dolfin/P5tri.h

src/kernel/fem/dolfin/dolfin_fem.h

src/kernel/form/BilinearForm.cpp

src/kernel/form/Form.cpp

src/kernel/form/Functional.cpp

src/kernel/form/LinearForm.cpp

src/kernel/form/Makefile.in

src/kernel/form/dolfin/BilinearForm.h

src/kernel/form/dolfin/Form.h

src/kernel/form/dolfin/Functional.h

src/kernel/form/dolfin/LinearForm.h

src/kernel/form/dolfin/Makefile.in

src/kernel/form/dolfin/dolfin_form.h

src/kernel/function/ConstantFunction.cpp

src/kernel/function/DiscreteFunction.cpp

src/kernel/function/Function.cpp

src/kernel/function/FunctionPointerFunction.cpp

src/kernel/function/GenericFunction.cpp

src/kernel/function/LocalFunctionData.cpp

src/kernel/function/Makefile.in

src/kernel/function/UserFunction.cpp

src/kernel/function/dolfin/ConstantFunction.h

src/kernel/function/dolfin/DiscreteFunction.h

src/kernel/function/dolfin/Function.h

src/kernel/function/dolfin/FunctionPointerFunction.h

src/kernel/function/dolfin/GenericFunction.h

src/kernel/function/dolfin/LocalFunctionData.h

src/kernel/function/dolfin/Makefile.in

src/kernel/function/dolfin/SpecialFunctions.h

src/kernel/function/dolfin/UserFunction.h

src/kernel/io/File.cpp

src/kernel/io/GenericFile.cpp

src/kernel/io/GiDFile.cpp

src/kernel/io/MFile.cpp

src/kernel/io/MTXFile.cpp

src/kernel/io/Makefile.am

src/kernel/io/Makefile.in

src/kernel/io/MatlabFile.cpp

src/kernel/io/OctaveFile.cpp

src/kernel/io/OpenDXFile.cpp

src/kernel/io/PythonFile.cpp

src/kernel/io/TecplotFile.cpp

src/kernel/io/VTKFile.cpp

src/kernel/io/XMLBLASFormData.cpp

src/kernel/io/XMLFile.cpp

src/kernel/io/XMLFiniteElementSpec.cpp

src/kernel/io/XMLFunction.cpp

src/kernel/io/XMLMatrix.cpp

src/kernel/io/XMLMesh.cpp

src/kernel/io/XMLObject.cpp

src/kernel/io/XMLParameterList.cpp

src/kernel/io/XMLVector.cpp

src/kernel/io/dolfin/File.h

src/kernel/io/dolfin/GenericFile.h

src/kernel/io/dolfin/MFile.h

src/kernel/io/dolfin/MTXFile.h

src/kernel/io/dolfin/Makefile.am

src/kernel/io/dolfin/Makefile.in

src/kernel/io/dolfin/MatlabFile.h

src/kernel/io/dolfin/OctaveFile.h

src/kernel/io/dolfin/VTKFile.h

src/kernel/io/dolfin/XMLBLASFormData.h

src/kernel/io/dolfin/XMLFile.h

src/kernel/io/dolfin/XMLFiniteElementSpec.h

src/kernel/io/dolfin/XMLFunction.h

src/kernel/io/dolfin/XMLMatrix.h

src/kernel/io/dolfin/XMLMesh.h

src/kernel/io/dolfin/XMLObject.h

src/kernel/io/dolfin/XMLParameterList.h

src/kernel/io/dolfin/XMLVector.h

src/kernel/la/GMRES.cpp

src/kernel/la/LU.cpp

src/kernel/la/Makefile.am

src/kernel/la/Makefile.in

src/kernel/la/PETScManager.cpp

src/kernel/la/dolfin/GMRES.h

src/kernel/la/dolfin/KrylovSolver.h

src/kernel/la/dolfin/LU.h

src/kernel/la/dolfin/LinearSolver.h

src/kernel/la/dolfin/Makefile.am

src/kernel/la/dolfin/Makefile.in

src/kernel/la/dolfin/Matrix.h

src/kernel/la/dolfin/PETScManager.h

src/kernel/la/dolfin/Preconditioner.h

src/kernel/la/dolfin/Vector.h

src/kernel/la/dolfin/dolfin_la.h

src/kernel/log/CursesLogger.cpp

src/kernel/log/LogStream.cpp

src/kernel/log/Logger.cpp

src/kernel/log/Makefile.in

src/kernel/log/Progress.cpp

src/kernel/log/dolfin/CursesLogger.h

src/kernel/log/dolfin/Makefile.in

src/kernel/main/Makefile.in

src/kernel/main/dolfin.h

src/kernel/main/dolfin/Makefile.in

src/kernel/main/dolfin_doc.h

src/kernel/main/init.cpp

src/kernel/math/Lagrange.cpp

src/kernel/math/Makefile.in

src/kernel/math/basic.cpp

src/kernel/math/dolfin/Makefile.in

src/kernel/math/dolfin/basic.h

src/kernel/mesh/Cell.cpp

src/kernel/mesh/Edge.cpp

src/kernel/mesh/Face.cpp

src/kernel/mesh/Makefile.am

src/kernel/mesh/Makefile.in

src/kernel/mesh/Mesh.cpp

src/kernel/mesh/MeshData.cpp

src/kernel/mesh/MeshHierarchy.cpp

src/kernel/mesh/Point.cpp

src/kernel/mesh/Tetrahedron.cpp

src/kernel/mesh/Triangle.cpp

src/kernel/mesh/UnitCube.cpp

src/kernel/mesh/UnitSquare.cpp

src/kernel/mesh/Vertex.cpp

src/kernel/mesh/dolfin/Cell.h

src/kernel/mesh/dolfin/Edge.h

src/kernel/mesh/dolfin/Face.h

src/kernel/mesh/dolfin/Makefile.am

src/kernel/mesh/dolfin/Makefile.in

src/kernel/mesh/dolfin/Mesh.h

src/kernel/mesh/dolfin/MeshData.h

src/kernel/mesh/dolfin/MeshHierarchy.h

src/kernel/mesh/dolfin/Point.h

src/kernel/mesh/dolfin/Tetrahedron.h

src/kernel/mesh/dolfin/Triangle.h

src/kernel/mesh/dolfin/UnitCube.h

src/kernel/mesh/dolfin/UnitSquare.h

src/kernel/mesh/dolfin/Vertex.h

src/kernel/mesh/dolfin/dolfin_mesh.h

src/kernel/nls/Makefile.in

src/kernel/nls/NewtonSolver.cpp

src/kernel/nls/NonlinearProblem.cpp

src/kernel/nls/dolfin/Makefile.in

src/kernel/nls/dolfin/NewtonSolver.h

src/kernel/nls/dolfin/NonlinearProblem.h

src/kernel/ode/Adaptivity.cpp

src/kernel/ode/ComplexODE.cpp

src/kernel/ode/Dependencies.cpp

src/kernel/ode/Homotopy.cpp

src/kernel/ode/HomotopyJacobian.cpp

src/kernel/ode/HomotopyODE.cpp

src/kernel/ode/Makefile.am

src/kernel/ode/Makefile.in

src/kernel/ode/Method.cpp

src/kernel/ode/MonoAdaptiveFixedPointSolver.cpp

src/kernel/ode/MonoAdaptiveJacobian.cpp

src/kernel/ode/MonoAdaptiveNewtonSolver.cpp

src/kernel/ode/MonoAdaptiveTimeSlab.cpp

src/kernel/ode/MonoAdaptivity.cpp

src/kernel/ode/MultiAdaptiveFixedPointSolver.cpp

src/kernel/ode/MultiAdaptiveJacobian.cpp

src/kernel/ode/MultiAdaptiveNewtonSolver.cpp

src/kernel/ode/MultiAdaptivePreconditioner.cpp

src/kernel/ode/MultiAdaptiveTimeSlab.cpp

src/kernel/ode/MultiAdaptivity.cpp

src/kernel/ode/ODE.cpp

src/kernel/ode/ODESolver.cpp

src/kernel/ode/Partition.cpp

src/kernel/ode/ReducedModel.cpp

src/kernel/ode/TimeSlab.cpp

src/kernel/ode/TimeSlabJacobian.cpp

src/kernel/ode/TimeSlabSolver.cpp

src/kernel/ode/TimeStepper.cpp

src/kernel/ode/UpdatedMultiAdaptiveJacobian.cpp

src/kernel/ode/cGqMethod.cpp

src/kernel/ode/dGqMethod.cpp

src/kernel/ode/dolfin/ComplexODE.h

src/kernel/ode/dolfin/Dependencies.h

src/kernel/ode/dolfin/Homotopy.h

src/kernel/ode/dolfin/HomotopyJacobian.h

src/kernel/ode/dolfin/HomotopyODE.h

src/kernel/ode/dolfin/Makefile.am

src/kernel/ode/dolfin/Makefile.in

src/kernel/ode/dolfin/Method.h

src/kernel/ode/dolfin/MonoAdaptiveFixedPointSolver.h

src/kernel/ode/dolfin/MonoAdaptiveJacobian.h

src/kernel/ode/dolfin/MonoAdaptiveNewtonSolver.h

src/kernel/ode/dolfin/MonoAdaptiveTimeSlab.h

src/kernel/ode/dolfin/MultiAdaptiveFixedPointSolver.h

src/kernel/ode/dolfin/MultiAdaptiveJacobian.h

src/kernel/ode/dolfin/MultiAdaptiveNewtonSolver.h

src/kernel/ode/dolfin/MultiAdaptivePreconditioner.h

src/kernel/ode/dolfin/MultiAdaptiveTimeSlab.h

src/kernel/ode/dolfin/MultiAdaptivity.h

src/kernel/ode/dolfin/ODE.h

src/kernel/ode/dolfin/ReducedModel.h

src/kernel/ode/dolfin/TimeSlab.h

src/kernel/ode/dolfin/TimeSlabJacobian.h

src/kernel/ode/dolfin/TimeSlabSolver.h

src/kernel/ode/dolfin/UpdatedMultiAdaptiveJacobian.h

src/kernel/ode/dolfin/cGqMethod.h

src/kernel/ode/dolfin/dGqMethod.h

src/kernel/ode/dolfin/dolfin_ode.h

src/kernel/parameter/Makefile.in

src/kernel/parameter/ParameterSystem.cpp

src/kernel/parameter/Parametrized.cpp

src/kernel/parameter/dolfin/DefaultParameters.h

src/kernel/parameter/dolfin/Makefile.in

src/kernel/parameter/dolfin/Parameter.h

src/kernel/parameter/dolfin/Parametrized.h

src/kernel/pde/LinearPDE.cpp

src/kernel/pde/Makefile.am

src/kernel/pde/Makefile.in

src/kernel/pde/NonlinearPDE.cpp

src/kernel/pde/PDE.cpp

src/kernel/pde/dolfin/GenericPDE.h

src/kernel/pde/dolfin/LinearPDE.h

src/kernel/pde/dolfin/Makefile.in

src/kernel/pde/dolfin/NonlinearPDE.h

src/kernel/pde/dolfin/PDE.h

src/kernel/pde/dolfin/dolfin_pde.h

src/kernel/quadrature/GaussQuadrature.cpp

src/kernel/quadrature/GaussianQuadrature.cpp

src/kernel/quadrature/LobattoQuadrature.cpp

src/kernel/quadrature/Makefile.in

src/kernel/quadrature/Quadrature.cpp

src/kernel/quadrature/RadauQuadrature.cpp

src/kernel/quadrature/dolfin/GaussianQuadrature.h

src/kernel/quadrature/dolfin/Makefile.in

src/kernel/quadrature/dolfin/Quadrature.h

src/lib/Makefile.in

src/modules/Makefile.am

src/modules/Makefile.in

src/modules/convdiff/ConvectionDiffusionSolver.cpp

src/modules/convdiff/Makefile.in

src/modules/convdiff/dolfin/ConvectionDiffusion2D.form

src/modules/convdiff/dolfin/ConvectionDiffusion2D.h

src/modules/convdiff/dolfin/ConvectionDiffusion3D.form

src/modules/convdiff/dolfin/ConvectionDiffusion3D.h

src/modules/convdiff/dolfin/ConvectionDiffusionSolver.h

src/modules/convdiff/dolfin/Makefile.in

src/modules/dolfin/Makefile.in

src/modules/dolfin/Solver.h

src/modules/dolfin/dolfin_modules.h

src/modules/elasticity-updated/ElasticityUpdatedSolver.cpp

src/modules/elasticity-updated/Makefile.in

src/modules/elasticity-updated/dolfin/ElasticityUpdated.form

src/modules/elasticity-updated/dolfin/ElasticityUpdated.h

src/modules/elasticity-updated/dolfin/ElasticityUpdatedMass.form

src/modules/elasticity-updated/dolfin/ElasticityUpdatedMass.h

src/modules/elasticity-updated/dolfin/ElasticityUpdatedSigma.form

src/modules/elasticity-updated/dolfin/ElasticityUpdatedSigma.h

src/modules/elasticity-updated/dolfin/ElasticityUpdatedSolver.h

src/modules/elasticity-updated/dolfin/Makefile.in

src/modules/elasticity/ElasticitySolver.cpp

src/modules/elasticity/Makefile.in

src/modules/elasticity/dolfin/Elasticity-alternative.form

src/modules/elasticity/dolfin/Elasticity.form

src/modules/elasticity/dolfin/Elasticity.h

src/modules/elasticity/dolfin/ElasticityMass.form

src/modules/elasticity/dolfin/ElasticityMass.h

src/modules/elasticity/dolfin/ElasticitySolver.h

src/modules/elasticity/dolfin/Makefile.in

src/modules/heat/HeatSolver.cpp

src/modules/heat/Makefile.in

src/modules/heat/dolfin/Heat.form

src/modules/heat/dolfin/Heat.h

src/modules/heat/dolfin/HeatSolver.h

src/modules/heat/dolfin/Makefile.in

src/modules/navierstokes/Makefile.in

src/modules/navierstokes/NSESolver.cpp

src/modules/navierstokes/dolfin/Makefile.in

src/modules/navierstokes/dolfin/NSEContinuity2D.form

src/modules/navierstokes/dolfin/NSEContinuity2D.h

src/modules/navierstokes/dolfin/NSEContinuity3D.form

src/modules/navierstokes/dolfin/NSEContinuity3D.h

src/modules/navierstokes/dolfin/NSEMomentum2D.form

src/modules/navierstokes/dolfin/NSEMomentum2D.h

src/modules/navierstokes/dolfin/NSEMomentum3D.form

src/modules/navierstokes/dolfin/NSEMomentum3D.h

src/modules/navierstokes/dolfin/NSESolver.h

src/modules/poisson/Makefile.in

src/modules/poisson/PoissonSolver.cpp

src/modules/poisson/dolfin/Makefile.in

src/modules/poisson/dolfin/Poisson2D.form

src/modules/poisson/dolfin/Poisson2D.h

src/modules/poisson/dolfin/Poisson3D.form

src/modules/poisson/dolfin/Poisson3D.h

src/modules/poisson/dolfin/PoissonSolver.h

src/modules/stokes/Makefile.in

src/modules/stokes/StokesSolver.cpp

src/modules/stokes/dolfin/L2Error.form

src/modules/stokes/dolfin/L2Error.h

src/modules/stokes/dolfin/Makefile.in

src/modules/stokes/dolfin/Stokes2D.form

src/modules/stokes/dolfin/Stokes2D.h

src/modules/stokes/dolfin/Stokes3D.form

src/modules/stokes/dolfin/Stokes3D.h

src/modules/stokes/dolfin/StokesSolver.h

src/pydolfin/Makefile.am

src/pydolfin/Makefile.in

src/pydolfin/README

src/pydolfin/__init__.py

src/pydolfin/dolfin.i

src/pydolfin/import_form.py

src/utils/Makefile.am

src/utils/Makefile.in

src/utils/convert/Makefile.in

src/utils/convert/dolfin-convert

src/utils/inp2dx/Makefile.in

Show diffs side-by-side

added added

removed removed

src/demo/fem/convergence/Poisson2D_3.h

1

// Automatically generated by FFC, the FEniCS Form Compiler, version 0.2.5.

2

// For further information, go to http://www/fenics.org/ffc/.

3

// Licensed under the GNU GPL Version 2.

4

5

#ifndef __POISSON2D_3_H

6

#define __POISSON2D_3_H

7

8

#include <dolfin/Mesh.h>

9

#include <dolfin/Cell.h>

10

#include <dolfin/Point.h>

11

#include <dolfin/Vector.h>

12

#include <dolfin/AffineMap.h>

13

#include <dolfin/FiniteElement.h>

14

#include <dolfin/FiniteElementSpec.h>

15

#include <dolfin/LinearForm.h>

16

#include <dolfin/BilinearForm.h>

17

18

namespace dolfin { namespace Poisson2D_3 {

19

20

/// This class contains the form to be evaluated, including

21

/// contributions from the interior and boundary of the domain.

22

23

class BilinearForm : public dolfin::BilinearForm

24

{

25

public:

26

27

class TestElement : public dolfin::FiniteElement

28

{

29

public:

30

31

TestElement() : dolfin::FiniteElement(), tensordims(0), subelements(0)

32

{

33

// Element is scalar, don't need to initialize tensordims

34

35

// Element is simple, don't need to initialize subelements

36

}

37

38

~TestElement()

39

{

40

if ( tensordims ) delete [] tensordims;

41

if ( subelements )

42

{

43

for (unsigned int i = 0; i < elementdim(); i++)

44

delete subelements[i];

45

delete [] subelements;

46

}

47

}

48

49

inline unsigned int spacedim() const

50

{

51

return 10;

52

}

53

54

inline unsigned int shapedim() const

55

{

56

return 2;

57

}

58

59

inline unsigned int tensordim(unsigned int i) const

60

{

61

dolfin_error("Element is scalar.");

62

return 0;

63

}

64

65

inline unsigned int elementdim() const

66

{

67

return 1;

68

}

69

70

inline unsigned int rank() const

71

{

72

return 0;

73

}

74

75

void nodemap(int nodes[], const Cell& cell, const Mesh& mesh) const

76

{

77

static unsigned int edge_reordering_0[2][2] = {{0, 1}, {1, 0}};

78

nodes[0] = cell.vertexID(0);

79

nodes[1] = cell.vertexID(1);

80

nodes[2] = cell.vertexID(2);

81

int alignment = cell.edgeAlignment(0);

82

int offset = mesh.numVertices();

83

nodes[3] = offset + 2*cell.edgeID(0) + edge_reordering_0[alignment][0];

84

nodes[4] = offset + 2*cell.edgeID(0) + edge_reordering_0[alignment][1];

85

alignment = cell.edgeAlignment(1);

86

nodes[5] = offset + 2*cell.edgeID(1) + edge_reordering_0[alignment][0];

87

nodes[6] = offset + 2*cell.edgeID(1) + edge_reordering_0[alignment][1];

88

alignment = cell.edgeAlignment(2);

89

nodes[7] = offset + 2*cell.edgeID(2) + edge_reordering_0[alignment][0];

90

nodes[8] = offset + 2*cell.edgeID(2) + edge_reordering_0[alignment][1];

91

offset = offset + 2*mesh.numEdges();

92

nodes[9] = offset + cell.id();

93

}

94

95

void pointmap(Point points[], unsigned int components[], const AffineMap& map) const

96

{

97

points[0] = map(0.000000000000000e+00, 0.000000000000000e+00);

98

points[1] = map(1.000000000000000e+00, 0.000000000000000e+00);

99

points[2] = map(0.000000000000000e+00, 1.000000000000000e+00);

100

points[3] = map(6.666666666666667e-01, 3.333333333333333e-01);

101

points[4] = map(3.333333333333334e-01, 6.666666666666666e-01);

102

points[5] = map(0.000000000000000e+00, 6.666666666666667e-01);

103

points[6] = map(0.000000000000000e+00, 3.333333333333334e-01);

104

points[7] = map(3.333333333333333e-01, 0.000000000000000e+00);

105

points[8] = map(6.666666666666666e-01, 0.000000000000000e+00);

106

points[9] = map(3.333333333333333e-01, 3.333333333333333e-01);

107

components[0] = 0;

108

components[1] = 0;

109

components[2] = 0;

110

components[3] = 0;

111

components[4] = 0;

112

components[5] = 0;

113

components[6] = 0;

114

components[7] = 0;

115

components[8] = 0;

116

components[9] = 0;

117

}

118

119

void vertexeval(real values[], unsigned int vertex, const real x[], const Mesh& mesh) const

120

{

121

// FIXME: Temporary fix for Lagrange elements

122

values[0] = x[vertex];

123

}

124

125

const FiniteElement& operator[] (unsigned int i) const

126

{

127

return *this;

128

}

129

130

FiniteElement& operator[] (unsigned int i)

131

{

132

return *this;

133

}

134

135

FiniteElementSpec spec() const

136

{

137

FiniteElementSpec s("Lagrange", "triangle", 3);

138

return s;

139

}

140

141

private:

142

143

unsigned int* tensordims;

144

FiniteElement** subelements;

145

146

};

147

148

class TrialElement : public dolfin::FiniteElement

149

{

150

public:

151

152

TrialElement() : dolfin::FiniteElement(), tensordims(0), subelements(0)

153

{

154

// Element is scalar, don't need to initialize tensordims

155

156

// Element is simple, don't need to initialize subelements

157

}

158

159

~TrialElement()

160

{

161

if ( tensordims ) delete [] tensordims;

162

if ( subelements )

163

{

164

for (unsigned int i = 0; i < elementdim(); i++)

165

delete subelements[i];

166

delete [] subelements;

167

}

168

}

169

170

inline unsigned int spacedim() const

171

{

172

return 10;

173

}

174

175

inline unsigned int shapedim() const

176

{

177

return 2;

178

}

179

180

inline unsigned int tensordim(unsigned int i) const

181

{

182

dolfin_error("Element is scalar.");

183

return 0;

184

}

185

186

inline unsigned int elementdim() const

187

{

188

return 1;

189

}

190

191

inline unsigned int rank() const

192

{

193

return 0;

194

}

195

196

void nodemap(int nodes[], const Cell& cell, const Mesh& mesh) const

197

{

198

static unsigned int edge_reordering_0[2][2] = {{0, 1}, {1, 0}};

199

nodes[0] = cell.vertexID(0);

200

nodes[1] = cell.vertexID(1);

201

nodes[2] = cell.vertexID(2);

202

int alignment = cell.edgeAlignment(0);

203

int offset = mesh.numVertices();

204

nodes[3] = offset + 2*cell.edgeID(0) + edge_reordering_0[alignment][0];

205

nodes[4] = offset + 2*cell.edgeID(0) + edge_reordering_0[alignment][1];

206

alignment = cell.edgeAlignment(1);

207

nodes[5] = offset + 2*cell.edgeID(1) + edge_reordering_0[alignment][0];

208

nodes[6] = offset + 2*cell.edgeID(1) + edge_reordering_0[alignment][1];

209

alignment = cell.edgeAlignment(2);

210

nodes[7] = offset + 2*cell.edgeID(2) + edge_reordering_0[alignment][0];

211

nodes[8] = offset + 2*cell.edgeID(2) + edge_reordering_0[alignment][1];

212

offset = offset + 2*mesh.numEdges();

213

nodes[9] = offset + cell.id();

214

}

215

216

void pointmap(Point points[], unsigned int components[], const AffineMap& map) const

217

{

218

points[0] = map(0.000000000000000e+00, 0.000000000000000e+00);

219

points[1] = map(1.000000000000000e+00, 0.000000000000000e+00);

220

points[2] = map(0.000000000000000e+00, 1.000000000000000e+00);

221

points[3] = map(6.666666666666667e-01, 3.333333333333333e-01);

222

points[4] = map(3.333333333333334e-01, 6.666666666666666e-01);

223

points[5] = map(0.000000000000000e+00, 6.666666666666667e-01);

224

points[6] = map(0.000000000000000e+00, 3.333333333333334e-01);

225

points[7] = map(3.333333333333333e-01, 0.000000000000000e+00);

226

points[8] = map(6.666666666666666e-01, 0.000000000000000e+00);

227

points[9] = map(3.333333333333333e-01, 3.333333333333333e-01);

228

components[0] = 0;

229

components[1] = 0;

230

components[2] = 0;

231

components[3] = 0;

232

components[4] = 0;

233

components[5] = 0;

234

components[6] = 0;

235

components[7] = 0;

236

components[8] = 0;

237

components[9] = 0;

238

}

239

240

void vertexeval(real values[], unsigned int vertex, const real x[], const Mesh& mesh) const

241

{

242

// FIXME: Temporary fix for Lagrange elements

243

values[0] = x[vertex];

244

}

245

246

const FiniteElement& operator[] (unsigned int i) const

247

{

248

return *this;

249

}

250

251

FiniteElement& operator[] (unsigned int i)

252

{

253

return *this;

254

}

255

256

FiniteElementSpec spec() const

257

{

258

FiniteElementSpec s("Lagrange", "triangle", 3);

259

return s;

260

}

261

262

private:

263

264

unsigned int* tensordims;

265

FiniteElement** subelements;

266

267

};

268

269

BilinearForm() : dolfin::BilinearForm(0)

270

{

271

// Create finite element for test space

272

_test = new TestElement();

273

274

// Create finite element for trial space

275

_trial = new TrialElement();

276

}

277

278

void eval(real block[], const AffineMap& map) const

279

{

280

// Compute geometry tensors

281

const real G0_0_0 = map.det*(map.g00*map.g00 + map.g01*map.g01);

282

const real G0_0_1 = map.det*(map.g00*map.g10 + map.g01*map.g11);

283

const real G0_1_0 = map.det*(map.g10*map.g00 + map.g11*map.g01);

284

const real G0_1_1 = map.det*(map.g10*map.g10 + map.g11*map.g11);

285

286

// Compute element tensor

287

block[0] = 4.249999999999998e-01*G0_0_0 + 4.249999999999995e-01*G0_0_1 + 4.249999999999996e-01*G0_1_0 + 4.249999999999993e-01*G0_1_1;

288

block[1] = -8.750000000000012e-02*G0_0_0 - 8.750000000000020e-02*G0_1_0;

289

block[2] = -8.750000000000001e-02*G0_0_1 - 8.750000000000005e-02*G0_1_1;

290

block[3] = -3.749999999999987e-02*G0_0_0 - 3.750000000000012e-02*G0_0_1 - 3.749999999999996e-02*G0_1_0 - 3.750000000000038e-02*G0_1_1;

291

block[4] = -3.750000000000003e-02*G0_0_0 - 3.749999999999992e-02*G0_0_1 - 3.750000000000013e-02*G0_1_0 - 3.749999999999980e-02*G0_1_1;

292

block[5] = 3.750000000000010e-02*G0_0_0 + 3.374999999999996e-01*G0_0_1 + 3.750000000000017e-02*G0_1_0 + 3.374999999999997e-01*G0_1_1;

293

block[6] = 3.750000000000024e-02*G0_0_0 - 6.749999999999989e-01*G0_0_1 + 3.749999999999991e-02*G0_1_0 - 6.749999999999988e-01*G0_1_1;

294

block[7] = -6.749999999999997e-01*G0_0_0 + 3.749999999999988e-02*G0_0_1 - 6.749999999999996e-01*G0_1_0 + 3.749999999999964e-02*G0_1_1;

295

block[8] = 3.375000000000001e-01*G0_0_0 + 3.750000000000022e-02*G0_0_1 + 3.375000000000001e-01*G0_1_0 + 3.750000000000051e-02*G0_1_1;

296

block[9] = 0.000000000000000e+00;

297

block[10] = -8.750000000000013e-02*G0_0_0 - 8.750000000000022e-02*G0_0_1;

298

block[11] = 4.249999999999997e-01*G0_0_0;

299

block[12] = 8.750000000000011e-02*G0_0_1;

300

block[13] = 3.750000000000030e-02*G0_0_0 + 7.124999999999994e-01*G0_0_1;

301

block[14] = 3.750000000000041e-02*G0_0_0 - 2.999999999999993e-01*G0_0_1;

302

block[15] = -3.750000000000044e-02*G0_0_0;

303

block[16] = -3.749999999999944e-02*G0_0_0;

304

block[17] = 3.375000000000000e-01*G0_0_0 + 3.000000000000002e-01*G0_0_1;

305

block[18] = -6.749999999999996e-01*G0_0_0 - 7.124999999999994e-01*G0_0_1;

306

block[19] = 0.000000000000000e+00;

307

block[20] = -8.750000000000001e-02*G0_1_0 - 8.750000000000005e-02*G0_1_1;

308

block[21] = 8.750000000000011e-02*G0_1_0;

309

block[22] = 4.249999999999997e-01*G0_1_1;

310

block[23] = -2.999999999999994e-01*G0_1_0 + 3.750000000000030e-02*G0_1_1;

311

block[24] = 7.124999999999992e-01*G0_1_0 + 3.749999999999998e-02*G0_1_1;

312

block[25] = -7.124999999999992e-01*G0_1_0 - 6.749999999999997e-01*G0_1_1;

313

block[26] = 2.999999999999999e-01*G0_1_0 + 3.375000000000001e-01*G0_1_1;

314

block[27] = -3.749999999999994e-02*G0_1_1;

315

block[28] = -3.750000000000028e-02*G0_1_1;

316

block[29] = 0.000000000000000e+00;

317

block[30] = -3.749999999999987e-02*G0_0_0 - 3.749999999999996e-02*G0_0_1 - 3.750000000000010e-02*G0_1_0 - 3.750000000000037e-02*G0_1_1;

318

block[31] = 3.750000000000024e-02*G0_0_0 + 7.124999999999992e-01*G0_1_0;

319

block[32] = -2.999999999999994e-01*G0_0_1 + 3.750000000000028e-02*G0_1_1;

320

block[33] = 1.687499999999997e+00*G0_0_0 + 8.437499999999991e-01*G0_0_1 + 8.437499999999991e-01*G0_1_0 + 1.687499999999998e+00*G0_1_1;

321

block[34] = -3.374999999999996e-01*G0_0_0 + 8.437499999999984e-01*G0_0_1 - 1.687499999999994e-01*G0_1_0 - 3.374999999999997e-01*G0_1_1;

322

block[35] = 3.374999999999995e-01*G0_0_0 + 1.687500000000000e-01*G0_0_1 + 1.687499999999993e-01*G0_1_0;

323

block[36] = 3.374999999999992e-01*G0_0_0 + 1.687499999999993e-01*G0_0_1 + 1.687500000000002e-01*G0_1_0;

324

block[37] = 1.687499999999996e-01*G0_0_1 + 1.687499999999996e-01*G0_1_0 + 3.374999999999994e-01*G0_1_1;

325

block[38] = -8.437499999999994e-01*G0_0_1 - 8.437499999999989e-01*G0_1_0 - 1.687499999999998e+00*G0_1_1;

326

block[39] = -2.024999999999996e+00*G0_0_0 - 1.012499999999998e+00*G0_0_1 - 1.012500000000000e+00*G0_1_0;

327

block[40] = -3.750000000000003e-02*G0_0_0 - 3.750000000000012e-02*G0_0_1 - 3.749999999999992e-02*G0_1_0 - 3.749999999999982e-02*G0_1_1;

328

block[41] = 3.750000000000041e-02*G0_0_0 - 2.999999999999993e-01*G0_1_0;

329

block[42] = 7.124999999999992e-01*G0_0_1 + 3.749999999999996e-02*G0_1_1;

330

block[43] = -3.374999999999995e-01*G0_0_0 - 1.687499999999994e-01*G0_0_1 + 8.437499999999984e-01*G0_1_0 - 3.374999999999997e-01*G0_1_1;

331

block[44] = 1.687499999999998e+00*G0_0_0 + 8.437499999999987e-01*G0_0_1 + 8.437499999999984e-01*G0_1_0 + 1.687499999999998e+00*G0_1_1;

332

block[45] = -1.687499999999998e+00*G0_0_0 - 8.437499999999996e-01*G0_0_1 - 8.437499999999987e-01*G0_1_0;

333

block[46] = 3.374999999999999e-01*G0_0_0 + 1.687500000000003e-01*G0_0_1 + 1.687499999999994e-01*G0_1_0;

334

block[47] = 1.687499999999998e-01*G0_0_1 + 1.687499999999996e-01*G0_1_0 + 3.374999999999996e-01*G0_1_1;

335

block[48] = 1.687499999999992e-01*G0_0_1 + 1.687499999999996e-01*G0_1_0 + 3.374999999999992e-01*G0_1_1;

336

block[49] = -1.012499999999998e+00*G0_0_1 - 1.012499999999998e+00*G0_1_0 - 2.024999999999997e+00*G0_1_1;

337

block[50] = 3.750000000000009e-02*G0_0_0 + 3.750000000000017e-02*G0_0_1 + 3.374999999999996e-01*G0_1_0 + 3.374999999999996e-01*G0_1_1;

338

block[51] = -3.750000000000045e-02*G0_0_0;

339

block[52] = -7.124999999999994e-01*G0_0_1 - 6.749999999999996e-01*G0_1_1;

340

block[53] = 3.374999999999995e-01*G0_0_0 + 1.687499999999993e-01*G0_0_1 + 1.687500000000000e-01*G0_1_0;

341

block[54] = -1.687499999999998e+00*G0_0_0 - 8.437499999999987e-01*G0_0_1 - 8.437499999999996e-01*G0_1_0;

342

block[55] = 1.687499999999998e+00*G0_0_0 + 8.437499999999996e-01*G0_0_1 + 8.437499999999996e-01*G0_1_0 + 1.687499999999999e+00*G0_1_1;

343

block[56] = -3.374999999999999e-01*G0_0_0 - 1.687500000000005e-01*G0_0_1 - 1.181249999999999e+00*G0_1_0 - 1.349999999999999e+00*G0_1_1;

344

block[57] = -1.687499999999997e-01*G0_0_1 - 1.687500000000000e-01*G0_1_0;

345

block[58] = -1.687499999999992e-01*G0_0_1 - 1.687499999999993e-01*G0_1_0;

346

block[59] = 1.012499999999998e+00*G0_0_1 + 1.012500000000000e+00*G0_1_0;

347

block[60] = 3.750000000000023e-02*G0_0_0 + 3.749999999999988e-02*G0_0_1 - 6.749999999999989e-01*G0_1_0 - 6.749999999999988e-01*G0_1_1;

348

block[61] = -3.749999999999945e-02*G0_0_0;

349

block[62] = 2.999999999999999e-01*G0_0_1 + 3.375000000000001e-01*G0_1_1;

350

block[63] = 3.374999999999992e-01*G0_0_0 + 1.687500000000003e-01*G0_0_1 + 1.687499999999994e-01*G0_1_0;

351

block[64] = 3.374999999999999e-01*G0_0_0 + 1.687499999999994e-01*G0_0_1 + 1.687500000000004e-01*G0_1_0;

352

block[65] = -3.374999999999999e-01*G0_0_0 - 1.181249999999999e+00*G0_0_1 - 1.687500000000005e-01*G0_1_0 - 1.349999999999999e+00*G0_1_1;

353

block[66] = 1.687499999999999e+00*G0_0_0 + 8.437499999999993e-01*G0_0_1 + 8.437499999999993e-01*G0_1_0 + 1.687499999999998e+00*G0_1_1;

354

block[67] = 8.437499999999999e-01*G0_0_1 + 8.437499999999993e-01*G0_1_0;

355

block[68] = -1.687500000000006e-01*G0_0_1 - 1.687500000000005e-01*G0_1_0;

356

block[69] = -2.024999999999998e+00*G0_0_0 - 1.012499999999999e+00*G0_0_1 - 1.012499999999999e+00*G0_1_0;

357

block[70] = -6.749999999999997e-01*G0_0_0 - 6.749999999999995e-01*G0_0_1 + 3.749999999999985e-02*G0_1_0 + 3.749999999999964e-02*G0_1_1;

358

block[71] = 3.375000000000000e-01*G0_0_0 + 3.000000000000002e-01*G0_1_0;

359

block[72] = -3.749999999999996e-02*G0_1_1;

360

block[73] = 1.687499999999996e-01*G0_0_1 + 1.687499999999996e-01*G0_1_0 + 3.374999999999994e-01*G0_1_1;

361

block[74] = 1.687499999999996e-01*G0_0_1 + 1.687499999999998e-01*G0_1_0 + 3.374999999999996e-01*G0_1_1;

362

block[75] = -1.687500000000000e-01*G0_0_1 - 1.687499999999997e-01*G0_1_0;

363

block[76] = 8.437499999999993e-01*G0_0_1 + 8.437499999999999e-01*G0_1_0;

364

block[77] = 1.687499999999999e+00*G0_0_0 + 8.437500000000001e-01*G0_0_1 + 8.437500000000002e-01*G0_1_0 + 1.687500000000000e+00*G0_1_1;

365

block[78] = -1.350000000000000e+00*G0_0_0 - 1.687500000000001e-01*G0_0_1 - 1.181250000000000e+00*G0_1_0 - 3.375000000000000e-01*G0_1_1;

366

block[79] = -1.012499999999999e+00*G0_0_1 - 1.012500000000000e+00*G0_1_0 - 2.024999999999999e+00*G0_1_1;

367

block[80] = 3.375000000000001e-01*G0_0_0 + 3.375000000000001e-01*G0_0_1 + 3.750000000000021e-02*G0_1_0 + 3.750000000000050e-02*G0_1_1;

368

block[81] = -6.749999999999996e-01*G0_0_0 - 7.124999999999992e-01*G0_1_0;

369

block[82] = -3.750000000000027e-02*G0_1_1;

370

block[83] = -8.437499999999990e-01*G0_0_1 - 8.437499999999994e-01*G0_1_0 - 1.687499999999998e+00*G0_1_1;

371

block[84] = 1.687499999999996e-01*G0_0_1 + 1.687499999999992e-01*G0_1_0 + 3.374999999999991e-01*G0_1_1;

372

block[85] = -1.687499999999992e-01*G0_0_1 - 1.687499999999992e-01*G0_1_0;

373

block[86] = -1.687500000000005e-01*G0_0_1 - 1.687500000000006e-01*G0_1_0;

374

block[87] = -1.350000000000000e+00*G0_0_0 - 1.181250000000000e+00*G0_0_1 - 1.687500000000001e-01*G0_1_0 - 3.375000000000000e-01*G0_1_1;

375

block[88] = 1.687499999999999e+00*G0_0_0 + 8.437499999999996e-01*G0_0_1 + 8.437499999999993e-01*G0_1_0 + 1.687499999999998e+00*G0_1_1;

376

block[89] = 1.012500000000000e+00*G0_0_1 + 1.012500000000000e+00*G0_1_0;

377

block[90] = 0.000000000000000e+00;

378

block[91] = 0.000000000000000e+00;

379

block[92] = 0.000000000000000e+00;

380

block[93] = -2.024999999999997e+00*G0_0_0 - 1.012500000000000e+00*G0_0_1 - 1.012499999999998e+00*G0_1_0;

381

block[94] = -1.012499999999998e+00*G0_0_1 - 1.012499999999998e+00*G0_1_0 - 2.024999999999997e+00*G0_1_1;

382

block[95] = 1.012500000000000e+00*G0_0_1 + 1.012499999999998e+00*G0_1_0;

383

block[96] = -2.024999999999999e+00*G0_0_0 - 1.012499999999999e+00*G0_0_1 - 1.012499999999999e+00*G0_1_0;

384

block[97] = -1.012500000000000e+00*G0_0_1 - 1.012499999999999e+00*G0_1_0 - 2.024999999999999e+00*G0_1_1;

385

block[98] = 1.012500000000000e+00*G0_0_1 + 1.012500000000000e+00*G0_1_0;

386

block[99] = 4.049999999999995e+00*G0_0_0 + 2.024999999999998e+00*G0_0_1 + 2.024999999999997e+00*G0_1_0 + 4.049999999999995e+00*G0_1_1;

387

}

388

389

};

390

391

/// This class contains the form to be evaluated, including

392

/// contributions from the interior and boundary of the domain.

393

394

class LinearForm : public dolfin::LinearForm

395

{

396

public:

397

398

class TestElement : public dolfin::FiniteElement

399

{

400

public:

401

402

TestElement() : dolfin::FiniteElement(), tensordims(0), subelements(0)

403

{

404

// Element is scalar, don't need to initialize tensordims

405

406

// Element is simple, don't need to initialize subelements

407

}

408

409

~TestElement()

410

{

411

if ( tensordims ) delete [] tensordims;

412

if ( subelements )

413

{

414

for (unsigned int i = 0; i < elementdim(); i++)

415

delete subelements[i];

416

delete [] subelements;

417

}

418

}

419

420

inline unsigned int spacedim() const

421

{

422

return 10;

423

}

424

425

inline unsigned int shapedim() const

426

{

427

return 2;

428

}

429

430

inline unsigned int tensordim(unsigned int i) const

431

{

432

dolfin_error("Element is scalar.");

433

return 0;

434

}

435

436

inline unsigned int elementdim() const

437

{

438

return 1;

439

}

440

441

inline unsigned int rank() const

442

{

443

return 0;

444

}

445

446

void nodemap(int nodes[], const Cell& cell, const Mesh& mesh) const

447

{

448

static unsigned int edge_reordering_0[2][2] = {{0, 1}, {1, 0}};

449

nodes[0] = cell.vertexID(0);

450

nodes[1] = cell.vertexID(1);

451

nodes[2] = cell.vertexID(2);

452

int alignment = cell.edgeAlignment(0);

453

int offset = mesh.numVertices();

454

nodes[3] = offset + 2*cell.edgeID(0) + edge_reordering_0[alignment][0];

455

nodes[4] = offset + 2*cell.edgeID(0) + edge_reordering_0[alignment][1];

456

alignment = cell.edgeAlignment(1);

457

nodes[5] = offset + 2*cell.edgeID(1) + edge_reordering_0[alignment][0];

458

nodes[6] = offset + 2*cell.edgeID(1) + edge_reordering_0[alignment][1];

459

alignment = cell.edgeAlignment(2);

460

nodes[7] = offset + 2*cell.edgeID(2) + edge_reordering_0[alignment][0];

461

nodes[8] = offset + 2*cell.edgeID(2) + edge_reordering_0[alignment][1];

462

offset = offset + 2*mesh.numEdges();

463

nodes[9] = offset + cell.id();

464

}

465

466

void pointmap(Point points[], unsigned int components[], const AffineMap& map) const

467

{

468

points[0] = map(0.000000000000000e+00, 0.000000000000000e+00);

469

points[1] = map(1.000000000000000e+00, 0.000000000000000e+00);

470

points[2] = map(0.000000000000000e+00, 1.000000000000000e+00);

471

points[3] = map(6.666666666666667e-01, 3.333333333333333e-01);

472

points[4] = map(3.333333333333334e-01, 6.666666666666666e-01);

473

points[5] = map(0.000000000000000e+00, 6.666666666666667e-01);

474

points[6] = map(0.000000000000000e+00, 3.333333333333334e-01);

475

points[7] = map(3.333333333333333e-01, 0.000000000000000e+00);

476

points[8] = map(6.666666666666666e-01, 0.000000000000000e+00);

477

points[9] = map(3.333333333333333e-01, 3.333333333333333e-01);

478

components[0] = 0;

479

components[1] = 0;

480

components[2] = 0;

481

components[3] = 0;

482

components[4] = 0;

483

components[5] = 0;

484

components[6] = 0;

485

components[7] = 0;

486

components[8] = 0;

487

components[9] = 0;

488

}

489

490

void vertexeval(real values[], unsigned int vertex, const real x[], const Mesh& mesh) const

491

{

492

// FIXME: Temporary fix for Lagrange elements

493

values[0] = x[vertex];

494

}

495

496

const FiniteElement& operator[] (unsigned int i) const

497

{

498

return *this;

499

}

500

501

FiniteElement& operator[] (unsigned int i)

502

{

503

return *this;

504

}

505

506

FiniteElementSpec spec() const

507

{

508

FiniteElementSpec s("Lagrange", "triangle", 3);

509

return s;

510

}

511

512

private:

513

514

unsigned int* tensordims;

515

FiniteElement** subelements;

516

517

};

518

519

class FunctionElement_0 : public dolfin::FiniteElement

520

{

521

public:

522

523

FunctionElement_0() : dolfin::FiniteElement(), tensordims(0), subelements(0)

524

{

525

// Element is scalar, don't need to initialize tensordims

526

527

// Element is simple, don't need to initialize subelements

528

}

529

530

~FunctionElement_0()

531

{

532

if ( tensordims ) delete [] tensordims;

533

if ( subelements )

534

{

535

for (unsigned int i = 0; i < elementdim(); i++)

536

delete subelements[i];

537

delete [] subelements;

538

}

539

}

540

541

inline unsigned int spacedim() const

542

{

543

return 10;

544

}

545

546

inline unsigned int shapedim() const

547

{

548

return 2;

549

}

550

551

inline unsigned int tensordim(unsigned int i) const

552

{

553

dolfin_error("Element is scalar.");

554

return 0;

555

}

556

557

inline unsigned int elementdim() const

558

{

559

return 1;

560

}

561

562

inline unsigned int rank() const

563

{

564

return 0;

565

}

566

567

void nodemap(int nodes[], const Cell& cell, const Mesh& mesh) const

568

{

569

static unsigned int edge_reordering_0[2][2] = {{0, 1}, {1, 0}};

570

nodes[0] = cell.vertexID(0);

571

nodes[1] = cell.vertexID(1);

572

nodes[2] = cell.vertexID(2);

573

int alignment = cell.edgeAlignment(0);

574

int offset = mesh.numVertices();

575

nodes[3] = offset + 2*cell.edgeID(0) + edge_reordering_0[alignment][0];

576

nodes[4] = offset + 2*cell.edgeID(0) + edge_reordering_0[alignment][1];

577

alignment = cell.edgeAlignment(1);

578

nodes[5] = offset + 2*cell.edgeID(1) + edge_reordering_0[alignment][0];

579

nodes[6] = offset + 2*cell.edgeID(1) + edge_reordering_0[alignment][1];

580

alignment = cell.edgeAlignment(2);

581

nodes[7] = offset + 2*cell.edgeID(2) + edge_reordering_0[alignment][0];

582

nodes[8] = offset + 2*cell.edgeID(2) + edge_reordering_0[alignment][1];

583

offset = offset + 2*mesh.numEdges();

584

nodes[9] = offset + cell.id();

585

}

586

587

void pointmap(Point points[], unsigned int components[], const AffineMap& map) const

588

{

589

points[0] = map(0.000000000000000e+00, 0.000000000000000e+00);

590

points[1] = map(1.000000000000000e+00, 0.000000000000000e+00);

591

points[2] = map(0.000000000000000e+00, 1.000000000000000e+00);

592

points[3] = map(6.666666666666667e-01, 3.333333333333333e-01);

593

points[4] = map(3.333333333333334e-01, 6.666666666666666e-01);

594

points[5] = map(0.000000000000000e+00, 6.666666666666667e-01);

595

points[6] = map(0.000000000000000e+00, 3.333333333333334e-01);

596

points[7] = map(3.333333333333333e-01, 0.000000000000000e+00);

597

points[8] = map(6.666666666666666e-01, 0.000000000000000e+00);

598

points[9] = map(3.333333333333333e-01, 3.333333333333333e-01);

599

components[0] = 0;

600

components[1] = 0;

601

components[2] = 0;

602

components[3] = 0;

603

components[4] = 0;

604

components[5] = 0;

605

components[6] = 0;

606

components[7] = 0;

607

components[8] = 0;

608

components[9] = 0;

609

}

610

611

void vertexeval(real values[], unsigned int vertex, const real x[], const Mesh& mesh) const

612

{

613

// FIXME: Temporary fix for Lagrange elements

614

values[0] = x[vertex];

615

}

616

617

const FiniteElement& operator[] (unsigned int i) const

618

{

619

return *this;

620

}

621

622

FiniteElement& operator[] (unsigned int i)

623

{

624

return *this;

625

}

626

627

FiniteElementSpec spec() const

628

{

629

FiniteElementSpec s("Lagrange", "triangle", 3);

630

return s;

631

}

632

633

private:

634

635

unsigned int* tensordims;

636

FiniteElement** subelements;

637

638

};

639

640

LinearForm(Function& w0) : dolfin::LinearForm(1)

641

{

642

// Create finite element for test space

643

_test = new TestElement();

644

645

// Add functions

646

add(w0, new FunctionElement_0());

647

}

648

649

void eval(real block[], const AffineMap& map) const

650

{

651

// Compute coefficients

652

const real c0_0 = c[0][0];

653

const real c0_1 = c[0][1];

654

const real c0_2 = c[0][2];

655

const real c0_3 = c[0][3];

656

const real c0_4 = c[0][4];

657

const real c0_5 = c[0][5];

658

const real c0_6 = c[0][6];

659

const real c0_7 = c[0][7];

660

const real c0_8 = c[0][8];

661

const real c0_9 = c[0][9];

662

663

// Compute geometry tensors

664

const real G0_0 = map.det*c0_0;

665

const real G0_1 = map.det*c0_1;

666

const real G0_2 = map.det*c0_2;

667

const real G0_3 = map.det*c0_3;

668

const real G0_4 = map.det*c0_4;

669

const real G0_5 = map.det*c0_5;

670

const real G0_6 = map.det*c0_6;

671

const real G0_7 = map.det*c0_7;

672

const real G0_8 = map.det*c0_8;

673

const real G0_9 = map.det*c0_9;

674

675

// Compute element tensor

676

block[0] = 5.654761904761898e-03*G0_0 + 8.184523809523806e-04*G0_1 + 8.184523809523796e-04*G0_2 + 2.008928571428569e-03*G0_3 + 2.008928571428570e-03*G0_4 + 1.339285714285713e-03*G0_6 + 1.339285714285708e-03*G0_7 + 2.678571428571429e-03*G0_9;

677

block[1] = 8.184523809523809e-04*G0_0 + 5.654761904761895e-03*G0_1 + 8.184523809523795e-04*G0_2 + 1.339285714285700e-03*G0_3 + 2.008928571428572e-03*G0_5 + 2.008928571428566e-03*G0_6 + 1.339285714285706e-03*G0_8 + 2.678571428571416e-03*G0_9;

678

block[2] = 8.184523809523796e-04*G0_0 + 8.184523809523796e-04*G0_1 + 5.654761904761899e-03*G0_2 + 1.339285714285713e-03*G0_4 + 1.339285714285709e-03*G0_5 + 2.008928571428571e-03*G0_7 + 2.008928571428569e-03*G0_8 + 2.678571428571422e-03*G0_9;

679

block[3] = 2.008928571428569e-03*G0_0 + 1.339285714285700e-03*G0_1 + 4.017857142857136e-02*G0_3 - 1.406249999999998e-02*G0_4 - 1.004464285714285e-02*G0_5 - 4.017857142857147e-03*G0_6 - 1.004464285714285e-02*G0_7 + 2.008928571428570e-02*G0_8 + 1.205357142857144e-02*G0_9;

680

block[4] = 2.008928571428570e-03*G0_0 + 1.339285714285713e-03*G0_2 - 1.406249999999998e-02*G0_3 + 4.017857142857139e-02*G0_4 + 2.008928571428570e-02*G0_5 - 1.004464285714285e-02*G0_6 - 4.017857142857136e-03*G0_7 - 1.004464285714285e-02*G0_8 + 1.205357142857142e-02*G0_9;

681

block[5] = 2.008928571428572e-03*G0_1 + 1.339285714285710e-03*G0_2 - 1.004464285714285e-02*G0_3 + 2.008928571428570e-02*G0_4 + 4.017857142857139e-02*G0_5 - 1.406249999999999e-02*G0_6 - 1.004464285714284e-02*G0_7 - 4.017857142857138e-03*G0_8 + 1.205357142857142e-02*G0_9;

682

block[6] = 1.339285714285713e-03*G0_0 + 2.008928571428566e-03*G0_1 - 4.017857142857146e-03*G0_3 - 1.004464285714285e-02*G0_4 - 1.406249999999999e-02*G0_5 + 4.017857142857138e-02*G0_6 + 2.008928571428570e-02*G0_7 - 1.004464285714285e-02*G0_8 + 1.205357142857141e-02*G0_9;

683

block[7] = 1.339285714285707e-03*G0_0 + 2.008928571428571e-03*G0_2 - 1.004464285714285e-02*G0_3 - 4.017857142857136e-03*G0_4 - 1.004464285714284e-02*G0_5 + 2.008928571428570e-02*G0_6 + 4.017857142857141e-02*G0_7 - 1.406249999999999e-02*G0_8 + 1.205357142857142e-02*G0_9;

684

block[8] = 1.339285714285706e-03*G0_1 + 2.008928571428569e-03*G0_2 + 2.008928571428570e-02*G0_3 - 1.004464285714285e-02*G0_4 - 4.017857142857138e-03*G0_5 - 1.004464285714285e-02*G0_6 - 1.406249999999999e-02*G0_7 + 4.017857142857140e-02*G0_8 + 1.205357142857143e-02*G0_9;

685

block[9] = 2.678571428571429e-03*G0_0 + 2.678571428571416e-03*G0_1 + 2.678571428571422e-03*G0_2 + 1.205357142857144e-02*G0_3 + 1.205357142857142e-02*G0_4 + 1.205357142857142e-02*G0_5 + 1.205357142857141e-02*G0_6 + 1.205357142857142e-02*G0_7 + 1.205357142857143e-02*G0_8 + 1.446428571428571e-01*G0_9;

686

}

687

688

};

689

690

} }

691

692

#endif

Older »