diff --git a/Apps/Common/SIMSemi3D.h b/Apps/Common/SIMSemi3D.h
index 34e26405..5e7eca8d 100644
--- a/Apps/Common/SIMSemi3D.h
+++ b/Apps/Common/SIMSemi3D.h
@@ -126,7 +126,7 @@ public:
   std::string getName() const override { return "Semi3D"; }
 
   //! \brief Adds fields to a data exporter.
-  void registerFields(DataExporter& exporter)
+  void registerFields(const DataExporter& exporter)
   {
     std::string name = exporter.getName();
     int plane = 1 + startCtx;
diff --git a/src/ASM/ASMs2DC1.C b/src/ASM/ASMs2DC1.C
index 4875268c..d986a77e 100644
--- a/src/ASM/ASMs2DC1.C
+++ b/src/ASM/ASMs2DC1.C
@@ -326,7 +326,7 @@ bool ASMs2DC1::addRigidCpl (int lindx, int ldim, int basis,
   \brief Computes coupling coefficients for a 2-master constraint in a C1-patch.
 */
 
-static void initMPC2 (MPC* mpc, std::vector<Vec3>& X)
+static void initMPC2 (MPC* mpc, const std::vector<Vec3>& X)
 {
   double s1 = (X[1]-X[0]).length();
   double s2 = (X[2]-X[0]).length();
@@ -339,7 +339,7 @@ static void initMPC2 (MPC* mpc, std::vector<Vec3>& X)
   \brief Computes coupling coefficients for a 3-master constraint in a C1-patch.
 */
 
-static void initMPC3 (MPC* mpc, std::vector<Vec3>& X)
+static void initMPC3 (MPC* mpc, const std::vector<Vec3>& X)
 {
   Vec3 V10(X[0]-X[1]);
   Vec3 V12(X[2]-X[1]);
diff --git a/src/ASM/DomainDecomposition.C b/src/ASM/DomainDecomposition.C
index 556f5aab..62e57504 100644
--- a/src/ASM/DomainDecomposition.C
+++ b/src/ASM/DomainDecomposition.C
@@ -207,7 +207,7 @@ static void getNumEdges(void* mesh, int sizeGID, int sizeLID,
                         int numCells, ZOLTAN_ID_PTR globalID,
                         ZOLTAN_ID_PTR localID, int* numEdges, int* err)
 {
-  IntMat& neigh = *static_cast<IntMat*>(mesh);
+  const IntMat& neigh = *static_cast<const IntMat*>(mesh);
   int* ne = numEdges;
   for (const std::vector<int>& n : neigh)
     *ne++ = std::accumulate(n.begin(), n.end(), 0,
@@ -225,7 +225,7 @@ static void getEdges(void* mesh, int sizeGID, int sizeLID, int numCells,
                      int* numEdges, ZOLTAN_ID_PTR nborGID, int* nborProc,
                      int wgtDim, float* egts, int* err)
 {
-  IntMat& neigh = *static_cast<IntMat*>(mesh);
+  const IntMat& neigh = *static_cast<const IntMat*>(mesh);
 
   memset(nborProc, 0, numCells*sizeof(int));
   for (const std::vector<int>& elm : neigh)
diff --git a/src/IFEM.C b/src/IFEM.C
index fa2fe7e8..5b7ee66b 100644
--- a/src/IFEM.C
+++ b/src/IFEM.C
@@ -42,7 +42,7 @@ int IFEM::Init (int arg_c, char** arg_v, const char* title)
 {
   argc = arg_c;
   argv = arg_v;
-  LinAlgInit& linalg = LinAlgInit::Init(argc,argv);
+  const LinAlgInit& linalg = LinAlgInit::Init(argc,argv);
   LinAlgInit::increfs();
 
   bool enableController = false;