]> https://gitweb.dealii.org/ - dealii.git/log
dealii.git
3 years agoMerge pull request #12709 from peterrum/is_artificial_on_level
Wolfgang Bangerth [Fri, 27 Aug 2021 03:08:05 +0000 (21:08 -0600)]
Merge pull request #12709 from peterrum/is_artificial_on_level

Introduce CellAccessor::is_artificial_on_level()

3 years agoMerge pull request #12707 from peterrum/range-based-loop
Wolfgang Bangerth [Fri, 27 Aug 2021 03:07:24 +0000 (21:07 -0600)]
Merge pull request #12707 from peterrum/range-based-loop

Replace loop by range-based loop

3 years agoReplace loop by range-based loop 12707/head
Peter Munch [Wed, 25 Aug 2021 17:15:38 +0000 (19:15 +0200)]
Replace loop by range-based loop

3 years agoIntroduce CellAccessor::is_artificial_on_level() 12709/head
Peter Munch [Wed, 25 Aug 2021 17:42:46 +0000 (19:42 +0200)]
Introduce CellAccessor::is_artificial_on_level()

3 years agoMerge pull request #12710 from bangerth/identities-5
Martin Kronbichler [Thu, 26 Aug 2021 07:27:13 +0000 (09:27 +0200)]
Merge pull request #12710 from bangerth/identities-5

Short-circuit compute_hp_dof_identities() when possible.

3 years agoMerge pull request #12702 from peterrum/setup_constraints_edge_multiple_components
Martin Kronbichler [Thu, 26 Aug 2021 05:10:11 +0000 (07:10 +0200)]
Merge pull request #12702 from peterrum/setup_constraints_edge_multiple_components

Fix HangingNodes::setup_constraints() for multiple components

3 years agoMerge pull request #12684 from bangerth/tensor-array-view
David Wells [Wed, 25 Aug 2021 23:49:50 +0000 (19:49 -0400)]
Merge pull request #12684 from bangerth/tensor-array-view

Introduce ParticleAccessor::set_properties() taking a Tensor.

3 years agoShort-circuit compute_hp_dof_identities() when possible. 12710/head
Wolfgang Bangerth [Wed, 25 Aug 2021 19:09:45 +0000 (13:09 -0600)]
Short-circuit compute_hp_dof_identities() when possible.

3 years agoAdd a changelog entry. 12684/head
Wolfgang Bangerth [Fri, 20 Aug 2021 01:38:58 +0000 (19:38 -0600)]
Add a changelog entry.

3 years agoFix HangingNodes::setup_constraints() for multiple components 12702/head
Peter Munch [Tue, 24 Aug 2021 08:34:09 +0000 (10:34 +0200)]
Fix HangingNodes::setup_constraints() for multiple components

3 years agoMerge pull request #12706 from bangerth/identities-2
Marc Fehling [Wed, 25 Aug 2021 15:42:41 +0000 (09:42 -0600)]
Merge pull request #12706 from bangerth/identities-2

Minor changes to internal function arguments.

3 years agoMerge pull request #12704 from peterrum/pft_create_description_from_triangulation_res...
Peter Munch [Wed, 25 Aug 2021 05:14:39 +0000 (07:14 +0200)]
Merge pull request #12704 from peterrum/pft_create_description_from_triangulation_restructure

Restructure TriangulationDescription::Utilities::create_description_from_triangulation()

3 years agoUpdate source/dofs/dof_handler_policy.cc 12706/head
Wolfgang Bangerth [Wed, 25 Aug 2021 03:18:17 +0000 (21:18 -0600)]
Update source/dofs/dof_handler_policy.cc

Co-authored-by: Marc Fehling <mafehling.git@gmail.com>
3 years agoMinor changes to internal function arguments.
Wolfgang Bangerth [Tue, 24 Aug 2021 20:41:07 +0000 (14:41 -0600)]
Minor changes to internal function arguments.

3 years agoMerge pull request #12690 from bangerth/identities
Marc Fehling [Wed, 25 Aug 2021 01:34:44 +0000 (19:34 -0600)]
Merge pull request #12690 from bangerth/identities

Implement multiway DoF identities on vertices/lines/quads.

3 years agoAdd more tests. 12690/head
Wolfgang Bangerth [Tue, 24 Aug 2021 20:02:18 +0000 (14:02 -0600)]
Add more tests.

3 years agoAdd a changelog entry.
Wolfgang Bangerth [Fri, 20 Aug 2021 19:22:15 +0000 (13:22 -0600)]
Add a changelog entry.

3 years agoAdd a number of tests.
Wolfgang Bangerth [Fri, 20 Aug 2021 19:18:43 +0000 (13:18 -0600)]
Add a number of tests.

3 years agoRestrcuture TriangulationDescription::Utilities::create_description_from_triangulation() 12704/head
Peter Munch [Tue, 24 Aug 2021 19:37:31 +0000 (21:37 +0200)]
Restrcuture TriangulationDescription::Utilities::create_description_from_triangulation()

3 years agoImplement multiway DoF identities on vertices/lines/quads.
Wolfgang Bangerth [Fri, 20 Aug 2021 19:18:22 +0000 (13:18 -0600)]
Implement multiway DoF identities on vertices/lines/quads.

3 years agoMerge pull request #12696 from peterrum/p_mg_any_level
Martin Kronbichler [Tue, 24 Aug 2021 08:40:45 +0000 (10:40 +0200)]
Merge pull request #12696 from peterrum/p_mg_any_level

Extend MGTwoLevelTransfer::reinit_polynomial_transfer() so that it works on any refinement level without HN

3 years agoMerge pull request #12700 from gfcas/typename_number
Martin Kronbichler [Mon, 23 Aug 2021 20:03:59 +0000 (22:03 +0200)]
Merge pull request #12700 from gfcas/typename_number

Fix typename number

3 years agoMerge pull request #12698 from gfcas/dimension_mismatch
Daniel Arndt [Mon, 23 Aug 2021 18:39:47 +0000 (14:39 -0400)]
Merge pull request #12698 from gfcas/dimension_mismatch

3 years agoUpdate include/deal.II/numerics/vector_tools_rhs.templates.h 12700/head
Fabian Castelli [Mon, 23 Aug 2021 18:08:04 +0000 (20:08 +0200)]
Update include/deal.II/numerics/vector_tools_rhs.templates.h

Co-authored-by: Daniel Arndt <arndtd@ornl.gov>
3 years agoUpdate include/deal.II/numerics/vector_tools_rhs.templates.h
Fabian Castelli [Mon, 23 Aug 2021 18:07:56 +0000 (20:07 +0200)]
Update include/deal.II/numerics/vector_tools_rhs.templates.h

Co-authored-by: Daniel Arndt <arndtd@ornl.gov>
3 years agoFix typename number
Fabian Castelli [Mon, 23 Aug 2021 17:27:40 +0000 (19:27 +0200)]
Fix typename number

3 years agoMerge pull request #12676 from peterrum/pft_manifold_2
Martin Kronbichler [Mon, 23 Aug 2021 07:49:17 +0000 (09:49 +0200)]
Merge pull request #12676 from peterrum/pft_manifold_2

Remove special treatment of boundary cells

3 years agoIndent 12698/head
Fabian Castelli [Mon, 23 Aug 2021 07:33:31 +0000 (09:33 +0200)]
Indent

3 years agoMerge pull request #12685 from bangerth/array-view-2
Martin Kronbichler [Mon, 23 Aug 2021 07:26:42 +0000 (09:26 +0200)]
Merge pull request #12685 from bangerth/array-view-2

Avoid the use of ArrayView for Tensors in MPI::sum().

3 years agoMerge pull request #12634 from tjhei/mg_transfer_bug
Martin Kronbichler [Mon, 23 Aug 2021 07:25:58 +0000 (09:25 +0200)]
Merge pull request #12634 from tjhei/mg_transfer_bug

Fix MGTransferMatrixFree plain_copy bug

3 years agoMerge pull request #12697 from lethe-cfd/tolerance_subdivided_cylinder
Martin Kronbichler [Mon, 23 Aug 2021 07:22:10 +0000 (09:22 +0200)]
Merge pull request #12697 from lethe-cfd/tolerance_subdivided_cylinder

Add tolerance to subdivided_cylinder generator

3 years agoAssertDimension macro in examples
Fabian Castelli [Sat, 21 Aug 2021 15:23:26 +0000 (17:23 +0200)]
AssertDimension macro in examples

3 years agoAdd tolerance to subdivided_cylinder generator 12697/head
Bruno [Sun, 22 Aug 2021 19:16:47 +0000 (15:16 -0400)]
Add tolerance to subdivided_cylinder generator

3 years agoactually fix the problem 12634/head
Timo Heister [Fri, 20 Aug 2021 21:19:08 +0000 (17:19 -0400)]
actually fix the problem

3 years agoExtend MGTwoLevelTransfer::reinit_polynomial_transfer() so that it works on any refin... 12696/head
Peter Munch [Sat, 21 Aug 2021 12:16:18 +0000 (14:16 +0200)]
Extend MGTwoLevelTransfer::reinit_polynomial_transfer() so that it works on any refinement level without HN

3 years agoMerge pull request #12493 from peterrum/MGTwoLevelTransfer_ls
Peter Munch [Sat, 21 Aug 2021 11:15:33 +0000 (13:15 +0200)]
Merge pull request #12493 from peterrum/MGTwoLevelTransfer_ls

Enable MGTwoLevelTransfer::reinit_geometric_transfer() for refinement levels

3 years agoEnable MGTwoLevelTransfer::reinit_geometric_transfer() for refinment levels 12493/head
Peter Munch [Tue, 17 Aug 2021 07:19:51 +0000 (09:19 +0200)]
Enable MGTwoLevelTransfer::reinit_geometric_transfer() for refinment levels

3 years agoMerge pull request #12695 from bangerth/stacktrace
Martin Kronbichler [Sat, 21 Aug 2021 06:19:07 +0000 (08:19 +0200)]
Merge pull request #12695 from bangerth/stacktrace

Do not store unnecessary variables.

3 years agoRemove special treatment of boundary cells 12676/head
Peter Munch [Thu, 19 Aug 2021 06:49:05 +0000 (08:49 +0200)]
Remove special treatment of boundary cells

3 years agoDo not store unnecessary variables. 12695/head
Wolfgang Bangerth [Fri, 20 Aug 2021 23:11:54 +0000 (17:11 -0600)]
Do not store unnecessary variables.

3 years agoMerge pull request #12694 from bangerth/exception
Marc Fehling [Sat, 21 Aug 2021 02:24:55 +0000 (20:24 -0600)]
Merge pull request #12694 from bangerth/exception

Make sure that copying exceptions doesn't imply dropping the backtrace.

3 years agoAdd changelog entry. 12694/head
Wolfgang Bangerth [Fri, 20 Aug 2021 22:54:43 +0000 (16:54 -0600)]
Add changelog entry.

3 years agoMake sure that copying exceptions doesn't imply dropping the backtrace.
Wolfgang Bangerth [Fri, 20 Aug 2021 21:57:54 +0000 (15:57 -0600)]
Make sure that copying exceptions doesn't imply dropping the backtrace.

3 years agoAvoid the use of ArrayView for Tensors in MPI::sum(). 12685/head
Wolfgang Bangerth [Fri, 20 Aug 2021 02:19:55 +0000 (20:19 -0600)]
Avoid the use of ArrayView for Tensors in MPI::sum().

3 years agoMerge pull request #12687 from bangerth/sort
Martin Kronbichler [Fri, 20 Aug 2021 08:00:48 +0000 (10:00 +0200)]
Merge pull request #12687 from bangerth/sort

Sort member functions in .cc file.

3 years agoMerge pull request #12688 from marcfehling/fe-index
Martin Kronbichler [Fri, 20 Aug 2021 08:00:15 +0000 (10:00 +0200)]
Merge pull request #12688 from marcfehling/fe-index

Unify spelling 2.

3 years agoUnify spelling 2. 12688/head
Marc Fehling [Fri, 20 Aug 2021 02:42:38 +0000 (20:42 -0600)]
Unify spelling 2.

3 years agoMerge pull request #12683 from bangerth/doc-5
Marc Fehling [Fri, 20 Aug 2021 02:42:20 +0000 (20:42 -0600)]
Merge pull request #12683 from bangerth/doc-5

Unify spelling.

3 years agoSort member functions in .cc file. 12687/head
Wolfgang Bangerth [Fri, 20 Aug 2021 02:12:53 +0000 (20:12 -0600)]
Sort member functions in .cc file.

3 years agoMerge pull request #12681 from bangerth/1d
Marc Fehling [Fri, 20 Aug 2021 02:15:23 +0000 (20:15 -0600)]
Merge pull request #12681 from bangerth/1d

Allow for the default construction of 1d face iterators.

3 years agoMerge pull request #12682 from bangerth/doc-4
Marc Fehling [Fri, 20 Aug 2021 02:15:06 +0000 (20:15 -0600)]
Merge pull request #12682 from bangerth/doc-4

Group member functions in hp::FECollection for doxygen.

3 years agoMerge pull request #12680 from marcfehling/fix-75
Wolfgang Bangerth [Fri, 20 Aug 2021 02:14:17 +0000 (20:14 -0600)]
Merge pull request #12680 from marcfehling/fix-75

Prevent excess timer output.

3 years agoUse ParticleAccessor::set_properties() taking a Tensor.
Wolfgang Bangerth [Fri, 20 Aug 2021 01:37:32 +0000 (19:37 -0600)]
Use ParticleAccessor::set_properties() taking a Tensor.

3 years agoIntroduce ParticleAccessor::set_properties() taking a Tensor.
Wolfgang Bangerth [Fri, 20 Aug 2021 01:36:59 +0000 (19:36 -0600)]
Introduce ParticleAccessor::set_properties() taking a Tensor.

3 years agoUnify spelling. 12683/head
Wolfgang Bangerth [Fri, 20 Aug 2021 00:33:20 +0000 (18:33 -0600)]
Unify spelling.

3 years agoGroup member functions in hp::FECollection for doxygen. 12682/head
Wolfgang Bangerth [Thu, 19 Aug 2021 23:59:49 +0000 (17:59 -0600)]
Group member functions in hp::FECollection for doxygen.

3 years agoAdd changelog entry. 12681/head
Wolfgang Bangerth [Thu, 19 Aug 2021 23:19:42 +0000 (17:19 -0600)]
Add changelog entry.

3 years agoAdd a test.
Wolfgang Bangerth [Thu, 19 Aug 2021 23:19:34 +0000 (17:19 -0600)]
Add a test.

3 years agoAllow for the default construction of 1d face iterators.
Wolfgang Bangerth [Thu, 19 Aug 2021 23:18:46 +0000 (17:18 -0600)]
Allow for the default construction of 1d face iterators.

3 years agoPrevent excess timer output. 12680/head
Marc Fehling [Thu, 19 Aug 2021 22:29:10 +0000 (16:29 -0600)]
Prevent excess timer output.

3 years agoMerge pull request #12546 from TLJB/gdb-pprinter-fix
David Wells [Thu, 19 Aug 2021 20:58:32 +0000 (16:58 -0400)]
Merge pull request #12546 from TLJB/gdb-pprinter-fix

fix: Updated dotgdbinit to work with version 9.3

3 years agofix: Updated dotgdbinit to work with version 9.3 12546/head
Till Budde [Mon, 5 Jul 2021 11:48:12 +0000 (13:48 +0200)]
fix: Updated dotgdbinit to work with version 9.3

Fixed VectorPrinters and AlignedVectorPrinters dependency on
deprecated members data_begin and data_end.

3 years agoMerge pull request #12609 from bangerth/move-replicate
Timo Heister [Thu, 19 Aug 2021 17:06:50 +0000 (13:06 -0400)]
Merge pull request #12609 from bangerth/move-replicate

Fix AlignedVector memory management after replicate_across_communicator()

3 years agoMerge pull request #12677 from kronbichler/minor_cleanup
Daniel Arndt [Thu, 19 Aug 2021 14:28:39 +0000 (10:28 -0400)]
Merge pull request #12677 from kronbichler/minor_cleanup

tria.cc: Make exception and comment easier to read

3 years agoMerge pull request #12678 from amartinhuertas/patch-1
Peter Munch [Thu, 19 Aug 2021 12:12:06 +0000 (14:12 +0200)]
Merge pull request #12678 from amartinhuertas/patch-1

Typo fix in step-38 intro.dox

3 years agoTypo fix in step-38 intro.dox 12678/head
Alberto F. Martin [Thu, 19 Aug 2021 10:36:04 +0000 (20:36 +1000)]
Typo fix in step-38 intro.dox

3 years agotria.cc: Make exception and comment easier to read 12677/head
Martin Kronbichler [Thu, 19 Aug 2021 09:48:28 +0000 (11:48 +0200)]
tria.cc: Make exception and comment easier to read

3 years agoMerge pull request #12674 from bangerth/const
Marc Fehling [Thu, 19 Aug 2021 02:38:47 +0000 (20:38 -0600)]
Merge pull request #12674 from bangerth/const

Make TriaAccessor::set_*_boundary() const.

3 years agoMerge pull request #12675 from marcfehling/tria-faces
Wolfgang Bangerth [Thu, 19 Aug 2021 02:34:15 +0000 (20:34 -0600)]
Merge pull request #12675 from marcfehling/tria-faces

Avoid serializing lines on TriaFaces twice.

3 years agoMerge pull request #12670 from drwells/avoid-quadrature-copy
Wolfgang Bangerth [Thu, 19 Aug 2021 02:33:46 +0000 (20:33 -0600)]
Merge pull request #12670 from drwells/avoid-quadrature-copy

Avoid copying some Quadrature rules in FESystem.

3 years agoAvoid serializing lines on TriaFaces twice. 12675/head
Marc Fehling [Thu, 19 Aug 2021 00:33:37 +0000 (18:33 -0600)]
Avoid serializing lines on TriaFaces twice.

3 years agoMerge pull request #12671 from drwells/mark-as-protected
Wolfgang Bangerth [Thu, 19 Aug 2021 00:14:31 +0000 (18:14 -0600)]
Merge pull request #12671 from drwells/mark-as-protected

Mark some things as protected.

3 years agoMerge pull request #12672 from drwells/use-get-vertices
Wolfgang Bangerth [Thu, 19 Aug 2021 00:14:10 +0000 (18:14 -0600)]
Merge pull request #12672 from drwells/use-get-vertices

Use Mapping::get_vertices().

3 years agoMerge pull request #12673 from drwells/remove-unused-field
Wolfgang Bangerth [Thu, 19 Aug 2021 00:13:37 +0000 (18:13 -0600)]
Merge pull request #12673 from drwells/remove-unused-field

Remove an unused member variable.

3 years agoMerge pull request #12601 from marcfehling/fix-load-3
Wolfgang Bangerth [Thu, 19 Aug 2021 00:02:40 +0000 (18:02 -0600)]
Merge pull request #12601 from marcfehling/fix-load-3

Deprecated autopartition parameter in parallel::DistributedTriangulationBase::load()

3 years agoAdd test. 12674/head
Wolfgang Bangerth [Wed, 18 Aug 2021 22:42:24 +0000 (16:42 -0600)]
Add test.

3 years agoAdd a changelog entry.
Wolfgang Bangerth [Wed, 18 Aug 2021 22:37:29 +0000 (16:37 -0600)]
Add a changelog entry.

3 years agoMake TriaAccessor::set_*_boundary() const.
Wolfgang Bangerth [Wed, 18 Aug 2021 22:20:53 +0000 (16:20 -0600)]
Make TriaAccessor::set_*_boundary() const.

3 years agoMark some things as protected. 12671/head
David Wells [Wed, 18 Aug 2021 20:08:50 +0000 (16:08 -0400)]
Mark some things as protected.

These are protected, not private, in the base class.

3 years agoUse Mapping::get_vertices(). 12672/head
David Wells [Wed, 18 Aug 2021 20:07:59 +0000 (16:07 -0400)]
Use Mapping::get_vertices().

This is much faster for MappingFEField and about the same for other mapping
classes.

3 years agoRemove an unused member variable. 12673/head
David Wells [Wed, 18 Aug 2021 20:04:19 +0000 (16:04 -0400)]
Remove an unused member variable.

Now that we have hp::Collection (see 84f55c7706) this field is no longer used.

3 years agoAvoid copying some Quadrature rules in FESystem. 12670/head
David Wells [Wed, 18 Aug 2021 19:47:36 +0000 (15:47 -0400)]
Avoid copying some Quadrature rules in FESystem.

3 years agoDeprecated autopartition parameter. 12601/head
Marc Fehling [Sun, 25 Jul 2021 23:16:11 +0000 (17:16 -0600)]
Deprecated autopartition parameter.

3 years agoMerge pull request #12666 from drwells/fix-compilation-gcc-5.4
Wolfgang Bangerth [Wed, 18 Aug 2021 16:55:10 +0000 (10:55 -0600)]
Merge pull request #12666 from drwells/fix-compilation-gcc-5.4

Fix compilation with GCC 5.4.

3 years agoMerge pull request #12667 from zjiaqi2018/fix-matrixfree-typo
Daniel Arndt [Wed, 18 Aug 2021 15:17:57 +0000 (11:17 -0400)]
Merge pull request #12667 from zjiaqi2018/fix-matrixfree-typo

Fix documentation

3 years agofix documentation 12667/head
Jiaqi Zhang [Wed, 18 Aug 2021 14:29:02 +0000 (14:29 +0000)]
fix documentation

3 years agoFix compilation with GCC 5.4. 12666/head
David Wells [Wed, 18 Aug 2021 13:58:06 +0000 (09:58 -0400)]
Fix compilation with GCC 5.4.

Using constexpr here causes an internal compiler error when optimizations are
turned on.

3 years agoMerge pull request #12514 from simonsticko/mesh_classifier
Martin Kronbichler [Wed, 18 Aug 2021 08:41:10 +0000 (10:41 +0200)]
Merge pull request #12514 from simonsticko/mesh_classifier

Add class determining how cells/faces relate to a level set function

3 years agoAdd class determining how cells/faces relate to a level set function 12514/head
Simon Sticko [Mon, 28 Jun 2021 09:43:29 +0000 (11:43 +0200)]
Add class determining how cells/faces relate to a level set function

Add a class NonMatching::MeshClassifier that determines how the sign of
a level set function vary over the active cells/faces of a
triangulation. Each cell/face are associated with one of values of enum
LocationToLevelSet. If the level set is positive/negative over the
cell/face it is associated with LocationToLevelSet::inside/outside. If
the zero contour intersects the cell/face it is associated with
LocationToLevelSet::intersected.

3 years agoMerge pull request #12661 from bangerth/doc-3
David Wells [Tue, 17 Aug 2021 22:55:43 +0000 (18:55 -0400)]
Merge pull request #12661 from bangerth/doc-3

Update documentation of FE_Q constructor.

3 years agoUpdate documentation of FE_Q constructor. 12661/head
Wolfgang Bangerth [Mon, 16 Aug 2021 23:53:09 +0000 (17:53 -0600)]
Update documentation of FE_Q constructor.

3 years agoMerge pull request #12016 from konsim83/improve_conformity_tests
Martin Kronbichler [Tue, 17 Aug 2021 12:21:36 +0000 (14:21 +0200)]
Merge pull request #12016 from konsim83/improve_conformity_tests

Improve conformity tests of 2D vector elements

3 years agoMerge pull request #12662 from bangerth/doc-4
Martin Kronbichler [Tue, 17 Aug 2021 10:22:47 +0000 (12:22 +0200)]
Merge pull request #12662 from bangerth/doc-4

Minor doc update.

3 years agoMerge pull request #12612 from drwells/fix-mass-diagonal
Martin Kronbichler [Tue, 17 Aug 2021 10:01:44 +0000 (12:01 +0200)]
Merge pull request #12612 from drwells/fix-mass-diagonal

Make MassOperator's diagonals true diagonals, not lumped ones.

3 years agoMerge pull request #12596 from marcfehling/fix-load
Wolfgang Bangerth [Tue, 17 Aug 2021 03:05:15 +0000 (21:05 -0600)]
Merge pull request #12596 from marcfehling/fix-load

Load: align dealii and p4est meshes before repartitioning.

3 years agoMinor doc update. 12662/head
Wolfgang Bangerth [Tue, 17 Aug 2021 00:14:00 +0000 (18:14 -0600)]
Minor doc update.

3 years agoExplicitly use mass lumping in VectorTools::project_matrix_free(). 12612/head
David Wells [Sat, 14 Aug 2021 21:39:02 +0000 (17:39 -0400)]
Explicitly use mass lumping in VectorTools::project_matrix_free().

This restores the behavior present before the last commit and its generally the
right choice.

3 years agoImplement explicit mass lumping in MassOperator.
David Wells [Sat, 31 Jul 2021 16:14:17 +0000 (12:14 -0400)]
Implement explicit mass lumping in MassOperator.

3 years agoModify a test to check the diagonal values (not lumped).
David Wells [Sat, 31 Jul 2021 14:31:56 +0000 (10:31 -0400)]
Modify a test to check the diagonal values (not lumped).

3 years agoMake MassOperator's diagonals true diagonals, not lumped ones.
David Wells [Thu, 29 Jul 2021 21:40:18 +0000 (17:40 -0400)]
Make MassOperator's diagonals true diagonals, not lumped ones.

The previous implementation used row sums (i.e., mass lumping) instead of
computing the diagonal. This should be changed for two reasons:
1. It's not really the diagonal (and its inconsistent with LaplaceOperator).
2. Lumped mass matrices are only positive for either low-order elements (where
   the basis functions are all positive) or very well-behaved elements (like
   FE_Q). In particular, the values on the diagonal are either zero or negative
   for TET10 if we use lumping, which isn't going to work.


In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.