subroutine, public mesh_torus(msh, h, r, center, radius, center_2, radius_2)
Mesh the unit torus.
subroutine, public fefunc_visualise(uh, Xh, type, title, fic, pause, verb)
Visualise a single finite element function (with gmsh)
real(rp) function, public l2dist_scalfunc_scalfe(uh, u, Xh, qdm, verb)
Given the vertexes of an edge and of a triangle.
integer, parameter, public fe_p3_1d
integer, parameter, public fe_p1_3d
DERIVED TYPE fe: finite element description
subroutine, public mesh_ball(msh, h, r, center, radius)
Mesh a ball of given center and radius Default = unit ball.
Module for feMeth outputs:
subroutine, public scalfe_write(Xh, uh, fic, fmt, name, time, idx, verb)
Write a scakar finite element function to file fic with given file format.
subroutine, public sourceterm_scalfe_lift(V, Xh, qdm, g, lift, k_ER, b, Db, verb)
Scalar source term with a lift
integer, parameter, public fe_p2_2d
subroutine, public mesh_ring(msh, h, r, center, radius, center_2, radius_2)
Mesh ring of given centers and radius Default = unit.
integer, parameter, public fe_p1_1d
DERIVED TYPE cell: geometrical description of mesh cells
integer, parameter, public cell_trg_2
Quadratic triangle.
Definition of for a finite element space.
DERIVED TYPE quadMesh: integration methods on meshes
real(rp) function, public l2dist_vectfunc_grad_scalfe_lift(uh, phi, Xh, qdm, lift, k_ER, b, Db, verb)
Case 1. Returns
Tools to build meshes using GMSH.
real(rp) function, public l2dist_vectfunc_grad_scalfe(uh, phi, Xh, qdm, verb)
Returns .
subroutine, public interpolate_scalfunc(uh, Xh, u, verb)
Finite element IInterpolant of a scalar function .
subroutine, public mesh_tag_cells(tag, m, dim, verb, f)
Tag the cells of mesh 'm': tag(cl) = .TRUE. iif:
real(rp) function, public integrate_grad_scalfe_lift(uh, phi, E, Xh, qdm, lift, k_ER, b, Db, verb)
Case 1. Returns ,
integer, parameter, public fe_p1_2d
subroutine, public mesh_write(m, fic, fmt, verb)
Write a mesh in a file with given file format
integer, parameter, public quad_gauss_trg_46
real(rp) function, dimension(n), public integrate_vectfunc_lift(func, n, qdm, lift, k_ER, b, Db, verb)
integer function, public pk_method_index(dim, deg)
Pk finite element method index.
subroutine, public vectfe_write(Xh, uh, fic, fmt, name, time, idx, verb)
Write a vector finite element function to file 'fic' with file format 'fmt'
integer, parameter, public quad_gauss_tet_1
subroutine, public gmsh_std_mesh(msh, N, r, geo)
Mesh cosructor: simple geometry used for integration tests.
DERIVED TYPE quad: quadrature rules description
integer, parameter, public quad_gauss_trg_16
integer, parameter, public fe_p0_2d
subroutine, public massmat_scalfe(mass, rho, Xh, qdm, pattern, verb, OMP)
Mass matrix:
real(rp) function, dimension(n, n), public l2grammatrix_vectfuncset_lift(f, n, qdm, lift, K_er, b, Db, verb)
integer, parameter, public cell_trg_3
Quadratic triangle.
Scalar finite element spaces.
real(rp) function, public scalar_1(x)
The function .
integer, parameter, public cell_vtx
Vertex.
integer, parameter, public quad_gauss_tet_15
integer, parameter, public lift_ortho_proj
real(rp) function, public integrate_grad_scalfe(uh, E, Xh, qdm, verb)
Returns .
integer, parameter, public quad_gauss_trg_19
subroutine, public linearform_grad_scalfe(B, Xh, v, Yh, qdm, func_f, verb)
Assembles the matrix of the linear form
assemble the output (finalise construction)
subroutine, public scan_xh_boundary(interior_cell_type, interior_cell_fe_type, Xh, verb)
Informations on boundary cells.
integer, parameter, public cell_trg
Triangle.
real(rp) function, public l2dist_scalfuncspace_scalfe_lift(uh, func_f, dim, Xh, qdm, lift, K_er, b, Db, verb)
Returns : L2-Distance betwween:
real(rp) function, dimension(n, n), public l2grammatrix_scalfuncset(f, n, qdm, verb)
Compute the L -product Gram matrix of a set of scalar functions:
real(rp) function, dimension(n), public l2prod_scalfuncset_scalfe_lift(uh, func_f, N, Xh, qdm, lift, k_ER, b, Db, verb)
Returns with
integer, parameter, public quad_gauss_trg_6
Vector finite element spaces.
subroutine, public stiffmat_elasticity(stiff, lambda, mu, Xh, qdm, pattern, verb, OMP)
Stiffness matrix in linear elasticity:
real(rp) function, public isotropic_1(x, u, v)
Euclidian metric on : This is the isotropic homogeneous metric equal to the usual scalar product ...
TOP-LEVEL MODULE for the 'feMeth' project.
integer, parameter, public quad_gauss_trg_3
real(rp) function, public l2dist_scalfunc_scalfe_lift(uh, u, Xh, qdm, lift, k_ER, b, Db, verb)
Returns .
integer, parameter, public fe_p0_1d
Finite element IInterpolant of a vector function
Derived type feSpace: finite element space on a mesh.
integer, parameter, public fe_p2_1d
Vector Dirichlet boundary condition on
subroutine, public mesh_disk(msh, h, r, center, radius)
Mesh a disk of given center and radius Default = unit disk.
real(rp) function, dimension(n), public l2prod_scalfuncset_scalfe(uh, func, n, Xh, qdm, verb)
Compute the -product between a finite element function and a set of scalar functions.
integer, parameter, public fe_p0_3d
integer, parameter, public quad_gauss_edg_3
subroutine, public mmg2d_remesh_levelset(m, val, verb, hausd, hMax)
Remesh respectively to a levelSet.
The type quadMesh defines integration methods on meshes.
integer, parameter, public fe_p2_3d
subroutine, public massmat_vectfe(mass, rho, Xh, qdm, pattern, verb, OMP)
Mass matrix:
subroutine, public fespace_write(Xh, fic, fmt, verb)
Write the finite element space 'Xh' in a file with given file format
real(rp) function, dimension(n, n), public l2grammatrix_scalfuncset_lift(f, N, qdm, lift, k_ER, b, Db, verb)
Similar to L2GramMatrix_scalFuncSet using a lift.
integer, parameter, public cell_tet_2
Quadratic tetrahedron.
subroutine, public tag_boundary_cell(msh, verb)
Let a mesh 'msh' with domain let its boundary.
integer, parameter, public quad_gauss_edg_4
integer, parameter, public lift_elliott_etal
subroutine, public dirichlet_cond_scalfe(K, rhs, Xh, g, rho, verb, f)
Scalar Dirichlet boundary condition (using penalisation) on
Vector source term: matrix of linear forms
subroutine, public mesh_jordan_curve_interior(msh, N, r, g, g_inv, verb)
See mesh_Jordan_curve_0, case dim =2.
finite element operator sparsity pattern
real(rp) function, public l2dist_vectfuncspace_grad_scalfe(uh, gard_func, dim, Xh, qdm, verb)
Returns : L2-Distance betwween:
integer function, public pk_max_implemented_degree(dim)
Maximal implemented degree for the Pk finite elements in dimension 'dim'.
integer, parameter, public quad_gauss_edg_2
set the output parameters
integer, parameter, public quad_gauss_trg_12
DERIVED TYPE feSpace: finite element spaces
Computation of integrals using meshes.
integer, parameter, public fe_p3_2d
subroutine, public mesh_box(msh, h, r, dx, dy, dz)
Mesh a box.
integer, parameter, public lift_demlow
integer, parameter, public fe_0d
Transformation between comutational and physical domains.
integer, parameter, public lift_modified_elliott_etal
integer, parameter, public cell_edg_2
Quadratic edge.
subroutine, public mesh_circle(msh, h, r, center, radius)
Mesh a circle of given center and radius Default = unit circle.
integer, parameter, public cell_edg_3
Cubic edge.
integer, parameter, public cell_tet
Tetrahedron.
integer, parameter, public cell_edg
Edge (line segment)
real(rp) function, dimension(n), public l2prod_vectfuncset_grad_scalfe(uh, func, n, Xh, qdm, verb)
Compute the -product between the gradient of a finite element function and a set of vector functions...
real(rp) function, public integrate_scalfe_lift(uh, u, E, Xh, qdm, lift, k_ER, b, Db, verb)
Returns .
subroutine, public mesh_disk_with_inclusion(msh, r, c1, r1, h1, c2, r2, h2)
Mesg the domain with an inclusion .
real(rp) function, public l2dist_scalfuncspace_scalfe(uh, func, dim, Xh, qdm, verb)
Returns : L2-Distance betwween:
subroutine, public cellvals_write(m, vals, fic, fmt, name, time, idx, verb)
Write scakar values on the mesh cells ( data)
subroutine, public mesh_describe(m)
Display a description of the mesh 'm'. Including the cells, edges and faces types.
real(rp) function, dimension(n), public l2prod_vectfuncset_grad_scalfe_lift(uh, func, n, Xh, qdm, lift, k_ER, b, Db, verb)
real(rp) function, public l2dist_vectfuncspace_grad_scalfe_lift(uh, gard_func, dim, Xh, qdm, lift, K_er, b, Db, verb)
subroutine, public mesh_line(msh, h, r, P1, P2)
Mesh a line segment.
real(rp) function, dimension(n, n), public l2grammatrix_vectfuncset(f, n, qdm, verb)
Compute the L -product Gram matrix of a set of vector functions:
real(rp) function, public l2norm_scalfe(uh, Xh, qdm, verb)
Returns .
real(rp) function, public integrate_scalfunc(func, qdm, verb)
Returns
DERIVED TYPE mesh: definition of meshes
integer, parameter, public quad_gauss_tet_4
subroutine, public mesh_rectangle(msh, h, r, dx, dy)
Mesh a rectangle.
real(rp) function, public integrate_scalfunc_lift(u, qdm, lift, k_ER, b, Db, verb)
Returns .
integer, parameter, public quad_gauss_edg_1
subroutine, public mesh_sphere(msh, h, r, center, radius)
Mesh a sphere of given center and radius Default = unit sphere.
real(rp) function, public integrate_dnu_scalfe(uh, E, Xh, qdm, verb)
Returns .
integer, parameter, public quad_gauss_trg_1
integer, parameter, public quad_0d
subroutine, public mesh_jordan_curve(msh, N, r, g, g_inv, verb)
See mesh_Jordan_curve_0, case dim =1.
