]> https://gitweb.dealii.org/ - dealii.git/log
dealii.git
5 months agoMerge pull request #15773 from luca-heltai/fe-coupling-values
Wolfgang Bangerth [Thu, 11 Apr 2024 17:37:05 +0000 (11:37 -0600)]
Merge pull request #15773 from luca-heltai/fe-coupling-values

FECouplingValues.

5 months agoFECouplingValues. 15773/head
Luca Heltai [Thu, 20 Jul 2023 22:54:48 +0000 (00:54 +0200)]
FECouplingValues.

5 months agoMerge pull request #16876 from drwells/hdf5-floating-point-exceptions
Wolfgang Bangerth [Wed, 10 Apr 2024 20:34:40 +0000 (14:34 -0600)]
Merge pull request #16876 from drwells/hdf5-floating-point-exceptions

Explicitly open HDF5 and then clear its floating-point exceptions.

5 months agoExplicitly open HDF5 and then clear its floating-point exceptions. 16876/head
David Wells [Tue, 9 Apr 2024 23:13:57 +0000 (19:13 -0400)]
Explicitly open HDF5 and then clear its floating-point exceptions.

New versions of HDF5 detect the floating-point environment by doing several
things which trigger floating point exceptions, i.e., programs in the test suite
now abort in H5open(). Work around that by explicitly setting up HDF5 before we
turn on the 'abort on floating point exceptions' code.

5 months agoMerge pull request #16869 from peterrum/extract_level_constant_modes
Bruno Blais [Tue, 9 Apr 2024 15:11:47 +0000 (17:11 +0200)]
Merge pull request #16869 from peterrum/extract_level_constant_modes

Add DoFTools::extract_level_constant_modes()

5 months agoMerge pull request #15572 from Rombur/init_finalize
Daniel Arndt [Tue, 9 Apr 2024 12:46:39 +0000 (08:46 -0400)]
Merge pull request #15572 from Rombur/init_finalize

Add InitFinalize base class of MPI_InitFinalize

5 months agoMerge pull request #16875 from bangerth/typos
Martin Kronbichler [Tue, 9 Apr 2024 07:41:02 +0000 (09:41 +0200)]
Merge pull request #16875 from bangerth/typos

Fix a whole lot of typos.

5 months agoMerge pull request #16870 from bangerth/distort-random
David Wells [Tue, 9 Apr 2024 01:41:28 +0000 (21:41 -0400)]
Merge pull request #16870 from bangerth/distort-random

Update documentation for GridTools::distort_random().

5 months agoUpdate documentation for GridTools::distort_random(). 16870/head
Wolfgang Bangerth [Sun, 7 Apr 2024 22:31:53 +0000 (16:31 -0600)]
Update documentation for GridTools::distort_random().

In particular also remove a reference to the function in the documentation of class
Triangulation. The function used to be part of that class, but has not been for a long
time.

5 months agoFix a whole lot of typos. 16875/head
Wolfgang Bangerth [Mon, 8 Apr 2024 20:32:28 +0000 (14:32 -0600)]
Fix a whole lot of typos.

5 months agoMerge pull request #16874 from bangerth/exc
David Wells [Mon, 8 Apr 2024 22:25:31 +0000 (18:25 -0400)]
Merge pull request #16874 from bangerth/exc

Remove two exceptions that are apparently no longer used.

5 months agoRemove two exceptions that are apparently no longer used. 16874/head
Wolfgang Bangerth [Mon, 8 Apr 2024 20:27:49 +0000 (14:27 -0600)]
Remove two exceptions that are apparently no longer used.

5 months agoMerge pull request #16871 from tjhei/maybe-uninit-orthogonal
Wolfgang Bangerth [Mon, 8 Apr 2024 17:29:48 +0000 (11:29 -0600)]
Merge pull request #16871 from tjhei/maybe-uninit-orthogonal

fix maybe-unused warning in grid_tools_dofhandler

5 months agoMerge pull request #16868 from masterleinad/remove_redundant_copy_assignment_vectoriz...
Wolfgang Bangerth [Mon, 8 Apr 2024 17:25:17 +0000 (11:25 -0600)]
Merge pull request #16868 from masterleinad/remove_redundant_copy_assignment_vectorizedarrayiterator

Remove redundant defaulted copy assignment operator for VectorizedArrayIterator

5 months agofix maybe-unused warning 16871/head
Timo Heister [Sun, 7 Apr 2024 23:37:12 +0000 (19:37 -0400)]
fix maybe-unused warning

```
./include/deal.II/grid/reference_cell.h:2799:22: warning:
‘face1_vertices.std::array<unsigned[763/12022] M_elems[0]’ may be used
uninitialized in this function [-Wmaybe-uninitialized]
./source/grid/grid_tools_dof_handlers.cc:2438:7: note:
‘face1_vertices.std::array<unsigned int, 1>::_M_elems[0]’ was declared
here
```

5 months agoAdd DoFTools::extract_level_constant_modes() 16869/head
Peter Munch [Wed, 3 Apr 2024 11:49:21 +0000 (13:49 +0200)]
Add DoFTools::extract_level_constant_modes()

5 months agoMerge pull request #15819 from luca-heltai/reordered-view
Wolfgang Bangerth [Sun, 7 Apr 2024 19:04:49 +0000 (13:04 -0600)]
Merge pull request #15819 from luca-heltai/reordered-view

FEValuesViews::RenumberedView

5 months agoRemove redundant defaulted copy assignment operator for VectorizedArrayIterator 16868/head
Daniel Arndt [Sun, 7 Apr 2024 16:16:08 +0000 (12:16 -0400)]
Remove redundant defaulted copy assignment operator for VectorizedArrayIterator

5 months agoMerge pull request #16866 from peterrum/step-89-fixes
Martin Kronbichler [Sun, 7 Apr 2024 09:01:35 +0000 (11:01 +0200)]
Merge pull request #16866 from peterrum/step-89-fixes

Small fixes in step-89

5 months agoMerge pull request #16864 from bangerth/simplify
Matthias Maier [Sun, 7 Apr 2024 03:25:00 +0000 (03:25 +0000)]
Merge pull request #16864 from bangerth/simplify

Simplify some functions by just deferring to member functions.

5 months agoMerge pull request #16863 from bangerth/40-2
Matthias Maier [Sun, 7 Apr 2024 03:21:27 +0000 (03:21 +0000)]
Merge pull request #16863 from bangerth/40-2

Unclutter step-40's run() function.

5 months agoMerge pull request #16862 from bangerth/40-1
Matthias Maier [Sun, 7 Apr 2024 03:21:09 +0000 (03:21 +0000)]
Merge pull request #16862 from bangerth/40-1

In step-40, move timer section into the function being timed.

5 months agoMerge pull request #16835 from drwells/fix-trailing-newline-mac
Matthias Maier [Sun, 7 Apr 2024 03:20:27 +0000 (03:20 +0000)]
Merge pull request #16835 from drwells/fix-trailing-newline-mac

Fix ensure_single_trailing_newline() on macs.

5 months agoSmall fixes in step-89 16866/head
Peter Munch [Sat, 6 Apr 2024 20:43:50 +0000 (22:43 +0200)]
Small fixes in step-89

5 months agoMerge pull request #16299 from jh66637/step-89-tutorial
Bruno Blais [Sat, 6 Apr 2024 17:36:33 +0000 (10:36 -0700)]
Merge pull request #16299 from jh66637/step-89-tutorial

Add tutorial on Nitsche-type mortaring

5 months agoIn step-40, move timer section into the function being timed. 16862/head
Wolfgang Bangerth [Sat, 6 Apr 2024 16:19:39 +0000 (10:19 -0600)]
In step-40, move timer section into the function being timed.

5 months agoSimplify some functions by just deferring to member functions. 16864/head
Wolfgang Bangerth [Sat, 6 Apr 2024 17:02:47 +0000 (11:02 -0600)]
Simplify some functions by just deferring to member functions.

5 months agoUnclutter step-40's run() function. 16863/head
Wolfgang Bangerth [Sat, 6 Apr 2024 16:22:37 +0000 (10:22 -0600)]
Unclutter step-40's run() function.

Move detailed output into the function that actually creates this information.
This ensures that the run() function itself remains rather clean and easy to
follow.

5 months agoMerge pull request #16858 from tamiko/fix_cgal_test_fpe
Matthias Maier [Sat, 6 Apr 2024 05:35:00 +0000 (05:35 +0000)]
Merge pull request #16858 from tamiko/fix_cgal_test_fpe

Test: cgal/cgal_mesh_criteria: disable floating point exceptions

5 months agoMerge pull request #16840 from tamiko/vectorization_type_trait
Wolfgang Bangerth [Sat, 6 Apr 2024 04:36:23 +0000 (22:36 -0600)]
Merge pull request #16840 from tamiko/vectorization_type_trait

Vectorization type trait

5 months agoMerge pull request #16859 from peterrum/remove_debug_output
Matthias Maier [Sat, 6 Apr 2024 03:34:16 +0000 (03:34 +0000)]
Merge pull request #16859 from peterrum/remove_debug_output

Remove debug output

5 months agoTest: cgal/cgal_mesh_criteria: disable floating point exceptions 16858/head
Matthias Maier [Fri, 5 Apr 2024 19:57:13 +0000 (14:57 -0500)]
Test: cgal/cgal_mesh_criteria: disable floating point exceptions

We trigger a (deliberate) floating point exception in the CGAL (header)
library:
```
    // The multiplications could produce some NaN, with 0 * inf. Replacing it with inf is safe.
    // min(x,y) (the order is essential) returns its second argument when the first is NaN.
    // An IEEE 754-2019 maximum could help.
    __m128d big = IA::largest().simd();
->  __m128d x1 = _mm_mul_pd(aa,bz);                       // {-ai*bi,as*bs}
    //x1 = _mm_min_pd(x1,big); // no NaN
    __m128d x2 = _mm_mul_pd(aa,c);                        // {-ai*bs,as*bi}
    x2 = _mm_min_pd(x2,big); // no NaN
    __m128d x3 = _mm_mul_pd(ap,bz);                       // {-as*bi,ai*bs}
    //x3 = _mm_min_pd(x3,big); // no NaN
    __m128d x4 = _mm_mul_pd(ap,c);                        // {-as*bs,ai*bi}
    x4 = _mm_min_pd(x4,big); // no NaN
```

Thus, disable floating point exceptions for the test.

5 months agoRemove debug output 16859/head
Peter Munch [Fri, 5 Apr 2024 20:22:55 +0000 (22:22 +0200)]
Remove debug output

5 months agoMerge pull request #16853 from kronbichler/avoid_allocations_2
Wolfgang Bangerth [Fri, 5 Apr 2024 17:37:00 +0000 (11:37 -0600)]
Merge pull request #16853 from kronbichler/avoid_allocations_2

Add Quadrature::initialize() function taking ArrayView

5 months agoMerge pull request #16624 from peterrum/matrixcreator_fullmatrix
Wolfgang Bangerth [Fri, 5 Apr 2024 17:29:29 +0000 (11:29 -0600)]
Merge pull request #16624 from peterrum/matrixcreator_fullmatrix

Enable MatrixCreator functions also for FullMatrix

5 months agoMerge both initialize functions 16853/head
Martin Kronbichler [Fri, 5 Apr 2024 14:54:37 +0000 (16:54 +0200)]
Merge both initialize functions

5 months agoMerge pull request #16852 from kronbichler/fix_docu
Wolfgang Bangerth [Fri, 5 Apr 2024 14:27:11 +0000 (08:27 -0600)]
Merge pull request #16852 from kronbichler/fix_docu

MappingQ: Remove outdated documentation

5 months agoMerge pull request #16855 from kronbichler/fix_pow
Daniel Arndt [Fri, 5 Apr 2024 13:31:49 +0000 (09:31 -0400)]
Merge pull request #16855 from kronbichler/fix_pow

Avoid std::pow in BarycentricPolynomials

5 months agoAvoid std::pow in BarycentricPolynomials 16855/head
Martin Kronbichler [Fri, 5 Apr 2024 10:31:30 +0000 (12:31 +0200)]
Avoid std::pow in BarycentricPolynomials

5 months agoNew test case
Martin Kronbichler [Fri, 5 Apr 2024 08:47:09 +0000 (10:47 +0200)]
New test case

5 months agoAdd Quadrature::initialize() function taking ArrayView
Martin Kronbichler [Fri, 5 Apr 2024 08:30:17 +0000 (10:30 +0200)]
Add Quadrature::initialize() function taking ArrayView

5 months agoRemove outdated documentation 16852/head
Martin Kronbichler [Fri, 5 Apr 2024 08:02:14 +0000 (10:02 +0200)]
Remove outdated documentation

5 months agoMerge pull request #16849 from tamiko/fix_fe_values_tests
Martin Kronbichler [Fri, 5 Apr 2024 07:05:18 +0000 (09:05 +0200)]
Merge pull request #16849 from tamiko/fix_fe_values_tests

Test fe/fe_values_view_{07|10|13|16}: make Assert more robust

5 months agoTest fe/fe_values_view_{07|10|13|16}: make Assert more robust 16849/head
Matthias Maier [Fri, 5 Apr 2024 03:19:25 +0000 (22:19 -0500)]
Test fe/fe_values_view_{07|10|13|16}: make Assert more robust

The Assert checks that we computed with a relative error of less than
1e-12. This fails however, if the vector in question is 0. Thus, also
check for an absolute error.

5 months agoEnable MatrixCreator functions also for FullMatrix 16624/head
Peter Munch [Sat, 10 Feb 2024 11:48:15 +0000 (12:48 +0100)]
Enable MatrixCreator functions also for FullMatrix

5 months agoMerge pull request #16847 from kronbichler/avoid_allocations
Martin Kronbichler [Thu, 4 Apr 2024 18:13:53 +0000 (20:13 +0200)]
Merge pull request #16847 from kronbichler/avoid_allocations

MappingQ: Favor get_vertices() over compute_mapping_support_points() if possible

5 months agoMerge pull request #16830 from kronbichler/mapping_cartesian
Daniel Arndt [Thu, 4 Apr 2024 16:39:19 +0000 (12:39 -0400)]
Merge pull request #16830 from kronbichler/mapping_cartesian

MappingCartesian: Some performance improvements

5 months agoMappingQ: Use get_vertices() instead of compute_mapping_support_points when possible 16847/head
Martin Kronbichler [Wed, 3 Apr 2024 20:46:41 +0000 (22:46 +0200)]
MappingQ: Use get_vertices() instead of compute_mapping_support_points when possible

5 months agoMerge pull request #16800 from kronbichler/fe_eval_docu
Wolfgang Bangerth [Thu, 4 Apr 2024 05:48:35 +0000 (23:48 -0600)]
Merge pull request #16800 from kronbichler/fe_eval_docu

FEEvaluation: Modernize documentation of get_xxx/submit_xxx

5 months agoMerge pull request #16843 from kronbichler/improve_vector_tools_interpolate
Matthias Maier [Thu, 4 Apr 2024 04:38:28 +0000 (04:38 +0000)]
Merge pull request #16843 from kronbichler/improve_vector_tools_interpolate

VectorTools::interpolate: Make common case of vector problem faster

5 months agoMerge pull request #16839 from bangerth/tensor-product
Matthias Maier [Thu, 4 Apr 2024 04:34:27 +0000 (04:34 +0000)]
Merge pull request #16839 from bangerth/tensor-product

Make tensor*tensor implementation easier to read.

5 months agoVectorizedArray: add constexpr boolean that records whether a type is vectorized 16840/head
Matthias Maier [Tue, 2 Apr 2024 20:55:11 +0000 (15:55 -0500)]
VectorizedArray: add constexpr boolean that records whether a type is vectorized

5 months agoVectorization: move assert into constructor
Matthias Maier [Tue, 2 Apr 2024 21:07:17 +0000 (16:07 -0500)]
Vectorization: move assert into constructor

5 months agoMerge pull request #16844 from kronbichler/eval_interface
Martin Kronbichler [Wed, 3 Apr 2024 18:03:54 +0000 (20:03 +0200)]
Merge pull request #16844 from kronbichler/eval_interface

Use ArrayView for interface to tensor-product point evaluator

5 months agoMerge pull request #16845 from kronbichler/pass_by_reference
Wolfgang Bangerth [Wed, 3 Apr 2024 17:11:43 +0000 (11:11 -0600)]
Merge pull request #16845 from kronbichler/pass_by_reference

Pass some SmartPointer variables by reference

5 months agoFEValuesViews::ReorderedView 15819/head
Luca Heltai [Wed, 2 Aug 2023 00:22:34 +0000 (02:22 +0200)]
FEValuesViews::ReorderedView

5 months agoFix ensure_single_trailing_newline() on macs. 16835/head
David Wells [Tue, 2 Apr 2024 02:16:18 +0000 (22:16 -0400)]
Fix ensure_single_trailing_newline() on macs.

I'm not sure exactly why mac's version of sed produces different output here but
it is easy enough to use something simpler. Just check if the last character is
a newline: if not, then append a newline character.

5 months agoUse ArrayView for interface to tensor-product point evaluator 16844/head
Martin Kronbichler [Wed, 3 Apr 2024 11:39:35 +0000 (13:39 +0200)]
Use ArrayView for interface to tensor-product point evaluator

5 months agoMerge pull request #16832 from bangerth/fix-tests
Martin Kronbichler [Wed, 3 Apr 2024 13:10:51 +0000 (15:10 +0200)]
Merge pull request #16832 from bangerth/fix-tests

Fix naming of output files.

5 months agoMerge pull request #16841 from bangerth/type-to-string
Martin Kronbichler [Wed, 3 Apr 2024 13:10:32 +0000 (15:10 +0200)]
Merge pull request #16841 from bangerth/type-to-string

Augment documentation of Utilities::type_to_string().

5 months agoChangelog 16843/head
Martin Kronbichler [Wed, 3 Apr 2024 10:52:07 +0000 (12:52 +0200)]
Changelog

5 months agoVectorTools::interpolate: Make common case of vector problem faster
Martin Kronbichler [Wed, 3 Apr 2024 10:50:35 +0000 (12:50 +0200)]
VectorTools::interpolate: Make common case of vector problem faster

5 months agoPass some SmartPointer variables by reference 16845/head
Martin Kronbichler [Wed, 3 Apr 2024 12:01:13 +0000 (14:01 +0200)]
Pass some SmartPointer variables by reference

5 months agoMerge pull request #16837 from bangerth/vec-array-2
Martin Kronbichler [Wed, 3 Apr 2024 07:28:19 +0000 (09:28 +0200)]
Merge pull request #16837 from bangerth/vec-array-2

Provide an implementation for VectorizedArray::dot_product().

5 months agoAugment documentation of Utilities::type_to_string(). 16841/head
Wolfgang Bangerth [Wed, 3 Apr 2024 03:09:14 +0000 (21:09 -0600)]
Augment documentation of Utilities::type_to_string().

5 months agoMerge pull request #16833 from kronbichler/fgmres_orth
Rene Gassmoeller [Tue, 2 Apr 2024 18:20:10 +0000 (14:20 -0400)]
Merge pull request #16833 from kronbichler/fgmres_orth

Switch SolverFGMRES to new delayed classical Gram-Schmidt

5 months agoImprove two assertions 16830/head
Martin Kronbichler [Tue, 2 Apr 2024 16:37:29 +0000 (18:37 +0200)]
Improve two assertions

5 months agoMake tensor*tensor implementation easier to read. 16839/head
Wolfgang Bangerth [Tue, 2 Apr 2024 16:11:03 +0000 (10:11 -0600)]
Make tensor*tensor implementation easier to read.

5 months agoChangelog
Martin Kronbichler [Mon, 1 Apr 2024 19:22:43 +0000 (21:22 +0200)]
Changelog

5 months agoTest case
Martin Kronbichler [Mon, 1 Apr 2024 19:20:36 +0000 (21:20 +0200)]
Test case

5 months agoMappingCartesian: Some performance improvements
Martin Kronbichler [Mon, 1 Apr 2024 19:20:23 +0000 (21:20 +0200)]
MappingCartesian: Some performance improvements

5 months agoMerge pull request #16828 from drwells/qprojector-same-orientation
Martin Kronbichler [Tue, 2 Apr 2024 15:10:12 +0000 (17:10 +0200)]
Merge pull request #16828 from drwells/qprojector-same-orientation

QProjector: internally use the combined orientation.

5 months agoProvide a default implementation for VectorizedArray::dot_product(). 16837/head
Wolfgang Bangerth [Tue, 2 Apr 2024 03:23:08 +0000 (21:23 -0600)]
Provide a default implementation for VectorizedArray::dot_product().

5 months agoQProjector: internally use the combined orientation. 16828/head
David Wells [Sun, 31 Mar 2024 20:02:27 +0000 (16:02 -0400)]
QProjector: internally use the combined orientation.

This replaces a unique indexing scheme for orientations with the standard one.

5 months agoFE_BDM: un-hardcode some QProjector properties.
David Wells [Tue, 2 Apr 2024 13:24:31 +0000 (09:24 -0400)]
FE_BDM: un-hardcode some QProjector properties.

5 months agoMerge pull request #16836 from bangerth/vec-array
Martin Kronbichler [Tue, 2 Apr 2024 07:13:37 +0000 (09:13 +0200)]
Merge pull request #16836 from bangerth/vec-array

Minor updates to VectorizedArray.

5 months agoMinor updates to VectorizedArray. 16836/head
Wolfgang Bangerth [Tue, 2 Apr 2024 03:11:14 +0000 (21:11 -0600)]
Minor updates to VectorizedArray.

5 months agoMerge pull request #16831 from kronbichler/fix_n_vertices
David Wells [Tue, 2 Apr 2024 02:35:01 +0000 (22:35 -0400)]
Merge pull request #16831 from kronbichler/fix_n_vertices

Fix two places where TriaAccessor::n_vertices is called repeatedly

5 months agoSwitch SolverFGMRES to new delayed classical Gram-Schmidt 16833/head
Martin Kronbichler [Mon, 1 Apr 2024 21:04:37 +0000 (23:04 +0200)]
Switch SolverFGMRES to new delayed classical Gram-Schmidt

5 months agoFix naming of output files. 16832/head
Wolfgang Bangerth [Mon, 1 Apr 2024 17:55:40 +0000 (11:55 -0600)]
Fix naming of output files.

5 months agoFix two places where CellAccessor::n_vertices is called repeatedly 16831/head
Martin Kronbichler [Mon, 1 Apr 2024 19:27:41 +0000 (21:27 +0200)]
Fix two places where CellAccessor::n_vertices is called repeatedly

5 months agoMerge pull request #16821 from drwells/qprojector-description-orientation
Martin Kronbichler [Mon, 1 Apr 2024 16:49:47 +0000 (18:49 +0200)]
Merge pull request #16821 from drwells/qprojector-description-orientation

Deprecate some more (bool, bool, bool) orientation functions

5 months agoImprove a comment. 16821/head
David Wells [Mon, 1 Apr 2024 13:12:28 +0000 (09:12 -0400)]
Improve a comment.

5 months agoRemove calls to said functions on the library.
David Wells [Sat, 30 Mar 2024 12:24:29 +0000 (08:24 -0400)]
Remove calls to said functions on the library.

5 months agoDeprecate (bool, bool, bool) QProjector functions.
David Wells [Sat, 30 Mar 2024 11:57:44 +0000 (07:57 -0400)]
Deprecate (bool, bool, bool) QProjector functions.

5 months agoMerge pull request #16804 from bangerth/comsol-tri
Daniel Arndt [Mon, 1 Apr 2024 12:33:36 +0000 (08:33 -0400)]
Merge pull request #16804 from bangerth/comsol-tri

Add a test for COMSOL .mphtxt triangle meshes.

5 months agoMerge pull request #16823 from bangerth/unalign
Martin Kronbichler [Sun, 31 Mar 2024 18:10:37 +0000 (20:10 +0200)]
Merge pull request #16823 from bangerth/unalign

Revert aligning tensors

5 months agoAdd a test for COMSOL .mphtxt triangle meshes. 16804/head
Wolfgang Bangerth [Thu, 28 Mar 2024 18:35:31 +0000 (12:35 -0600)]
Add a test for COMSOL .mphtxt triangle meshes.

5 months agoRevert "No longer output the size of tensors into output file of a test." 16823/head
Wolfgang Bangerth [Sat, 30 Mar 2024 15:04:10 +0000 (09:04 -0600)]
Revert "No longer output the size of tensors into output file of a test."

This reverts commit a492d57657769bdc38a65c2532f1711a55abab21.

5 months agoRevert "Add a changelog entry for the alignment of tensors."
Wolfgang Bangerth [Sat, 30 Mar 2024 15:03:43 +0000 (09:03 -0600)]
Revert "Add a changelog entry for the alignment of tensors."

This reverts commit d14a12db1cc7b91380a143350e613442362c876e.

5 months agoRevert "Align Tensor<rank,dim,Number> to make vectorization possible."
Wolfgang Bangerth [Sat, 30 Mar 2024 15:03:26 +0000 (09:03 -0600)]
Revert "Align Tensor<rank,dim,Number> to make vectorization possible."

This reverts commit d6a65a9b929b1862a78de843ceaa7848e5a588e3.

5 months agoMerge pull request #16820 from tamiko/fix_avx512_fpe
Wolfgang Bangerth [Sat, 30 Mar 2024 15:01:57 +0000 (09:01 -0600)]
Merge pull request #16820 from tamiko/fix_avx512_fpe

distributed/tria.cc: avoid FPE due to undefined behavior

5 months agoMerge pull request #16814 from simonsticko/codim_in_getname
Wolfgang Bangerth [Sat, 30 Mar 2024 14:59:56 +0000 (08:59 -0600)]
Merge pull request #16814 from simonsticko/codim_in_getname

Add codimension in get_name() for a few elements

5 months agodistributed/tria.cc: avoid FPE due to undefined behavior 16820/head
Matthias Maier [Sat, 30 Mar 2024 00:49:16 +0000 (19:49 -0500)]
distributed/tria.cc: avoid FPE due to undefined behavior

The following construct
```
static_cast<types::subdomain_id>(this_sc_point[dim])
```
triggers a floating point exception when the stored double is -1.0 (at
least with gcc-11 and avx512 instructions enabled).

Judging from cppreference [1] this is undefined behavior:
"""
A prvalue of floating-point type can be converted to a prvalue of any
integer type. The fractional part is truncated, that is, the fractional
part is discarded.

If the truncated value cannot fit into the destination type, the
behavior is undefined (even when the destination type is unsigned,
modulo arithmetic does not apply).
"""

Thus, work around the issue by explicitly setting negative values to
`numbers::invalid_subdomain_id` if we encounter -1.

[1] https://en.cppreference.com/w/cpp/language/implicit_conversion

5 months agoMerge pull request #16816 from bangerth/fix-test
Matthias Maier [Sat, 30 Mar 2024 00:28:21 +0000 (00:28 +0000)]
Merge pull request #16816 from bangerth/fix-test

No longer output the size of tensors into output file of a test.

5 months agoMerge pull request #16815 from bangerth/fix-particle-size-2
Matthias Maier [Sat, 30 Mar 2024 00:27:36 +0000 (00:27 +0000)]
Merge pull request #16815 from bangerth/fix-particle-size-2

Fix a second place where particles compute the size of a buffer.

5 months agoNo longer output the size of tensors into output file of a test. 16816/head
Wolfgang Bangerth [Fri, 29 Mar 2024 22:30:32 +0000 (16:30 -0600)]
No longer output the size of tensors into output file of a test.

5 months agoFix a second place where particles compute the size of a buffer. 16815/head
Wolfgang Bangerth [Fri, 29 Mar 2024 19:21:50 +0000 (13:21 -0600)]
Fix a second place where particles compute the size of a buffer.

5 months agoMerge pull request #16813 from bangerth/fix-particle-size
Martin Kronbichler [Fri, 29 Mar 2024 19:49:19 +0000 (20:49 +0100)]
Merge pull request #16813 from bangerth/fix-particle-size

Allocate the right amount of memory when serializing particles.

5 months agoMerge pull request #16788 from kronbichler/fix_asserts
Daniel Arndt [Fri, 29 Mar 2024 17:30:56 +0000 (13:30 -0400)]
Merge pull request #16788 from kronbichler/fix_asserts

FEEvaluation: Add static asserts to functions selected by enable_if


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.