luca-heltai · Oct 20, 2026 · Oct 20, 2026 · Oct 20, 2026 · Oct 20, 2026
diff --git a/.github/workflows/indentation.yml b/.github/workflows/indentation.yml
@@ -11,12 +11,12 @@ jobs:
 
     runs-on: ubuntu-latest
     container:
-      image: dealii/dealii:v9.6.0-jammy
+      image: dealii/dealii:v9.7.1-noble
       options: --user root --env OMPI_ALLOW_RUN_AS_ROOT=1 --env OMPI_ALLOW_RUN_AS_ROOT_CONFIRM=1
 
     steps:
     - uses: actions/checkout@v2
     - name: Check indentation
       run: |
         git config --global --add safe.directory /__w/pi-BEM/pi-BEM
-        ./scripts/check_indentation.sh 
+        ./scripts/check_indentation.sh
diff --git a/.github/workflows/tests.yml b/.github/workflows/tests.yml
@@ -12,7 +12,7 @@ jobs:
     runs-on: ubuntu-latest
 
     container:
-      image: dealii/dealii:v9.6.0-jammy
+      image: dealii/dealii:v9.7.1-noble
       options: --user root --env OMPI_ALLOW_RUN_AS_ROOT=1 --env OMPI_ALLOW_RUN_AS_ROOT_CONFIRM=1
 
     steps:
@@ -32,7 +32,7 @@ jobs:
     runs-on: ubuntu-latest
 
     container:
-      image: dealii/dealii:v9.6.0-jammy
+      image: dealii/dealii:v9.7.1-noble
       options: --user root --env OMPI_ALLOW_RUN_AS_ROOT=1 --env OMPI_ALLOW_RUN_AS_ROOT_CONFIRM=1
 
     steps:
@@ -45,4 +45,4 @@ jobs:
           cmake .. -GNinja -DCMAKE_CXX_FLAGS=-Werror -DCMAKE_BUILD_TYPE=Release
           ninja
           ctest -N
-          ctest --output-on-failure
+          ctest --output-on-failure
diff --git a/CMakeLists.txt b/CMakeLists.txt
@@ -39,7 +39,7 @@ ENDIF()
 CMAKE_MINIMUM_REQUIRED(VERSION 3.13.4)
 
 
-FIND_PACKAGE(deal.II 9.5.0 REQUIRED
+FIND_PACKAGE(deal.II 9.7.1 REQUIRED
   HINTS ${deal.II_DIR} ${DEAL_II_DIR} ../ ../../ $ENV{DEAL_II_DIR}
 )
 

diff --git a/README.md b/README.md
@@ -69,7 +69,7 @@ Just follow the [instructions](https://github.com/dealii/dealii/wiki/deal.II-in-
 We provide the possibility of using Docker as a tool to provide a fully operational environment for our library. To use such tool you need to install Docker following the [instructions](https://docs.docker.com/engine/installation/) provided by its authors. Then you can execute the following command line instruction
 
 ```
-docker run -v `pwd`:/pwd_to_your_own_directory/ -i -t dealii/dealii:v9.6.0-jammy bash
+docker run -v `pwd`:/pwd_to_your_own_directory/ -i -t dealii/dealii:v9.7.1-noble bash
 ```
 
 to retrieve the environment. In such a shell you can easily compile the pi-BEM library following its own instructions.

diff --git a/VERSION b/VERSION
@@ -1 +1 @@
-9.6.0
+9.7.1
diff --git a/docker/Dockerfile b/docker/Dockerfile
@@ -1,4 +1,4 @@
-FROM dealii/dealii:v9.6.0-jammy
+FROM dealii/dealii:v9.7.1-noble
 
 # pi-BEM master image
 RUN git clone https://github.com/mathLab/pi-BEM/ &&\

diff --git a/docs/complex-robin.md b/docs/complex-robin.md
@@ -1,16 +1,16 @@
 # Multiple potentials, complex problems and Robin conditions
 
-## Build and run with deal.II 9.6.0
+## Build and run with deal.II 9.7.1
 
 From the repository root, build both configurations in the supplied container:
 
 ```sh
 docker run --rm -v "$PWD:/workspace" -w /workspace \
-  dealii/dealii:v9.6.0-jammy bash -lc '
-    cmake -S . -B build-docker/9.6.0 -GNinja \
+  dealii/dealii:v9.7.1-noble bash -lc '
+    cmake -S . -B build-docker/9.7.1 -GNinja \
       -DCMAKE_BUILD_TYPE=DebugRelease -DCMAKE_CXX_FLAGS=-Werror &&
-    cmake --build build-docker/9.6.0 -j 4 &&
-    cd build-docker/9.6.0 && ctest --output-on-failure
+    cmake --build build-docker/9.7.1 -j 4 &&
+    cd build-docker/9.7.1 && ctest --output-on-failure
   '
 ```
 
@@ -24,13 +24,13 @@ The sphere examples use paths relative to `samples/`. Keep each run's output
 separate, and copy both the parameter file and its grid directory:
 
 ```sh
-mkdir -p build-docker/9.6.0/example-complex
-cp -R samples/grids build-docker/9.6.0/example-complex/
+mkdir -p build-docker/9.7.1/example-complex
+cp -R samples/grids build-docker/9.7.1/example-complex/
 cp samples/ref_params_sphere_complex_DR_direct.prm \
-  build-docker/9.6.0/example-complex/parameters_bem_3.prm
+  build-docker/9.7.1/example-complex/parameters_bem_3.prm
 docker run --rm -v "$PWD:/workspace" \
-  -w /workspace/build-docker/9.6.0/example-complex \
-  dealii/dealii:v9.6.0-jammy bash -lc '
+  -w /workspace/build-docker/9.7.1/example-complex \
+  dealii/dealii:v9.7.1-noble bash -lc '
     mpirun -np 2 ../bem_fma_3d 2
   '
 ```

diff --git a/include/ass_leg_function.h b/include/ass_leg_function.h
@@ -3,10 +3,10 @@
 #define __sissa__ass_leg_function_h
 
 #include <deal.II/base/convergence_table.h>
+#include <deal.II/base/observer_pointer.h>
 #include <deal.II/base/parsed_function.h>
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/base/quadrature_selector.h>
-#include <deal.II/base/smartpointer.h>
 #include <deal.II/base/utilities.h>
 
 #include <deal.II/dofs/dof_accessor.h>

diff --git a/include/bem_fma.h b/include/bem_fma.h
@@ -272,8 +272,8 @@ class BEMFMA : public ParameterAcceptor
 protected:
   /// Three pointers to the problem parameters to be set equal to
   /// the ones in the calling problem through the init function
-  SmartPointer<const DoFHandler<dim - 1, dim>> fma_dh;
-  SmartPointer<const Mapping<dim - 1, dim>>    fma_mapping;
+  ObserverPointer<const DoFHandler<dim - 1, dim>> fma_dh;
+  ObserverPointer<const Mapping<dim - 1, dim>>    fma_mapping;
 
   /// Truncation order for the multipole
   /// and local expansion series: it is
@@ -455,8 +455,8 @@ class BEMFMA : public ParameterAcceptor
   ConditionalOStream pcout;
 
   /// TODO parsed quadrature?
-  std::shared_ptr<Quadrature<dim - 1>> quadrature;
-  SmartPointer<const Vector<double>>   dirichlet_nodes;
+  std::shared_ptr<Quadrature<dim - 1>>  quadrature;
+  ObserverPointer<const Vector<double>> dirichlet_nodes;
   /// This should be erased by the usage of the constraint matrix.
   const std::vector<std::set<types::global_dof_index>> *double_nodes_set;
 

diff --git a/include/bem_problem.h b/include/bem_problem.h
@@ -22,11 +22,11 @@
 
 #include <deal.II/base/conditional_ostream.h>
 #include <deal.II/base/convergence_table.h>
+#include <deal.II/base/observer_pointer.h>
 #include <deal.II/base/parameter_acceptor.h>
 #include <deal.II/base/parsed_function.h>
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/base/quadrature_selector.h>
-#include <deal.II/base/smartpointer.h>
 #include <deal.II/base/types.h>
 #include <deal.II/base/utilities.h>
 

diff --git a/include/boundary_conditions.h b/include/boundary_conditions.h
@@ -15,11 +15,11 @@
 #ifndef boundary_conditions_h
 #define boundary_conditions_h
 #include <deal.II/base/convergence_table.h>
+#include <deal.II/base/observer_pointer.h>
 #include <deal.II/base/parameter_acceptor.h>
 #include <deal.II/base/parsed_function.h>
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/base/quadrature_selector.h>
-#include <deal.II/base/smartpointer.h>
 #include <deal.II/base/types.h>
 #include <deal.II/base/utilities.h>
 

diff --git a/include/computational_domain.h b/include/computational_domain.h
@@ -20,11 +20,11 @@
 
 #include <deal.II/base/conditional_ostream.h>
 #include <deal.II/base/convergence_table.h>
+#include <deal.II/base/observer_pointer.h>
 #include <deal.II/base/parameter_acceptor.h>
 #include <deal.II/base/parsed_function.h>
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/base/quadrature_selector.h>
-#include <deal.II/base/smartpointer.h>
 #include <deal.II/base/types.h>
 #include <deal.II/base/utilities.h>
 

diff --git a/include/constrained_matrix.h b/include/constrained_matrix.h
@@ -18,7 +18,7 @@
 #include <deal.II/base/config.h>
 
 #include <deal.II/base/memory_consumption.h>
-#include <deal.II/base/smartpointer.h>
+#include <deal.II/base/observer_pointer.h>
 #include <deal.II/base/thread_management.h>
 #include <deal.II/base/types.h>
 
@@ -76,7 +76,7 @@ class ConstrainedOperator
   unsigned int                     this_mpi_process;
 };
 
-/*@}*/
+/** @} */
 /*---------------------- Inline functions -----------------------------------*/
 
 //--------------------------------Iterators--------------------------------------//

diff --git a/include/constrained_matrix_complex.h b/include/constrained_matrix_complex.h
@@ -18,7 +18,7 @@
 #include <deal.II/base/config.h>
 
 #include <deal.II/base/memory_consumption.h>
-#include <deal.II/base/smartpointer.h>
+#include <deal.II/base/observer_pointer.h>
 #include <deal.II/base/thread_management.h>
 #include <deal.II/base/types.h>
 

diff --git a/include/driver.h b/include/driver.h
@@ -17,11 +17,11 @@
 
 #include <deal.II/base/conditional_ostream.h>
 #include <deal.II/base/convergence_table.h>
+#include <deal.II/base/observer_pointer.h>
 #include <deal.II/base/parameter_acceptor.h>
 #include <deal.II/base/parsed_function.h>
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/base/quadrature_selector.h>
-#include <deal.II/base/smartpointer.h>
 #include <deal.II/base/utilities.h>
 
 #include <deal.II/lac/full_matrix.h>

diff --git a/include/laplace_kernel.h b/include/laplace_kernel.h
@@ -1,8 +1,8 @@
 #ifndef laplace_kernel_h
 #define laplace_kernel_h
 
+#include <deal.II/base/observer_pointer.h>
 #include <deal.II/base/point.h>
-#include <deal.II/base/smartpointer.h>
 #include <deal.II/base/utilities.h>
 
 #include <cmath>

diff --git a/include/local_expansion.h b/include/local_expansion.h
@@ -146,8 +146,9 @@ class LocalExpansion
 
   /// Return a read-only reference to owned triangular coefficient storage,
   /// indexed by n*(n+1)/2+m. Copy the vector to retain an independent snapshot.
-  /// Use SetCoeff()/AddToCoeff() for changes; assignment invalidates
-  /// references.
+  /// Use SetCoeff()/AddToCoeff() for changes. Assignment can invalidate
+  /// references to individual coefficients; the vector object remains owned
+  /// here.
   inline const std::vector<std::complex<double>> &
   GetCoeffs() const
   {

diff --git a/include/multipole_expansion.h b/include/multipole_expansion.h
@@ -141,8 +141,9 @@ class MultipoleExpansion
 
   /// Return a read-only reference to owned triangular coefficient storage,
   /// indexed by n*(n+1)/2+m. Copy the vector to retain an independent snapshot.
-  /// Use SetCoeff()/AddToCoeff() for changes; assignment invalidates
-  /// references.
+  /// Use SetCoeff()/AddToCoeff() for changes. Assignment can invalidate
+  /// references to individual coefficients; the vector object remains owned
+  /// here.
   inline const std::vector<std::complex<double>> &
   GetCoeffs() const
   {

diff --git a/include/octree_block.h b/include/octree_block.h
@@ -16,11 +16,11 @@
 #define octree_block_h
 
 #include <deal.II/base/convergence_table.h>
+#include <deal.II/base/observer_pointer.h>
 #include <deal.II/base/parsed_function.h>
 #include <deal.II/base/point.h>
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/base/quadrature_selector.h>
-#include <deal.II/base/smartpointer.h>
 #include <deal.II/base/utilities.h>
 
 #include <deal.II/dofs/dof_accessor.h>

diff --git a/include/true_constrained_matrix.h b/include/true_constrained_matrix.h
@@ -19,7 +19,7 @@
 #  include <deal.II/base/config.h>
 
 #  include <deal.II/base/memory_consumption.h>
-#  include <deal.II/base/smartpointer.h>
+#  include <deal.II/base/observer_pointer.h>
 #  include <deal.II/base/thread_management.h>
 
 #  include <deal.II/lac/affine_constraints.h>
@@ -90,7 +90,7 @@ class ConstrainedOperator
   unsigned int                     this_mpi_process;
 };
 
-/*@}*/
+/** @} */
 /*---------------------- Inline functions -----------------------------------*/
 
 

diff --git a/include/vector_tools_integrate_difference.h b/include/vector_tools_integrate_difference.h
@@ -2,10 +2,10 @@
 #define vector_tools_integrate_difference_h
 
 #include <deal.II/base/convergence_table.h>
+#include <deal.II/base/observer_pointer.h>
 #include <deal.II/base/parsed_function.h>
 #include <deal.II/base/quadrature_lib.h>
 #include <deal.II/base/quadrature_selector.h>
-#include <deal.II/base/smartpointer.h>
 #include <deal.II/base/types.h>
 #include <deal.II/base/utilities.h>
 

diff --git a/scripts/docker.sh b/scripts/docker.sh
@@ -1 +1,3 @@
-docker run -P -t -i -v `pwd`:/home/dealii/app:rw mathlab/deal2lkit:v8.5.0-debugrelease  
+#!/bin/sh
+# Open the supported deal.II environment in the current checkout.
+exec docker run --rm -it -v "$PWD:/app" -w /app dealii/dealii:v9.7.1-noble bash
diff --git a/scripts/dr.sh b/scripts/dr.sh
@@ -1,2 +1,4 @@
-#!/bin/bash
-docker run  --user $(id -u):$(id -g) -i -t --rm -P -v `pwd`:/app:rw dealii/dealii:v9.6.0-jammy  /bin/sh -c "cd /app; $@"
+#!/bin/sh
+# Run a shell command in the supported deal.II environment.
+exec docker run --user "$(id -u):$(id -g)" -it --rm \
+  -v "$PWD:/app" -w /app dealii/dealii:v9.7.1-noble /bin/sh -c "$*"
diff --git a/scripts/run_docker.sh b/scripts/run_docker.sh
@@ -1,3 +1,4 @@
-#!/bin/bash
-CMD=`grep -e '- ' .travis.yml | sed 's/- //'`
-docker run  --user $(id -u):$(id -g) -i -t --rm -P -v `pwd`:/app:rw mathlab/deal2lkit:v9.1.1-debugrelease  /bin/sh -c "cd /app; $CMD"
+#!/bin/sh
+# Run the same warning-free Debug/Release checks used for releases.
+exec docker run --rm -v "$PWD:/app" -w /app dealii/dealii:v9.7.1-noble \
+  bash -lc 'set -e; cmake -S . -B build-docker/9.7.1 -GNinja -DCMAKE_BUILD_TYPE=DebugRelease -DCMAKE_CXX_FLAGS=-Werror; cmake --build build-docker/9.7.1 -j 4; cd build-docker/9.7.1; ctest --output-on-failure'