JuliaInv / JOcTree / 86 / 1

64%

master: 64%

Build Type

push

travis-ci

Committed by Christoph Schwarzbach

Commit Message

<a href="https://github.com/JuliaInv/JOcTree/commit/<a class=hub.com/JuliaInv/JOcTree/commit/8dbf4c0a04402d08102accf9c6539694c36601f1">8dbf4c0a0<a href="https://github.com/JuliaInv/JOcTree/commit/8dbf4c0a04402d08102accf9c6539694c36601f1">&quot;&gt;Unify mass matrix integration for nodes, edges and faces

The same 8-point quadrature rule is applied for the integration
of nodal, edge and face mass matrices. This commit extends the new
implementation for edge mass matrices from commit </a><a class="double-link" href="https://github.com/JuliaInv/JOcTree/commit/<a class="double-link" href="https://github.com/JuliaInv/JOcTree/commit/3048a0fbee5dd345892aa442a6a9b6b2c383d427">3048a0fbe</a>">3048a0fbe</a><a href="https://github.com/JuliaInv/JOcTree/commit/8dbf4c0a04402d08102accf9c6539694c36601f1">
to nodal and face mass matrices.

The new implementation of getNodalMassMatrix neglects hanging nodes. Nodal
mass matrices obtained with the 8-point quadrature differ from the previous
implementation (based on getNodalAverageMatrix) unless the OcTree is a regular
mesh.

Changes

- Use consistent integration method for all types of mass matrices and
  isotropic and anisotropic coefficients.

- Compute everything required for mass matrix integration once and store in new
  container (struct) MassMatrix. This accelerates the assembly for subsequent
  calls as well as the mass matrix derivatives significantly.

- Replace structure EdgeMassMatrix and its default constructor by
  MassMatrix and move to OcTreeMeshFV.jl.

- Move method merge!(a::Vector{Int64}, b::Vector{Int64}) to Utils.jl.

- Add methods getdNodalMassMatrix, dNodalMassMatrixTimesVector and
  dNodalMassMatrixTrTimesVector.

- Add methods dFaceMassMatrixTimesVector and dFaceMassMatrixTrTimesVector.

- Add tests for new methods.

- Add tests for mass matrix integration which check that constant and
  linear coefficients are integrated exactly.

Coverage Stats

1869 of 2904 relevant lines covered (64.36%)

18901.0 hits per line