JuliaInv / JOcTree / 96

64%

push

travis-ci

<a href="https://github.com/JuliaInv/JOcTree/commit/<a class=hub.com/JuliaInv/JOcTree/commit/<a class="double-link" href="https://git"><a class=hub.com/JuliaInv/JOcTree/commit/<a class="double-link" href="https://git"><a class=hub.com/JuliaInv/JOcTree/commit/<a class="double-link" href="https://git"><a class=hub.com/JuliaInv/JOcTree/commit/<a class="double-link" href="https://git"><a class=hub.com/JuliaInv/JOcTree/commit/<a class="double-link" href="https://git"><a class=hub.com/JuliaInv/JOcTree/commit/<a class="double-link" href="https://git"><a class=hub.com/JuliaInv/JOcTree/commit/f62d822d4bb044dc9b6fc856e73d19c4c266ea5f">f62d822d4<a href="https://github.com/JuliaInv/JOcTree/commit/f62d822d4bb044dc9b6fc856e73d19c4c266ea5f">&lt;a href=&quot;https://github.com/JuliaInv/JOcTree/commit/f62d822d4bb044dc9b6fc856e73d19c4c266ea5f&quot;&gt;&amp;quot;&amp;gt;&amp;amp;lt;a href=&amp;amp;quot;https://github.com/JuliaInv/JOcTree/commit/&amp;lt;/a&amp;gt;&amp;lt;a class=&amp;quot;double-link&amp;quot; href=&amp;quot;https://github.com/JuliaInv/JOcTree/commit/&amp;amp;lt;a class=&amp;amp;quot;double-link&amp;amp;quot; href=&amp;amp;quot;https://git&amp;quot;&amp;gt;&amp;amp;lt;a class=&amp;lt;/a&amp;gt;hub.com/JuliaInv/JOcTree/commit/&amp;amp;lt;a class=&amp;amp;quot;double-link&amp;amp;quot; href=&amp;amp;quot;https://git&amp;quot;&amp;gt;&amp;amp;lt;a class=&amp;lt;/a&amp;gt;hub.com/JuliaInv/JOcTree/commit/f62d822d4bb044dc9b6fc856e73d19c4c266ea5f&amp;quot;&amp;gt;f62d822d4&amp;lt;/a&amp;gt;&amp;lt;a href=&amp;quot;https://github.com/JuliaInv/JOcTree/commit/f62d822d4bb044dc9b6fc856e73d19c4c266ea5f&amp;quot;&amp;gt;&amp;amp;lt;a href=&amp;amp;quot;https://github.com/JuliaInv/JOcTree/commit/f62d822d4bb044dc9b6fc856e73d19c4c266ea5f&amp;amp;quot;&amp;amp;gt;&amp;amp;amp;lt;a href=&amp;amp;amp;quot;https://github.com/JuliaInv/JOcTree/commit/f62d822d4bb044dc9b6fc856e73d19c4c266ea5f&amp;amp;quot;&amp;amp;gt;&amp;amp;amp;quot;&amp;amp;amp;gt;Unify mass matrix integration for nodes, edges and faces (#25)

* Make getEdgeMassMatrix faster and consistent

Commit &amp;amp;amp;lt;/a&amp;amp;amp;gt;&amp;amp;amp;lt;a class=&amp;amp;amp;quot;double-link&amp;amp;amp;quot; href=&amp;amp;amp;quot;https://github.com/JuliaInv/JOcTree/commit/&amp;amp;lt;/a&amp;amp;gt;&amp;amp;lt;a class=&amp;amp;quot;double-link&amp;amp;quot; href=&amp;amp;quot;https://github.com/JuliaInv/JOcTree/commit/&lt;a class=&quot;double-link&quot; href=&quot;https://github.com/JuliaInv/JOcTree/commit/403ce28feca25f46f2bb8ebfd6aa8593806d548e&quot;&gt;403ce28fe&lt;/a&gt;&amp;amp;quot;&amp;amp;gt;403ce28fe&amp;amp;lt;/a&amp;amp;gt;&amp;amp;quot;&amp;amp;gt;403ce28fe&amp;amp;lt;/a&amp;amp;gt;&amp;amp;lt;a href=&amp;amp;quot;https://github.com/JuliaInv/JOcTree/commit/f62d822d4bb044dc9b6fc856e73d19c4c266ea5f&amp;amp;quot;&amp;amp;gt; reverted edge mass matrix
integration for an isotropic coefficient to edge to cell center averaging. This
gave inconsistent mass matrices with the integration for an anisotropic coefficient.

Major changes: edge mass matrix assembly
 - Use consistent integration method for isotropic and anisotropic coefficients.
 - Compute everything required for mass matrix integration once and store in new
   container (struct) EdgeMassMatrix. This accelerates the assembly for subsequent
   calls as well as the mass matrix derivatives significantly.

Minor changes: cleanup of OcTreeMeshFV
 - Store number of nodes.
 - Store nodal numbering and cell numbering.
 - Remove storage for transposes of Ae and Pe.
 - The methods clear!(::Array) and clear!(::SparseMatrixCSC) implemented in jInv.Utils
   are a misnomer. Instead of changing the contents, as implied by the exclamation mark,
   a new, empty object is returned. To make this evident, clear!(::OcTree&lt;/a&gt;MeshFV) now
   calls the constructors for empty objects directly and avoids calling clear!.
 - Make display(::OcTreeMeshFV) work for empty mesh (after calling clear!(::OcTreeMeshFV)).

* Cleanup clear!(S::S</a>parseArray3D)

Continuation of </a><a class="double-link" href="https://github.com/JuliaInv/JOcTree/commit/&lt;a class=&quot;double-link&quot; href=&quot;https://git">&lt;a class=</a>hub.com/JuliaInv/JOcTree/commit/<a class="double-link" href="https://github.com/JuliaInv/JOcTree/commit/&lt;a class=&quot;double-link" href="https://git"><a class=hub.com/JuliaInv/JOcTree/commit/<a class=hub.com/JuliaInv/JOcTree/commit/3048a0fbee5dd345892aa442a6a9b6b2c383d427">3048a0fbe">3048a0fbe</a><a href="https://github.com/JuliaInv/JOcTree/commit/f62d822d4bb044dc9b6fc856e73d19c4c266ea5f">&quot;&gt;3048a0fbe&lt;/a&gt;&lt;a href=&quot;https://github.com/JuliaInv/JOcTree/commit/f62d822d4bb044dc9b6fc856e73d19c4c266ea5f&quot;&gt;

* Improve performance of getEdgeMassMatrix

* Fix edge mass matrix derivative for complex fields

* Export method getVolumeVector

* Add test for consistency of anisotropic mass matrices

* Add documentation

- Location of code inside file getEdgeMassMatrix.jl has changed but
  not the contents.

* 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 3048a0fbee5dd345892aa442a6a9b6b2c383d427
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.

* Use integer arithmetics in findNonRegularBlocks

* Fix automated testing for optional package HDF5

* Fix automated testing for optional package HDF5&lt;/a&gt;</a>

1871 of 2906 relevant lines covered (64.38%)

37575.73 hits per line