diff --git a/examples/SFMdata.h b/examples/SFMdata.h
index 04d3c9e47..3031828f1 100644
--- a/examples/SFMdata.h
+++ b/examples/SFMdata.h
@@ -22,6 +22,8 @@
  * Passing function argument allows to specificy an initial position, a pose increment and step count.
  */
 
+#pragma once
+
 // As this is a full 3D problem, we will use Pose3 variables to represent the camera
 // positions and Point3 variables (x, y, z) to represent the landmark coordinates.
 // Camera observations of landmarks (i.e. pixel coordinates) will be stored as Point2 (x, y).
@@ -66,4 +68,4 @@ std::vector<gtsam::Pose3> createPoses(
   }
 
   return poses;
-}
\ No newline at end of file
+}
diff --git a/gtsam/discrete/DiscreteJunctionTree.h b/gtsam/discrete/DiscreteJunctionTree.h
index f5bc9be1d..c74ad3cc2 100644
--- a/gtsam/discrete/DiscreteJunctionTree.h
+++ b/gtsam/discrete/DiscreteJunctionTree.h
@@ -16,6 +16,8 @@
  * @author Richard Roberts
  */
 
+#pragma once
+
 #include <gtsam/discrete/DiscreteFactorGraph.h>
 #include <gtsam/discrete/DiscreteBayesTree.h>
 #include <gtsam/inference/JunctionTree.h>
diff --git a/gtsam/geometry/Cyclic.h b/gtsam/geometry/Cyclic.h
index 35578ffe0..065cd0140 100644
--- a/gtsam/geometry/Cyclic.h
+++ b/gtsam/geometry/Cyclic.h
@@ -15,6 +15,8 @@
  * @author Frank Dellaert
  **/
 
+#pragma once
+
 #include <gtsam/base/Group.h>
 #include <gtsam/base/Testable.h>
 
diff --git a/gtsam/linear/GaussianJunctionTree.h b/gtsam/linear/GaussianJunctionTree.h
index 67e5a374b..c7f13ea5c 100644
--- a/gtsam/linear/GaussianJunctionTree.h
+++ b/gtsam/linear/GaussianJunctionTree.h
@@ -16,6 +16,8 @@
  * @author Richard Roberts
  */
 
+#pragma once
+
 #include <gtsam/linear/GaussianFactorGraph.h>
 #include <gtsam/linear/GaussianBayesTree.h>
 #include <gtsam/inference/JunctionTree.h>
diff --git a/gtsam/linear/linearAlgorithms-inst.h b/gtsam/linear/linearAlgorithms-inst.h
index 811e07121..d19ac6de5 100644
--- a/gtsam/linear/linearAlgorithms-inst.h
+++ b/gtsam/linear/linearAlgorithms-inst.h
@@ -15,6 +15,8 @@
  * @author  Richard Roberts
  */
 
+#pragma once
+
 #include <gtsam/linear/VectorValues.h>
 #include <gtsam/linear/GaussianConditional.h>
 #include <gtsam/base/treeTraversal-inst.h>
diff --git a/gtsam/linear/tests/powerMethodExample.h b/gtsam/linear/tests/powerMethodExample.h
index f80299386..994fcc640 100644
--- a/gtsam/linear/tests/powerMethodExample.h
+++ b/gtsam/linear/tests/powerMethodExample.h
@@ -19,6 +19,8 @@
  * PowerMethod/AcceleratedPowerMethod
  */
 
+#pragma once
+
 #include <gtsam/inference/Symbol.h>
 
 #include <iostream>
diff --git a/gtsam/navigation/ConstantVelocityFactor.h b/gtsam/navigation/ConstantVelocityFactor.h
index ed68ac077..6ab4c2f02 100644
--- a/gtsam/navigation/ConstantVelocityFactor.h
+++ b/gtsam/navigation/ConstantVelocityFactor.h
@@ -15,6 +15,8 @@
  * @author  Asa Hammond
  */
 
+#pragma once
+
 #include <gtsam/navigation/NavState.h>
 #include <gtsam/nonlinear/NonlinearFactor.h>
 
diff --git a/gtsam/navigation/MagFactor.h b/gtsam/navigation/MagFactor.h
index 74e9177d5..895ac6512 100644
--- a/gtsam/navigation/MagFactor.h
+++ b/gtsam/navigation/MagFactor.h
@@ -16,6 +16,8 @@
  * @date   January 29, 2014
  */
 
+#pragma once
+
 #include <gtsam/nonlinear/NonlinearFactor.h>
 #include <gtsam/geometry/Rot2.h>
 #include <gtsam/geometry/Rot3.h>
diff --git a/gtsam/nonlinear/WhiteNoiseFactor.h b/gtsam/nonlinear/WhiteNoiseFactor.h
index 95f46ab6c..1cd117437 100644
--- a/gtsam/nonlinear/WhiteNoiseFactor.h
+++ b/gtsam/nonlinear/WhiteNoiseFactor.h
@@ -17,6 +17,8 @@
  * @date   September 2011
  */
 
+#pragma once
+
 #include <gtsam/nonlinear/NonlinearFactor.h>
 #include <gtsam/linear/HessianFactor.h>
 #include <cmath>
diff --git a/gtsam/nonlinear/internal/LevenbergMarquardtState.h b/gtsam/nonlinear/internal/LevenbergMarquardtState.h
index cee839540..75e5a5135 100644
--- a/gtsam/nonlinear/internal/LevenbergMarquardtState.h
+++ b/gtsam/nonlinear/internal/LevenbergMarquardtState.h
@@ -16,6 +16,8 @@
  * @date    April 2016
  */
 
+#pragma once
+
 #include "NonlinearOptimizerState.h"
 
 #include <gtsam/nonlinear/Values.h>
diff --git a/gtsam/sfm/TranslationRecovery.h b/gtsam/sfm/TranslationRecovery.h
index c99836853..620a1c130 100644
--- a/gtsam/sfm/TranslationRecovery.h
+++ b/gtsam/sfm/TranslationRecovery.h
@@ -16,6 +16,8 @@
  * @brief Recovering translations in an epipolar graph when rotations are given.
  */
 
+#pragma once
+
 #include <map>
 #include <set>
 #include <utility>
diff --git a/gtsam/slam/TriangulationFactor.h b/gtsam/slam/TriangulationFactor.h
index 40e9538e2..b6da02d55 100644
--- a/gtsam/slam/TriangulationFactor.h
+++ b/gtsam/slam/TriangulationFactor.h
@@ -15,6 +15,8 @@
  * @author Frank Dellaert
  */
 
+#pragma once
+
 #include <gtsam/nonlinear/NonlinearFactor.h>
 #include <gtsam/geometry/CalibratedCamera.h>
 #include <boost/make_shared.hpp>
diff --git a/gtsam/symbolic/SymbolicJunctionTree.h b/gtsam/symbolic/SymbolicJunctionTree.h
index 7a152e532..0dcfae541 100644
--- a/gtsam/symbolic/SymbolicJunctionTree.h
+++ b/gtsam/symbolic/SymbolicJunctionTree.h
@@ -16,6 +16,8 @@
  * @author Richard Roberts
  */
 
+#pragma once
+
 #include <gtsam/symbolic/SymbolicFactorGraph.h>
 #include <gtsam/symbolic/SymbolicBayesTree.h>
 #include <gtsam/inference/JunctionTree.h>
diff --git a/gtsam_unstable/base/BTree.h b/gtsam_unstable/base/BTree.h
index 9d854a169..94e27d6c4 100644
--- a/gtsam_unstable/base/BTree.h
+++ b/gtsam_unstable/base/BTree.h
@@ -17,6 +17,8 @@
  * @date Feb 3, 2010
  */
 
+#pragma once
+
 #include <stack>
 #include <sstream>
 #include <boost/shared_ptr.hpp>
diff --git a/gtsam_unstable/base/Dummy.h b/gtsam_unstable/base/Dummy.h
index a2f544de5..548bce344 100644
--- a/gtsam_unstable/base/Dummy.h
+++ b/gtsam_unstable/base/Dummy.h
@@ -17,6 +17,8 @@
  * @date June 14, 2012
  */
 
+#pragma once
+
 #include <gtsam/global_includes.h>
 #include <gtsam_unstable/dllexport.h>
 #include <string>
diff --git a/gtsam_unstable/linear/ActiveSetSolver-inl.h b/gtsam_unstable/linear/ActiveSetSolver-inl.h
index 12374ac76..350985cf4 100644
--- a/gtsam_unstable/linear/ActiveSetSolver-inl.h
+++ b/gtsam_unstable/linear/ActiveSetSolver-inl.h
@@ -17,6 +17,8 @@
  * @date     2/11/16
  */
 
+#pragma once
+
 #include <gtsam_unstable/linear/InfeasibleInitialValues.h>
 
 /******************************************************************************/
@@ -283,4 +285,4 @@ Template std::pair<VectorValues, VectorValues> This::optimize() const {
 }
 
 #undef Template
-#undef This
\ No newline at end of file
+#undef This
diff --git a/gtsam_unstable/linear/LPSolver.h b/gtsam_unstable/linear/LPSolver.h
index 460b4b7ee..f36462bda 100644
--- a/gtsam_unstable/linear/LPSolver.h
+++ b/gtsam_unstable/linear/LPSolver.h
@@ -17,6 +17,8 @@
  * @date     6/16/16
  */
 
+#pragma once
+
 #include <gtsam_unstable/linear/LP.h>
 #include <gtsam_unstable/linear/ActiveSetSolver.h>
 #include <gtsam_unstable/linear/LPInitSolver.h>
diff --git a/gtsam_unstable/linear/QPSolver.h b/gtsam_unstable/linear/QPSolver.h
index 3854d2a15..ae87b3ab7 100644
--- a/gtsam_unstable/linear/QPSolver.h
+++ b/gtsam_unstable/linear/QPSolver.h
@@ -17,6 +17,8 @@
  * @date     6/16/16
  */
 
+#pragma once
+
 #include <gtsam_unstable/linear/QP.h>
 #include <gtsam_unstable/linear/ActiveSetSolver.h>
 #include <gtsam_unstable/linear/QPInitSolver.h>
@@ -45,4 +47,4 @@ struct QPPolicy {
 
 using QPSolver = ActiveSetSolver<QP, QPPolicy, QPInitSolver>;
 
-}
\ No newline at end of file
+}
diff --git a/gtsam_unstable/partition/FindSeparator.h b/gtsam_unstable/partition/FindSeparator.h
index 42d971a82..f4342695b 100644
--- a/gtsam_unstable/partition/FindSeparator.h
+++ b/gtsam_unstable/partition/FindSeparator.h
@@ -6,6 +6,8 @@
  *  Description: find the separator of bisectioning for a given graph
  */
 
+#pragma once
+
 #include <map>
 #include <vector>
 #include <boost/optional.hpp>
diff --git a/timing/timeSFMBAL.h b/timing/timeSFMBAL.h
index e24b50089..7af798887 100644
--- a/timing/timeSFMBAL.h
+++ b/timing/timeSFMBAL.h
@@ -16,6 +16,8 @@
  * @date    July 5, 2015
  */
 
+#pragma once
+
 #include <gtsam/sfm/SfmData.h>
 #include <gtsam/slam/dataset.h>
 #include <gtsam/nonlinear/LevenbergMarquardtOptimizer.h>