Merge pull request #26462 from cticenhour/petsc-3203

Update PETSc to 3.20.3, update libmesh

Author: loganharbour

conda/libmesh/conda_build_config.yaml

@@ -1,5 +1,5 @@
 moose_petsc:
-  - moose-petsc 3.20.1 build_0
+  - moose-petsc 3.20.3 build_0
 
 moose_libmesh_vtk:
   - moose-libmesh-vtk 9.2.6 build_8

conda/libmesh/meta.yaml

@@ -5,7 +5,7 @@
 # As well as any directions pertaining to modifying those files.
 # ALSO: Follow the directions in scripts/tests/versioner_hashes.yaml
 {% set build = 0 %}
-{% set version = "2023.11.29" %}
+{% set version = "2024.01.23" %}
 {% set vtk_friendly_version = "9.2" %}
 
 # permanent

conda/moose-dev/conda_build_config.yaml

@@ -1,5 +1,5 @@
 moose_libmesh:
-  - moose-libmesh 2023.11.29 build_0
+  - moose-libmesh 2024.01.23 build_0
 
 moose_wasp:
   - moose-wasp 2023.11.29

conda/moose-dev/meta.yaml

@@ -6,7 +6,7 @@
 #   framework/doc/packages_config.yml
 #
 # As well as any directions pertaining to modifying those files.
-{% set version = "2023.12.20" %}
+{% set version = "2024.01.23" %}
 
 package:
   name: moose-dev

conda/moose/conda_build_config.yaml

@@ -1,5 +1,5 @@
 moose_dev:
-  - moose-dev 2023.12.20
+  - moose-dev 2024.01.23
 
 #### Darwin SDK SYSROOT
 CONDA_BUILD_SYSROOT:                                        # [osx]

conda/petsc/meta.yaml

@@ -8,7 +8,7 @@
 # As well as any directions pertaining to modifying those files.
 # ALSO: Follow the directions in scripts/tests/versioner_hashes.yaml
 {% set build = 0 %}
-{% set version = "3.20.1" %}
+{% set version = "3.20.3" %}
 
 # permanent
 {% set strbuild = "build_" + build|string %}

framework/doc/packages_config.yml

@@ -5,8 +5,8 @@ minimum_clang: 10.0.1
 minimum_python: 3.7
 mpich: 4.0.2
 petsc_alt: 3.11.4
-petsc: 3.20.1
-moose_dev: 2023.12.20
+petsc: 3.20.3
+moose_dev: 2024.01.23
 moose_tools: 2023.12.20
 apptainer_rockylinux: 8.8.20230518
 apptainer_mpich: 3.4.3

libmesh

@@ -9,6 +9,44 @@ for a complete description of all MOOSE changes.
 
 ## libMesh-level Changes
 
+### `2024.01.23` Update
+
+- New example illustrating a Hybridized Discontinuous Galerkin solve
+- Added new finite element support: the `HIERARCHIC_VEC` FE family,
+  for vector-valued hierarchic shape functions, and the
+  `L2_HIERARCHIC_VEC` FE family, for the broken/L2 version of that
+  space.  Added simple-vector-valued example and
+  Hybridized-Discontinuous-Galerkin example runs using these families.
+- Support for a pre-check object callback from PetscNonlinearSolver
+- Keep track of libMesh command line queries, and avoid spurious
+  unused-option warnings about them from PETSc
+- Support for upgrading the order of mixed-order input meshes in
+  `all_second_order()` and `all_complete_order()`
+- Support for neighbor-finding and mesh stitching with mixed-order
+  input meshes.  Mixed-order meshes must still be upgraded before
+  being used to support C0 finite element types in a linear solve
+  (where constraint equations would be needed between
+  incompatible-order neighbors), but mesh generation and modification
+  with pre-upgraded mixed-order meshes is now possible.
+- Support for more preconditioning options in EigenSparseLinearSolver;
+  this allows us to reenable a dozen examples that were previously
+  skipped when Eigen was configured as the default linear algebra
+  package.
+- Minor fixes to reduced basis code, string parsing of uncommon
+  quadrature rule names, configuration on systems without stdlib.h,
+- Minor clarifications and typo fixes in documentation, comments, and
+  runtime error messages
+- Fixes for VTK .pvtu reading, and preservation of libMesh node, element,
+  and subdomain ids written to VTK files
+- `Elem::volume()` no longer can throw an exception when called on an
+  element without an invertible mapping; hence `print_info()` and
+  `get_info()` are now safe to use when emitting debugging messages about
+  twisted geometric elements.
+- Support `interior_parent` in `Elem` order upgrades
+- Avoid creating duplicate entries in `BoundaryInfo` multimaps
+- Removes an overzealous assertion in variable name collection that
+  does not work with `NodeElem`-only meshes
+
 ## PETSc-level Changes
 
 ## Bug Fixes and Minor Enhancements

modules/doc/content/newsletter/2024/2024_01.md

@@ -2,6 +2,7 @@
 <VTKFile type="PUnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
   <PUnstructuredGrid GhostLevel="1">
     <PPointData RationalWeights="RationalWeights">
+      <PDataArray type="Int32" Name="libmesh_node_id"/>
       <PDataArray type="Float64" Name="RationalWeights"/>
       <PDataArray type="Float64" Name="von_mises"/>
       <PDataArray type="Float64" Name="Max_Princ"/>