Website changes for MFEM v4.9

src/about.md

@@ -1,10 +1,26 @@
 # About MFEM
 
-MFEM originates from previous research effort in the (unreleased) [AggieFEM/aFEM](https://www.math.tamu.edu/research/vigre/archive/2000c-Lazarov.html) project.
+MFEM originates from previous research effort in the (unreleased) [AggieFEM/aFEM](img/2000c-Lazarov.png) project.
+
+<!-- The AggieFEM/aFEM project used to be https://www.math.tamu.edu/research/vigre/archive/2000c-Lazarov.html. The above screenshot was taken in 2024. -->
 
 Please cite with:
 ```c
-@article{mfem,
+@article{mfem-2024,
+  title = {High-Performance Finite Elements with {MFEM}},
+  author = {J. Andrej and N. Atallah and J.-P. Bäcker and  J.-S. Camier and
+            D. Copeland and V. Dobrev and Y. Dudouit and T. Duswald and
+            B. Keith and D. Kim and T. Kolev and B. Lazarov and  K. Mittal and
+            W. Pazner and S. Petrides and S. Shiraiwa and M. Stowell and V. Tomov},
+  journal={The International Journal of High Performance Computing Applications},
+  volume={38},
+  number={5},
+  pages={447-467},
+  year={2024},
+  publisher={SAGE Publications Sage UK: London, England}
+}
+
+@article{mfem-2021,
   title = {{MFEM}: A Modular Finite Element Methods Library},
   author = {R. Anderson and J. Andrej and A. Barker and J. Bramwell and J.-S. Camier and
             J. Cerveny and V. Dobrev and Y. Dudouit and A. Fisher and Tz. Kolev and W. Pazner and
@@ -79,6 +95,7 @@ Please cite with:
 - [Dohyun Kim](https://github.com/dohyun-cse)
 - [Patrick Knupp](https://dblp.org/pid/38/3416.html)
 - [Tzanio Kolev](https://people.llnl.gov/kolev1) &mdash; Project Leader
+- [Anthony Kolshorn]()
 - [Chris Laganella](https://github.com/Laganellac)
 - [Justin Laughlin](https://github.com/justinlaughlin)
 - [Ilya Lashuk](https://scholar.google.com/citations?user=5KQiAUoAAAAJ&hl=en)

src/bilininteg.md

@@ -92,7 +92,7 @@ Note that any operators involving a derivative of the range function
 $v$ or $\vec\{v}$ are computed using integration by parts.  This leads
 to a boundary integral which can be used to apply Neumann boundary
 conditions.  Some of these operators are listed along with their
-boundary terms in section [Weak Operators](#weak-operators-and-their-boundary-integrals).
+boundary terms in section [Weak Operators and Their Boundary Integrals](#weak-operators-and-their-boundary-integrals).
 
 ## Scalar Field Operators
 
@@ -102,7 +102,7 @@ require a gradient operator should be used with H1.
 
 ### Square Operators
 
-These integrators are designed to be used with the BilinearForm object
+These integrators are designed to be used with a `BilinearForm` object
 to assemble square linear operators.
 
 | Class Name          | Spaces | Coef.| Operator                    | Continuous Op.          | Dimension  |
@@ -112,7 +112,7 @@ to assemble square linear operators.
 
 ### Mixed Operators
 
-These integrators are designed to be used with the MixedBilinearForm object to assemble square or rectangular linear operators.
+These integrators are designed to be used with a `MixedBilinearForm` object to assemble square or rectangular linear operators.
 
 | Class Name                            | Domain | Range  | Coef.   | Operator                                       | Continuous Op.                       | Dimension  |
 |---------------------------------------|--------|--------|:-------:|------------------------------------------------|--------------------------------------|:----------:|
@@ -150,7 +150,7 @@ functions but others require one or the other.
 
 ### Square Operators
 
-These integrators are designed to be used with the BilinearForm object
+These integrators are designed to be used with a `BilinearForm` object
 to assemble square linear operators.
 
 | Class Name             | Spaces | Coef.   | Operator                               | Continuous Op.                | Dimension  |
@@ -161,7 +161,7 @@ to assemble square linear operators.
 
 ### Mixed Operators
 
-These integrators are designed to be used with the MixedBilinearForm object to assemble square or rectangular linear operators.
+These integrators are designed to be used with a `MixedBilinearForm` object to assemble square or rectangular linear operators.
 
 | Class Name                           | Domain | Range  | Coef.   | Operator                                                 | Continuous Op.                                | Dimension  |
 |--------------------------------------|--------|--------|:-------:|----------------------------------------------------------|-----------------------------------------------|:----------:|
@@ -190,9 +190,9 @@ These integrators are designed to be used with the MixedBilinearForm object to a
 
 | Class Name                       | Domain | Range  | Coef. | Operator                                                               | Dimension | Notes |
 |----------------------------------|--------|--------|-------|------------------------------------------------------------------------|:---------:|-------|
-| VectorFEDivergenceIntegrator     | RT     | H1, L2 |   S   | $(\lambda\div\vec\{u}, v)$                                             | 2D, 3D    | Alternate implementation of MixedScalarDivergenceIntegrator. |
-| VectorFEWeakDivergenceIntegrator | ND     | H1     |   S   | $(-\lambda\vec\{u},\grad v)$                                           | 2D, 3D    | See MixedVectorWeakDivergenceIntegrator for a more general implementation. |
-| VectorFECurlIntegrator           | ND, RT | ND, RT |   S   | $(\lambda\curl\vec\{u},\vec\{v})$ or $(\lambda\vec\{u},\curl\vec\{v})$ | 3D        | If the domain is ND then the Curl operator is returned, if the range is ND then the weak Curl is returned, otherwise a failure is encountered. See MixedVectorCurlIntegrator and MixedVectorWeakCurlIntegrator for more general implementations. |
+| VectorFEDivergenceIntegrator     | RT     | H1, L2 |   S   | $(\lambda\div\vec\{u}, v)$                                             | 2D, 3D    | Alternate implementation of `MixedScalarDivergenceIntegrator`. |
+| VectorFEWeakDivergenceIntegrator | ND     | H1     |   S   | $(-\lambda\vec\{u},\grad v)$                                           | 2D, 3D    | See `MixedVectorWeakDivergenceIntegrator` for a more general implementation. |
+| VectorFECurlIntegrator           | ND, RT | ND, RT |   S   | $(\lambda\curl\vec\{u},\vec\{v})$ or $(\lambda\vec\{u},\curl\vec\{v})$ | 3D        | If the domain is ND then the Curl operator is returned, if the range is ND then the weak Curl is returned, otherwise a failure is encountered. See `MixedVectorCurlIntegrator` and `MixedVectorWeakCurlIntegrator` for more general implementations. |
 
 ## Vector Field Operators
 
@@ -299,9 +299,9 @@ situations rather than needing to reimplement their functionality.
 
 | Class Name          | Description                                                                                                            |
 |---------------------|------------------------------------------------------------------------------------------------------------------------|
-| TransposeIntegrator | Returns the transpose of the local matrix computed by another BilinearFormIntegrator                                   |
-| LumpedIntegrator    | Returns a diagonal local matrix where each entry is the sum of the corresponding row of a local matrix computed by another BilinearFormIntegrator (only implemented for square matrices) |
-| InverseIntegrator   | Returns the inverse of the local matrix computed by another BilinearFormIntegrator which produces a square local matrix |
+| TransposeIntegrator | Returns the transpose of the local matrix computed by another `BilinearFormIntegrator`                                 |
+| LumpedIntegrator    | Returns a diagonal local matrix where each entry is the sum of the corresponding row of a local matrix computed by another `BilinearFormIntegrator` (only implemented for square matrices) |
+| InverseIntegrator   | Returns the inverse of the local matrix computed by another `BilinearFormIntegrator` which produces a square local matrix |
 | SumIntegrator       | Returns the sum of a series of integrators with compatible dimensions (only implemented for square matrices)           |
 
 ## Weak Operators and Their Boundary Integrals

src/download.md

@@ -1,12 +1,12 @@
 ## Latest Release
 
-[New features](https://github.com/mfem/mfem/blob/v4.8/CHANGELOG)
+[New features](https://github.com/mfem/mfem/blob/v4.9/CHANGELOG)
 ┊ [Example codes](examples.md)
 ┊ [Code documentation](dox.md)
 ┊ [Sources](https://github.com/mfem/mfem)
 
 [<button type="button" class="btn btn-success">
-**Download mfem-4.8.tgz**
+**Download mfem-4.9.tgz**
 </button>](https://bit.ly/mfem-4-8)
 
 Please use the GitHub [issue tracker](https://github.com/mfem/mfem/issues)
@@ -17,6 +17,7 @@ or post [questions](https://github.com/mfem/mfem/issues/new?labels=question) or
 
  **Filename** | **Version** | **Release Date**  | **Size** | **[SLOC](https://github.com/AlDanial/cloc)** | **Notes** | **Documentation** |
  ------------ | ----------- | ----------------- | -------- | -------------------------------------------- | --------- | ----------------- |
+ [mfem-4.9.tgz](https://bit.ly/mfem-4-9) | v4.9 | Nov 2025 | X.XM | XXXK | | [docs/4.9](https://docs.mfem.org/4.9) |
  [mfem-4.8.tgz](https://bit.ly/mfem-4-8) | v4.8 | Apr 2025 | 4.1M | 458K | | [docs/4.8](https://docs.mfem.org/4.8) |
  [mfem-4.7.tgz](https://bit.ly/mfem-4-7) | v4.7 | May 2024 | 3.8M | 420K | | [docs/4.7](https://docs.mfem.org/4.7) |
  [mfem-4.6.tgz](https://bit.ly/mfem-4-6) | v4.6 | Sep 2023 | 3.6M | 397K | | [docs/4.6](https://docs.mfem.org/4.6) |

src/examples.md

@@ -2346,7 +2346,23 @@ _This miniapp has only a parallel ([contact-patch-test.cpp](https://github.com/m
 <div style="clear:both;"/></div>
 <br></div>
 
+<div id="contact" markdown="1">
+## Contact miniapp
 
+<a href="../img/gallery/workshop25/contact.mp4" target="_blank">
+<img class="floatright" width="250" src="../img/gallery/workshop25/contact.png">
+</a>
+
+This [miniapp](https://github.com/mfem/mfem/blob/master/miniapps/contact/contact.cpp) demonstrates how to use MFEM in combination with [Tribol](https://github.com/LLNL/Tribol)
+ to simulate frictionless contact between deformable solids. The contact gap function and its Jacobian are evaluated using Tribol’s mortar segment-to-segment method (see [Puso and Laursen, 2004](https://doi.org/10.1016/j.cma.2003.10.010)).  An **Interior-Point (IP)** optimization solver is used to solve an inequality-constrained minimization problem enforcing a non-penetration condition. Linear systems arising in the IP iterations are solved using **Conjugate Gradient (CG)** preconditioned with the [AMG with Filtering (AMGF)](https://arxiv.org/abs/2505.18576) solver.
+
+Several benchmark problems are provided, including the two-block compression, ironing and beam-sphere tests. These examples illustrate how MFEM and Tribol can be combined to build large-scale scalable simulations of contact mechanics problems.
+
+_This miniapp has only a parallel ([contact.cpp](https://github.com/mfem/mfem/blob/master/miniapps/contact/contact.cpp)) implementation. For more details, please see the documentation in [miniapps/contact/README.md](https://github.com/mfem/mfem/blob/master/miniapps/contact/README.md).
+**We recommend viewing Example 2 and the Tribol miniapp before viewing this miniapp.**_
+
+<div style="clear:both;"/></div>
+<br></div>
 
 <!-- ------------------------------------------------------------------------- -->
 
@@ -2492,6 +2508,7 @@ function update()
    + showElement("spde", (diffusion || nonlocal || stochastic) && h1 && galerkin && amg)
    + showElement("dpgminiapp", (diffusion || convectiondiffusion || maxwell || acoustics || wave) && (group2) && (dpg || amr) && (amg || ams || ads || pcg))
    + showElement("tribol", elasticity && h1 && galerkin && (superlu || minres || jacobi))
+   + showElement("contact", elasticity && h1 && galerkin && (pcg || amg))
 
    ; // ...end of expression
 

src/howto/findpts.md

@@ -8,7 +8,7 @@ a part of the gather-scatter library, [gslib](https://github.com/Nek5000/gslib).
 While `findpts` was originally developed for interpolation of grid functions in
 `H1` for meshes with quadrilateral or hexahedron elements, `FindPointsGSLIB`
 also enables interpolation of functions in `L2, H(div), H(curl)` on meshes with
-triangle and tetrahedral elements.
+triangle and tetrahedral elements. `FindPointsGSLIB` has also been extended to support surface meshes and effect the functionality on GPUs, and the core methodology is described in [Mittal et al., "General Field Evaluation in High-Order Meshes on GPUs"](https://www.sciencedirect.com/science/article/abs/pii/S0045793025000660).
 
 The key steps of using `FindPointsGSLIB`, as demonstrated in the
 [gslib](https://github.com/mfem/mfem/tree/master/miniapps/gslib) miniapps
@@ -29,7 +29,7 @@ These computational coordinates include the element number
 (e<sub>j</sub> in `mfem::Array<int> gsl_elem`) in which the point is found,
 the reference-space coordinates (**r**<sub>j</sub> in `mfem::Vector gsl_ref`) inside e<sub>j</sub>,
 and the MPI rank that the element is partitioned on (p<sub>j</sub> in `mfem::Array<int> gsl_proc`).
-`FindPoints` also returns a code (`mfem::Array<int> gsl_code`) to indicate weather
+`FindPoints(xyz)` also returns a code (`mfem::Array<int> gsl_code`) to indicate whether
 the point was found inside an element (`gsl_code[j] = 0`), on the edge/face of an
 element (`gsl_code[j] = 1`), or not found at all (`gsl_code[j] = 2`) for the case
 when the point is located outside the mesh.
@@ -48,9 +48,9 @@ we use `findpts` only for communicating computational coordinates of each point
 across MPI ranks, followed by MFEM's internal methods (`mfem::GridFunction::GetValues`)
 for interpolation.
 
-* **Note**, the `FindPointsGSLIB::FreeData()` method must be used before the
-program is terminated to free up the memory setup internally by `findpts` during
-the `setup` phase.
+    For custom interpolation (e.g., evaluating strain rate tensor), we have added functions that allow the user to first transfer computational coordinates for each point to ranks that overlap the corresponding point (`DistributePointInfoToOwningMPIRanks()`) and later transfer interpolated values back to ranks where the queries originated from (`DistributeInterpolatedValues()`).
+
+* **Note** All gslib allocations are freed if the `FindPointsGSLIB` object is destroyed before `MPI_Finalize()`. If the object might outlive `MPI_Finalize()`, use `FindPointsGSLIB::FreeData()` to avoid memory leaks.
 
 For convenience, `FindPointsGSLIB` class provides methods such as
 `FindPointsGSLIB::Interpolate(mesh, xyz, u, ui)` which combines the three steps
@@ -59,7 +59,7 @@ interpolation) into a single method. Please see the [class definition](https://g
 for more details.
 
 ## Application of FindPointsGSLIB
-The `gslib` miniapps demonstrate several application of `FindPointsGSLIB`:
+The `gslib` miniapps demonstrate several applications of `FindPointsGSLIB`:
 
 * [findpts/pfindpts](https://github.com/mfem/mfem/blob/master/miniapps/gslib/findpts.cpp)
 miniapps demonstrate high-order interpolation of a function in `H1`, `L2`, `H(div)`, or `H(curl)` at an
@@ -81,6 +81,8 @@ at the inter-domain boundaries.
 
 * [cht](https://github.com/mfem/mfem/blob/master/miniapps/navier/navier_cht.cpp)
 Navier miniapp demonstrates how a conjugate heat transfer problem can be
-solve with the incompressible Navier-Stokes equations and the unsteady heat
+solved with the incompressible Navier-Stokes equations and the unsteady heat
 equation solved on different grids. Here, `FindPointsGSLIB` is used to
 transfer the solution from one mesh to another to couple the two PDEs.
+
+`FindPointsGSLIB` is also used in TMOP-based r-adaptivity miniapps (e.g., [meshing/pmesh-optimizer](https://github.com/mfem/mfem/blob/master/miniapps/meshing/pmesh-optimizer.cpp)) to remap grid functions from the initial mesh to optimized mesh.

src/howto/howto-index.md

@@ -22,7 +22,7 @@ by our users in practice.
 ### Meshes
 - [Navigating Mesh Connectivity](nav-mesh-connectivity.md)
 - [Parallel Element Numbering](element-local-global-numbering.md)
-- [Finding Local Element Coordinates of Physical Points](findpts.md)
+- [Finding Local Element Coordinates of Physical Points with FindPointsGSLIB](findpts.md)
 - [Working with Nonconforming Meshes for AMR](ncmesh.md)
 
 ### Linear Algebra

src/index.md

@@ -92,22 +92,22 @@ See also our [Gallery](gallery.md), [Publications](publications.md), [Videos](vi
 
 Date         | Message
 ------------ | -----------------------------------------------------------------
+Nov XX, 2025 | Version 4.9 [released](https://github.com/mfem/mfem/blob/v4.9/CHANGELOG).
 Nov 13, 2025 | [Recap](https://computing.llnl.gov/about/newsroom/mfem-workshop-2025) of 2025 workshop.
-Oct 8, 2025 | [Postdoc position](https://jobs.smartrecruiters.com/LLNL/3743990009818416-casc-postdoctoral-research-staff-member) on the MFEM team. &nbsp;&nbsp;**[<i class="fa fa-envelope-o"></i>&nbsp;Apply](https://jobs.smartrecruiters.com/LLNL/3743990009818416-casc-postdoctoral-research-staff-member)**
+Oct 8, 2025  | [Postdoc position](https://jobs.smartrecruiters.com/LLNL/3743990009818416-casc-postdoctoral-research-staff-member) on the MFEM team. &nbsp;&nbsp;**[<i class="fa fa-envelope-o"></i>&nbsp;Apply](https://jobs.smartrecruiters.com/LLNL/3743990009818416-casc-postdoctoral-research-staff-member)**
 Sep 18, 2025 | [PyMFEM](https://github.com/mfem/PyMFEM) version 4.8 [released](https://github.com/mfem/PyMFEM/releases/tag/v_4.8.0).
 Sep 10, 2025 | 2025 [MFEM community workshop](workshop.md).
-Apr 9, 2025  | Version 4.8 [released](https://github.com/mfem/mfem/blob/v4.8/CHANGELOG).
 
 ## Latest Release
 
-[New features](https://github.com/mfem/mfem/blob/v4.8/CHANGELOG)
+[New features](https://github.com/mfem/mfem/blob/v4.9/CHANGELOG)
 ┊ [Examples](examples.md)
 ┊ [Code documentation](dox.md)
 ┊ [Sources](https://github.com/mfem/mfem)
 
 [<button type="button" class="btn btn-success">
-**Download mfem-4.8.tgz**
-</button>](https://bit.ly/mfem-4-8)
+**Download mfem-4.9.tgz**
+</button>](https://bit.ly/mfem-4-9)
 
 [Older releases](download.md) ┊ [Python wrapper](https://github.com/mfem/PyMFEM) ┊
 [![Launch Colab](https://colab.research.google.com/assets/colab-badge.svg){style="display:inline;margin:0"}](https://colab.research.google.com/github/GLVis/pyglvis/blob/main/examples/ex1.ipynb "Python Jupyter notebook")

src/news.md

@@ -2,6 +2,7 @@
 
 <img width=130px, style="margin:-10px"> | |
 ------------ | -----------------------------------------------------------------
+**Nov XX, 2025** | *Version 4.9 [released](https://github.com/mfem/mfem/blob/v4.9/CHANGELOG).*
 **Nov 13, 2025** | *[Recap](https://computing.llnl.gov/about/newsroom/mfem-workshop-2025) of the [2025 Workshop](workshop.md), held on September 10-11.*
 **Oct 8, 2025** | *[Postdoc position](https://jobs.smartrecruiters.com/LLNL/3743990009818416-casc-postdoctoral-research-staff-member) on the MFEM team. &nbsp;&nbsp;**[<i class="fa fa-envelope-o"></i>&nbsp;Apply](https://jobs.smartrecruiters.com/LLNL/3743990009818416-casc-postdoctoral-research-staff-member)***
 **Sep 18, 2025** | *[PyMFEM](https://github.com/mfem/PyMFEM) version 4.8 [released](https://github.com/mfem/PyMFEM/releases/tag/v_4.8.0)*