diff --git a/limap/base/camera.h b/limap/base/camera.h
index 1cc80260..88dbf53e 100644
--- a/limap/base/camera.h
+++ b/limap/base/camera.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_BASE_CAMERA_H_
-#define LIMAP_BASE_CAMERA_H_
+#pragma once
 
 #include <cmath>
 #include <fstream>
@@ -114,5 +113,3 @@ class CameraPose {
 };
 
 } // namespace limap
-
-#endif
diff --git a/limap/base/camera_models.h b/limap/base/camera_models.h
index 4c460226..89ca759f 100644
--- a/limap/base/camera_models.h
+++ b/limap/base/camera_models.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_BASE_CAMERA_MODELS_H_
-#define LIMAP_BASE_CAMERA_MODELS_H_
+#pragma once
 
 #include <colmap/sensor/models.h>
 
@@ -45,5 +44,3 @@ void ParamsToKvec(const colmap::CameraModelId model_id, const T *params,
 }
 
 } // namespace limap
-
-#endif
diff --git a/limap/base/camera_view.h b/limap/base/camera_view.h
index c1dbff37..466b20e4 100644
--- a/limap/base/camera_view.h
+++ b/limap/base/camera_view.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_BASE_CAMERA_VIEW_H_
-#define LIMAP_BASE_CAMERA_VIEW_H_
+#pragma once
 
 #include <cmath>
 #include <fstream>
@@ -89,5 +88,3 @@ class CameraView : public CameraImage {
 };
 
 } // namespace limap
-
-#endif
diff --git a/limap/base/graph.h b/limap/base/graph.h
index b73eae2a..6d73dfe1 100644
--- a/limap/base/graph.h
+++ b/limap/base/graph.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_BASE_GRAPH_H_
-#define LIMAP_BASE_GRAPH_H_
+#pragma once
 
 #include <algorithm>
 #include <chrono>
@@ -97,5 +96,3 @@ CountTrackEdges(const Graph &graph, const std::vector<size_t> &track_labels,
                 const std::vector<bool> is_root);
 
 } // namespace limap
-
-#endif
diff --git a/limap/base/image_collection.h b/limap/base/image_collection.h
index 3284804d..a26fe458 100644
--- a/limap/base/image_collection.h
+++ b/limap/base/image_collection.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_BASE_IMAGE_COLLECTION_H_
-#define LIMAP_BASE_IMAGE_COLLECTION_H_
+#pragma once
 
 #include <cmath>
 #include <fstream>
@@ -115,5 +114,3 @@ class ImageCollection {
 };
 
 } // namespace limap
-
-#endif
diff --git a/limap/base/infinite_line.h b/limap/base/infinite_line.h
index e5eaba0e..12c7eed7 100644
--- a/limap/base/infinite_line.h
+++ b/limap/base/infinite_line.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_BASE_INFINITE_LINE_H_
-#define LIMAP_BASE_INFINITE_LINE_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -86,5 +85,3 @@ Line3d GetLineSegmentFromInfiniteLine3d(const InfiniteLine3d &inf_line,
                                         const int num_outliers = 2);
 
 } // namespace limap
-
-#endif
diff --git a/limap/base/line_dists.h b/limap/base/line_dists.h
index d4b5c7ad..78d4b9cb 100644
--- a/limap/base/line_dists.h
+++ b/limap/base/line_dists.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_BASE_LINE_DISTS_H_
-#define LIMAP_BASE_LINE_DISTS_H_
+#pragma once
 
 #include <cmath>
 #include <pybind11/eigen.h>
@@ -221,5 +220,3 @@ double infinite_dist_perpendicular_scaleinv_line3dpp(const Line3d &l1,
                                                      const Line3d &l2);
 
 } // namespace limap
-
-#endif
diff --git a/limap/base/line_linker.h b/limap/base/line_linker.h
index ce167d4d..e10f3258 100644
--- a/limap/base/line_linker.h
+++ b/limap/base/line_linker.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_BASE_LINE_LINKER_H_
-#define LIMAP_BASE_LINE_LINKER_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -228,5 +227,3 @@ class LineLinker {
 };
 
 } // namespace limap
-
-#endif
diff --git a/limap/base/linebase.h b/limap/base/linebase.h
index f0ae238a..3a6683cc 100644
--- a/limap/base/linebase.h
+++ b/limap/base/linebase.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_BASE_LINEBASE_H_
-#define LIMAP_BASE_LINEBASE_H_
+#pragma once
 
 #include <cmath>
 #include <pybind11/eigen.h>
@@ -69,5 +68,3 @@ Line3d unprojection_line2d(const Line2d &line2d, const CameraView &view,
                            const std::pair<double, double> &depths);
 
 } // namespace limap
-
-#endif
diff --git a/limap/base/linetrack.h b/limap/base/linetrack.h
index 506d9f9d..ab4aabed 100644
--- a/limap/base/linetrack.h
+++ b/limap/base/linetrack.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_BASE_LINETRACK_H_
-#define LIMAP_BASE_LINETRACK_H_
+#pragma once
 
 #include <cmath>
 #include <map>
@@ -82,5 +81,3 @@ void ComputeFConsistencySamples(
     const std::pair<double, double> sample_range, const int n_samples);
 
 } // namespace limap
-
-#endif
diff --git a/limap/base/pointtrack.h b/limap/base/pointtrack.h
index 70f0be58..30e51dc5 100644
--- a/limap/base/pointtrack.h
+++ b/limap/base/pointtrack.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_BASE_POINTTRACK_H_
-#define LIMAP_BASE_POINTTRACK_H_
+#pragma once
 
 #include <cmath>
 #include <map>
@@ -53,5 +52,3 @@ class PointTrack {
 };
 
 } // namespace limap
-
-#endif
diff --git a/limap/base/pose.h b/limap/base/pose.h
index 56deb01f..e84a57bd 100644
--- a/limap/base/pose.h
+++ b/limap/base/pose.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_BASE_POSE_H_
-#define LIMAP_BASE_POSE_H_
+#pragma once
 
 #include <colmap/geometry/pose.h>
 
@@ -30,5 +29,3 @@ Eigen::Matrix3d QuaternionToRotationMatrix(const Eigen::Vector4d &qvec);
 Eigen::Vector4d NormalizeQuaternion(const Eigen::Vector4d &qvec);
 
 } // namespace limap
-
-#endif
diff --git a/limap/base/transforms.h b/limap/base/transforms.h
index 6ee4ea7d..68254821 100644
--- a/limap/base/transforms.h
+++ b/limap/base/transforms.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_BASE_TRANSFORMS_H_
-#define LIMAP_BASE_TRANSFORMS_H_
+#pragma once
 
 #include "base/camera.h"
 #include "base/pose.h"
@@ -30,5 +29,3 @@ CameraPose pose_similarity_transform(const CameraPose &pose,
                                      const SimilarityTransform3 &transform);
 
 } // namespace limap
-
-#endif
diff --git a/limap/ceresbase/ceres_extensions.h b/limap/ceresbase/ceres_extensions.h
index 852aaa84..e2bf22c3 100644
--- a/limap/ceresbase/ceres_extensions.h
+++ b/limap/ceresbase/ceres_extensions.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_CERESBASE_CERES_EXTENSIONS_H
-#define LIMAP_CERESBASE_CERES_EXTENSIONS_H
+#pragma once
 
 // Modified from the pixel-perfect-sfm project
 
@@ -49,5 +48,3 @@ inline void AssignSolverOptionsFromDict(ceres::Solver::Options &solver_options,
   ASSIGN_PYDICT_ITEM_TO_MEMBER(solver_options, dict, inner_iteration_tolerance,
                                double)
 }
-
-#endif
diff --git a/limap/ceresbase/interpolation.h b/limap/ceresbase/interpolation.h
index 8a9992e3..38820569 100644
--- a/limap/ceresbase/interpolation.h
+++ b/limap/ceresbase/interpolation.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_CERESBASE_INTERPOLATION_H
-#define LIMAP_CERESBASE_INTERPOLATION_H
+#pragma once
 
 // Modified from the pixel-perfect-sfm project
 // [LINK] https://github.com/cvg/pixel-perfect-sfm
@@ -745,5 +744,3 @@ template <typename Grid> class PixelInterpolator : public Interpolator<Grid> {
 };
 
 } // namespace limap
-
-#endif
diff --git a/limap/ceresbase/line_dists.h b/limap/ceresbase/line_dists.h
index 40ef3a5e..e46602c0 100644
--- a/limap/ceresbase/line_dists.h
+++ b/limap/ceresbase/line_dists.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_CERESBASE_LINE_DISTS_H
-#define LIMAP_CERESBASE_LINE_DISTS_H
+#pragma once
 
 #include <ceres/ceres.h>
 #include <ceres/rotation.h>
@@ -88,5 +87,3 @@ T CeresComputeDist3D_angle(const T dir1[3], const T dir2[3]) {
 }
 
 } // namespace limap
-
-#endif
diff --git a/limap/ceresbase/line_projection.h b/limap/ceresbase/line_projection.h
index 563a02c6..ae6f3149 100644
--- a/limap/ceresbase/line_projection.h
+++ b/limap/ceresbase/line_projection.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_CERESBASE_LINE_PROJECTION_H
-#define LIMAP_CERESBASE_LINE_PROJECTION_H
+#pragma once
 
 #include "ceresbase/line_transforms.h"
 #include "ceresbase/point_projection.h"
@@ -214,5 +213,3 @@ bool Ceres_GetIntersection2dFromInfiniteLine3d(
 }
 
 } // namespace limap
-
-#endif
diff --git a/limap/ceresbase/line_transforms.h b/limap/ceresbase/line_transforms.h
index b5a6be3d..a04b1425 100644
--- a/limap/ceresbase/line_transforms.h
+++ b/limap/ceresbase/line_transforms.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_CERESBASE_LINE_TRANSFORMS_H
-#define LIMAP_CERESBASE_LINE_TRANSFORMS_H
+#pragma once
 
 #include <ceres/ceres.h>
 #include <ceres/rotation.h>
@@ -73,5 +72,3 @@ bool Ceres_IntersectLineCoordinates(const T coor1[3], const T coor2[3],
 }
 
 } // namespace limap
-
-#endif
diff --git a/limap/ceresbase/loss_functions.h b/limap/ceresbase/loss_functions.h
index 34d51d1a..02edb868 100644
--- a/limap/ceresbase/loss_functions.h
+++ b/limap/ceresbase/loss_functions.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_CERESBASE_LOSS_FUNCTIONS_H
-#define LIMAP_CERESBASE_LOSS_FUNCTIONS_H
+#pragma once
 
 // Modified from the pixel-perfect-sfm project.
 
@@ -249,5 +248,3 @@ inline ceres::LossFunction *CreateLossFunction(py::dict dict) {
 }
 
 } // namespace limap
-
-#endif
diff --git a/limap/ceresbase/parameterization.h b/limap/ceresbase/parameterization.h
index 39b5b58e..ea5a8b74 100644
--- a/limap/ceresbase/parameterization.h
+++ b/limap/ceresbase/parameterization.h
@@ -1,7 +1,4 @@
-#ifndef LIMAP_CERESBASE_PARAMETERIZATION_H
-#define LIMAP_CERESBASE_PARAMETERIZATION_H
-
-// Inspired by COLMAP: https://colmap.github.io/
+#pragma once
 
 #include <ceres/ceres.h>
 #include <ceres/rotation.h>
@@ -38,5 +35,3 @@ inline void SetSphereManifold(ceres::Problem *problem, double *params) {
 }
 
 } // namespace limap
-
-#endif
diff --git a/limap/ceresbase/point_projection.h b/limap/ceresbase/point_projection.h
index 908f7e13..d952488d 100644
--- a/limap/ceresbase/point_projection.h
+++ b/limap/ceresbase/point_projection.h
@@ -1,7 +1,4 @@
-#ifndef LIMAP_CERESBASE_POINT_PROJECTION_H
-#define LIMAP_CERESBASE_POINT_PROJECTION_H
-
-// Modified from the pixel-perfect-sfm project
+#pragma once
 
 #include <ceres/ceres.h>
 #include <ceres/rotation.h>
@@ -64,5 +61,3 @@ template <typename T> inline bool IsInsideZeroL(const T &value, double L) {
 }
 
 } // namespace limap
-
-#endif
diff --git a/limap/estimators/absolute_pose/hybrid_pose_estimator.h b/limap/estimators/absolute_pose/hybrid_pose_estimator.h
index 45f1b7ee..a3491838 100644
--- a/limap/estimators/absolute_pose/hybrid_pose_estimator.h
+++ b/limap/estimators/absolute_pose/hybrid_pose_estimator.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_ESTIMATORS_POSE_HYBRID_POSE_ESTIMATOR_H_
-#define LIMAP_ESTIMATORS_POSE_HYBRID_POSE_ESTIMATOR_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include "base/camera.h"
@@ -139,5 +138,3 @@ EstimateAbsolutePose_PointLine_Hybrid(
 } // namespace estimators
 
 } // namespace limap
-
-#endif
diff --git a/limap/estimators/absolute_pose/joint_pose_estimator.h b/limap/estimators/absolute_pose/joint_pose_estimator.h
index 2aa41987..7275df93 100644
--- a/limap/estimators/absolute_pose/joint_pose_estimator.h
+++ b/limap/estimators/absolute_pose/joint_pose_estimator.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_ESTIMATORS_POSE_JOINT_POSE_ESTIMATOR_H_
-#define LIMAP_ESTIMATORS_POSE_JOINT_POSE_ESTIMATOR_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include "base/camera.h"
@@ -105,5 +104,3 @@ EstimateAbsolutePose_PointLine(const std::vector<Line3d> &l3ds,
 } // namespace estimators
 
 } // namespace limap
-
-#endif
diff --git a/limap/estimators/absolute_pose/pl_absolute_pose_hybrid_ransac.h b/limap/estimators/absolute_pose/pl_absolute_pose_hybrid_ransac.h
index b16c4a8c..1d15b564 100644
--- a/limap/estimators/absolute_pose/pl_absolute_pose_hybrid_ransac.h
+++ b/limap/estimators/absolute_pose/pl_absolute_pose_hybrid_ransac.h
@@ -28,8 +28,7 @@
 //
 // author: Torsten Sattler, torsten.sattler.de@googlemail.com
 
-#ifndef LIMAP_ESTIMATORS_POSE_PL_ABSOLUTE_POSE_HYBRID_RANSAC_H_
-#define LIMAP_ESTIMATORS_POSE_PL_ABSOLUTE_POSE_HYBRID_RANSAC_H_
+#pragma once
 
 #include <algorithm>
 #include <cmath>
@@ -622,5 +621,3 @@ class PointLineAbsolutePoseHybridRansac : public HybridRansacBase {
 } // namespace estimators
 
 } // namespace limap
-
-#endif
diff --git a/limap/estimators/absolute_pose/pl_absolute_pose_ransac.h b/limap/estimators/absolute_pose/pl_absolute_pose_ransac.h
index 0f93c8de..25c61bcd 100644
--- a/limap/estimators/absolute_pose/pl_absolute_pose_ransac.h
+++ b/limap/estimators/absolute_pose/pl_absolute_pose_ransac.h
@@ -1,4 +1,3 @@
-// Copyright (c) 2019, Torsten Sattler
 // All rights reserved.
 //
 // Redistribution and use in source and binary forms, with or without
@@ -28,8 +27,7 @@
 //
 // author: Torsten Sattler, torsten.sattler.de@googlemail.com
 
-#ifndef LIMAP_ESTIMATORS_POSE_PL_ABSOLUTE_POSE_RANSAC_H_
-#define LIMAP_ESTIMATORS_POSE_PL_ABSOLUTE_POSE_RANSAC_H_
+#pragma once
 
 #include <algorithm>
 #include <cmath>
@@ -457,5 +455,3 @@ class PointLineAbsolutePoseRansac : public RansacBase {
 } // namespace estimators
 
 } // namespace limap
-
-#endif
diff --git a/limap/estimators/extended_hybrid_ransac.h b/limap/estimators/extended_hybrid_ransac.h
index 31588b34..2556fbff 100644
--- a/limap/estimators/extended_hybrid_ransac.h
+++ b/limap/estimators/extended_hybrid_ransac.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_ESTIMATORS_EXTENDED_HYBRID_RANSAC_H_
-#define LIMAP_ESTIMATORS_EXTENDED_HYBRID_RANSAC_H_
+#pragma once
 
 #include <RansacLib/hybrid_ransac.h>
 
@@ -18,5 +17,3 @@ class ExtendedHybridLORansacOptions : public ransac_lib::HybridLORansacOptions {
 } // namespace estimators
 
 } // namespace limap
-
-#endif
\ No newline at end of file
diff --git a/limap/evaluation/base_evaluator.h b/limap/evaluation/base_evaluator.h
index 8d8d4d30..528bc460 100644
--- a/limap/evaluation/base_evaluator.h
+++ b/limap/evaluation/base_evaluator.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_EVALUATION_BASE_EVALUATOR_H_
-#define LIMAP_EVALUATION_BASE_EVALUATOR_H_
+#pragma once
 
 #include "base/linebase.h"
 #include "util/types.h"
@@ -37,5 +36,3 @@ class BaseEvaluator {
 } // namespace evaluation
 
 } // namespace limap
-
-#endif
diff --git a/limap/evaluation/mesh_evaluator.h b/limap/evaluation/mesh_evaluator.h
index 151afe1e..38a97272 100644
--- a/limap/evaluation/mesh_evaluator.h
+++ b/limap/evaluation/mesh_evaluator.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_EVALUATION_MESH_EVALUATOR_H_
-#define LIMAP_EVALUATION_MESH_EVALUATOR_H_
+#pragma once
 
 #include "base/linebase.h"
 #include "evaluation/base_evaluator.h"
@@ -30,5 +29,3 @@ class MeshEvaluator : public BaseEvaluator {
 } // namespace evaluation
 
 } // namespace limap
-
-#endif
diff --git a/limap/evaluation/point_cloud_evaluator.h b/limap/evaluation/point_cloud_evaluator.h
index 84b67406..cb2d0047 100644
--- a/limap/evaluation/point_cloud_evaluator.h
+++ b/limap/evaluation/point_cloud_evaluator.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_EVALUATION_POINT_CLOUD_EVALUATOR_H_
-#define LIMAP_EVALUATION_POINT_CLOUD_EVALUATOR_H_
+#pragma once
 
 #include "base/linebase.h"
 #include "evaluation/base_evaluator.h"
@@ -47,5 +46,3 @@ class PointCloudEvaluator : public BaseEvaluator {
 } // namespace evaluation
 
 } // namespace limap
-
-#endif
diff --git a/limap/evaluation/refline_evaluator.h b/limap/evaluation/refline_evaluator.h
index 1fca0978..b1ba9373 100644
--- a/limap/evaluation/refline_evaluator.h
+++ b/limap/evaluation/refline_evaluator.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_EVALUATION_REFLINE_EVALUATOR_H_
-#define LIMAP_EVALUATION_REFLINE_EVALUATOR_H_
+#pragma once
 
 #include "base/linebase.h"
 #include "util/types.h"
@@ -39,5 +38,3 @@ class RefLineEvaluator {
 } // namespace evaluation
 
 } // namespace limap
-
-#endif
diff --git a/limap/features/dense_sift.h b/limap/features/dense_sift.h
index 4c2f27e9..53359eb3 100644
--- a/limap/features/dense_sift.h
+++ b/limap/features/dense_sift.h
@@ -1,3 +1,4 @@
+#pragma once
 
 #include <pybind11/numpy.h>
 #include <pybind11/pybind11.h>
diff --git a/limap/features/featuremap.h b/limap/features/featuremap.h
index 3d047cb6..f95728bd 100644
--- a/limap/features/featuremap.h
+++ b/limap/features/featuremap.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_FEATURES_FEATUREMAP_H_
-#define LIMAP_FEATURES_FEATUREMAP_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -95,5 +94,3 @@ bool FeatureInterpolator<DTYPE, CHANNELS>::CheckBounds(const T *xy) const {
 } // namespace features
 
 } // namespace limap
-
-#endif
diff --git a/limap/features/featurepatch.h b/limap/features/featurepatch.h
index dd7122c5..e37d7760 100644
--- a/limap/features/featurepatch.h
+++ b/limap/features/featurepatch.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_FEATURES_FEATUREPATCH_H_
-#define LIMAP_FEATURES_FEATUREPATCH_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -140,5 +139,3 @@ bool PatchInterpolator<DTYPE, CHANNELS>::CheckBounds(const T *xy) const {
 } // namespace features
 
 } // namespace limap
-
-#endif
diff --git a/limap/features/line_patch_extractor.h b/limap/features/line_patch_extractor.h
index e8d3a3cf..ff4b4ec8 100644
--- a/limap/features/line_patch_extractor.h
+++ b/limap/features/line_patch_extractor.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_FEATURES_LINE_PATCH_EXTRACTOR_H_
-#define LIMAP_FEATURES_LINE_PATCH_EXTRACTOR_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -64,5 +63,3 @@ template <typename DTYPE, int CHANNELS> class LinePatchExtractor {
 } // namespace features
 
 } // namespace limap
-
-#endif
diff --git a/limap/fitting/line3d_estimator.h b/limap/fitting/line3d_estimator.h
index c5057f91..1ca96354 100644
--- a/limap/fitting/line3d_estimator.h
+++ b/limap/fitting/line3d_estimator.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_FITTING_LINE_ESTIMATOR_H_
-#define LIMAP_FITTING_LINE_ESTIMATOR_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include "base/infinite_line.h"
@@ -55,5 +54,3 @@ Fit3DPoints(const Eigen::Matrix3Xd points,
 } // namespace fitting
 
 } // namespace limap
-
-#endif
diff --git a/limap/merging/aggregator.h b/limap/merging/aggregator.h
index 8ce49e2e..801028e9 100644
--- a/limap/merging/aggregator.h
+++ b/limap/merging/aggregator.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_MERGING_AGGREGATOR_H_
-#define LIMAP_MERGING_AGGREGATOR_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -32,5 +31,3 @@ class Aggregator {
 } // namespace merging
 
 } // namespace limap
-
-#endif
diff --git a/limap/merging/merging.h b/limap/merging/merging.h
index 8f101858..f473b7e6 100644
--- a/limap/merging/merging.h
+++ b/limap/merging/merging.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_MERGING_MERGING_H_
-#define LIMAP_MERGING_MERGING_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -48,5 +47,3 @@ RemergeLineTracks(const std::vector<LineTrack> &linetracks,
 } // namespace merging
 
 } // namespace limap
-
-#endif
diff --git a/limap/merging/merging_utils.h b/limap/merging/merging_utils.h
index 2e75fe22..653d1e44 100644
--- a/limap/merging/merging_utils.h
+++ b/limap/merging/merging_utils.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_MERGING_MERGING_UTILS_H_
-#define LIMAP_MERGING_MERGING_UTILS_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -49,5 +48,3 @@ void FilterTracksByOverlap(std::vector<LineTrack> &new_linetracks,
 } // namespace merging
 
 } // namespace limap
-
-#endif
diff --git a/limap/optimize/global_pl_association/cost_functions.h b/limap/optimize/global_pl_association/cost_functions.h
index 79817188..eddf5a40 100644
--- a/limap/optimize/global_pl_association/cost_functions.h
+++ b/limap/optimize/global_pl_association/cost_functions.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_PL_TRIANGULATION_COST_FUNCTIONS_H_
-#define LIMAP_PL_TRIANGULATION_COST_FUNCTIONS_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <ceres/ceres.h>
@@ -204,5 +203,3 @@ struct VPCollinearityFunctor {
 } // namespace optimize
 
 } // namespace limap
-
-#endif
diff --git a/limap/optimize/global_pl_association/global_associator.h b/limap/optimize/global_pl_association/global_associator.h
index a6c8848b..c242504f 100644
--- a/limap/optimize/global_pl_association/global_associator.h
+++ b/limap/optimize/global_pl_association/global_associator.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_OPTIMIZE_GLOBAL_PL_ASSOCIATION_GLOBAL_ASSOCIATOR_H_
-#define LIMAP_OPTIMIZE_GLOBAL_PL_ASSOCIATION_GLOBAL_ASSOCIATOR_H_
+#pragma once
 
 #include "base/image_collection.h"
 #include "base/infinite_line.h"
@@ -160,5 +159,3 @@ class GlobalAssociator : public hybrid_bundle_adjustment::HybridBAEngine {
 } // namespace optimize
 
 } // namespace limap
-
-#endif
diff --git a/limap/optimize/hybrid_bundle_adjustment/cost_functions.h b/limap/optimize/hybrid_bundle_adjustment/cost_functions.h
index 13183119..1e5e0dfb 100644
--- a/limap/optimize/hybrid_bundle_adjustment/cost_functions.h
+++ b/limap/optimize/hybrid_bundle_adjustment/cost_functions.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_HYBRIDBA_COST_FUNCTIONS_H_
-#define LIMAP_HYBRIDBA_COST_FUNCTIONS_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <ceres/ceres.h>
@@ -90,5 +89,3 @@ template <typename CameraModel> struct PointGeometricRefinementFunctor {
 } // namespace optimize
 
 } // namespace limap
-
-#endif
diff --git a/limap/optimize/hybrid_bundle_adjustment/hybrid_bundle_adjustment.h b/limap/optimize/hybrid_bundle_adjustment/hybrid_bundle_adjustment.h
index 66b32e23..6528aa16 100644
--- a/limap/optimize/hybrid_bundle_adjustment/hybrid_bundle_adjustment.h
+++ b/limap/optimize/hybrid_bundle_adjustment/hybrid_bundle_adjustment.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_OPTIMIZE_HYBRIDBA_HYBRIDBA_H_
-#define LIMAP_OPTIMIZE_HYBRIDBA_HYBRIDBA_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -74,5 +73,3 @@ class HybridBAEngine {
 } // namespace optimize
 
 } // namespace limap
-
-#endif
diff --git a/limap/optimize/hybrid_bundle_adjustment/hybrid_bundle_adjustment_config.h b/limap/optimize/hybrid_bundle_adjustment/hybrid_bundle_adjustment_config.h
index 40546057..629d890b 100644
--- a/limap/optimize/hybrid_bundle_adjustment/hybrid_bundle_adjustment_config.h
+++ b/limap/optimize/hybrid_bundle_adjustment/hybrid_bundle_adjustment_config.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_OPTIMIZE_HYBRIDBA_HYBRIDBA_CONFIG_H_
-#define LIMAP_OPTIMIZE_HYBRIDBA_HYBRIDBA_CONFIG_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -54,5 +53,3 @@ class HybridBAConfig : public line_refinement::RefinementConfig {
 } // namespace optimize
 
 } // namespace limap
-
-#endif
diff --git a/limap/optimize/hybrid_localization/cost_functions.h b/limap/optimize/hybrid_localization/cost_functions.h
index 930aad48..048069d2 100644
--- a/limap/optimize/hybrid_localization/cost_functions.h
+++ b/limap/optimize/hybrid_localization/cost_functions.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_OPTIMIZE_HYBRID_LOCALIZATION_COST_FUNCTIONS_H_
-#define LIMAP_OPTIMIZE_HYBRID_LOCALIZATION_COST_FUNCTIONS_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -352,5 +351,3 @@ struct ReprojectionPointFunctor {
 } // namespace optimize
 
 } // namespace limap
-
-#endif
diff --git a/limap/optimize/hybrid_localization/hybrid_localization.h b/limap/optimize/hybrid_localization/hybrid_localization.h
index 53e62f23..e847a9da 100644
--- a/limap/optimize/hybrid_localization/hybrid_localization.h
+++ b/limap/optimize/hybrid_localization/hybrid_localization.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_OPTIMIZE_HYBRID_LOCALIZATION_HYBRID_LOCALIZATION_H_
-#define LIMAP_OPTIMIZE_HYBRID_LOCALIZATION_HYBRID_LOCALIZATION_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -131,5 +130,3 @@ class JointLocEngine : public LineLocEngine {
 } // namespace optimize
 
 } // namespace limap
-
-#endif
diff --git a/limap/optimize/hybrid_localization/hybrid_localization_config.h b/limap/optimize/hybrid_localization/hybrid_localization_config.h
index 3118bdf3..9e246fea 100644
--- a/limap/optimize/hybrid_localization/hybrid_localization_config.h
+++ b/limap/optimize/hybrid_localization/hybrid_localization_config.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_OPTIMIZE_HYBRID_LOCALIZATION_HYBRID_LOCALIZATION_CONFIG_H_
-#define LIMAP_OPTIMIZE_HYBRID_LOCALIZATION_HYBRID_LOCALIZATION_CONFIG_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -83,5 +82,3 @@ class LineLocConfig {
 } // namespace optimize
 
 } // namespace limap
-
-#endif
diff --git a/limap/optimize/line_refinement/cost_functions.h b/limap/optimize/line_refinement/cost_functions.h
index 3d053915..c7dcf27a 100644
--- a/limap/optimize/line_refinement/cost_functions.h
+++ b/limap/optimize/line_refinement/cost_functions.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_OPTIMIZE_LINE_REFINEMENT_COST_FUNCTIONS_H_
-#define LIMAP_OPTIMIZE_LINE_REFINEMENT_COST_FUNCTIONS_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -199,5 +198,3 @@ template <typename CameraModel> struct GeometricRefinementFunctor {
 } // namespace optimize
 
 } // namespace limap
-
-#endif
diff --git a/limap/optimize/line_refinement/pixel_cost_functions.h b/limap/optimize/line_refinement/pixel_cost_functions.h
index bb927f61..8d6df5b3 100644
--- a/limap/optimize/line_refinement/pixel_cost_functions.h
+++ b/limap/optimize/line_refinement/pixel_cost_functions.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_OPTIMIZE_LINE_REFINEMENT_PIXEL_COST_FUNCTIONS_H_
-#define LIMAP_OPTIMIZE_LINE_REFINEMENT_PIXEL_COST_FUNCTIONS_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -405,5 +404,3 @@ struct FeatureConsisTgtFunctor {
 } // namespace optimize
 
 } // namespace limap
-
-#endif
diff --git a/limap/optimize/line_refinement/refine.h b/limap/optimize/line_refinement/refine.h
index 40815703..4479a006 100644
--- a/limap/optimize/line_refinement/refine.h
+++ b/limap/optimize/line_refinement/refine.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_OPTIMIZE_LINE_REFINEMENT_REFINE_H_
-#define LIMAP_OPTIMIZE_LINE_REFINEMENT_REFINE_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -149,5 +148,3 @@ template <typename DTYPE, int CHANNELS> class RefinementEngine {
 } // namespace optimize
 
 } // namespace limap
-
-#endif
diff --git a/limap/optimize/line_refinement/refinement_config.h b/limap/optimize/line_refinement/refinement_config.h
index 05f2d5c1..1fdfbec1 100644
--- a/limap/optimize/line_refinement/refinement_config.h
+++ b/limap/optimize/line_refinement/refinement_config.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_OPTIMIZE_LINE_REFINEMENT_REFINEMENT_CONFIG_H_
-#define LIMAP_OPTIMIZE_LINE_REFINEMENT_REFINEMENT_CONFIG_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -95,5 +94,3 @@ class RefinementConfig {
 } // namespace optimize
 
 } // namespace limap
-
-#endif
diff --git a/limap/pointsfm/sfm_model.h b/limap/pointsfm/sfm_model.h
index b8cc0517..027e798d 100644
--- a/limap/pointsfm/sfm_model.h
+++ b/limap/pointsfm/sfm_model.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_POINTSFM_SFM_MODEL_H_
-#define LIMAP_POINTSFM_SFM_MODEL_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -68,5 +67,3 @@ class SfmModel : public colmap::mvs::Model {
 } // namespace pointsfm
 
 } // namespace limap
-
-#endif
diff --git a/limap/solvers/triangulation/triangulate_line_with_one_point.h b/limap/solvers/triangulation/triangulate_line_with_one_point.h
index e1978b67..17877caa 100644
--- a/limap/solvers/triangulation/triangulate_line_with_one_point.h
+++ b/limap/solvers/triangulation/triangulate_line_with_one_point.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_SOLVERS_TRIANGULATION_TRIANGULATE_LINE_WITH_ONE_POINT_H_
-#define LIMAP_SOLVERS_TRIANGULATION_TRIANGULATE_LINE_WITH_ONE_POINT_H_
+#pragma once
 
 #include <Eigen/Core>
 
@@ -40,5 +39,3 @@ std::pair<double, double> triangulate_line_with_one_point(
 } // namespace solvers
 
 } // namespace limap
-
-#endif
diff --git a/limap/structures/pl_bipartite.h b/limap/structures/pl_bipartite.h
index 0b588a02..38ac5844 100644
--- a/limap/structures/pl_bipartite.h
+++ b/limap/structures/pl_bipartite.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_STRUCTURES_PL_BIPARTITE_H
-#define LIMAP_STRUCTURES_PL_BIPARTITE_H
+#pragma once
 
 #include <pybind11/eigen.h>
 #include <pybind11/numpy.h>
@@ -76,5 +75,3 @@ class PL_Bipartite3d : public PL_Bipartite<PointTrack, LineTrack> {
 } // namespace structures
 
 } // namespace limap
-
-#endif
diff --git a/limap/structures/pl_bipartite_base.h b/limap/structures/pl_bipartite_base.h
index 1aadbb98..12249960 100644
--- a/limap/structures/pl_bipartite_base.h
+++ b/limap/structures/pl_bipartite_base.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_STRUCTURES_PL_BIPARTITE_BASE_H
-#define LIMAP_STRUCTURES_PL_BIPARTITE_BASE_H
+#pragma once
 
 #include <map>
 #include <set>
@@ -100,5 +99,3 @@ template <typename PTYPE, typename LTYPE> class PL_Bipartite {
 } // namespace structures
 
 } // namespace limap
-
-#endif
diff --git a/limap/structures/vpline_bipartite.h b/limap/structures/vpline_bipartite.h
index fe424a7c..86bc9f3e 100644
--- a/limap/structures/vpline_bipartite.h
+++ b/limap/structures/vpline_bipartite.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_STRUCTURES_VPLINE_BIPARTITE_H
-#define LIMAP_STRUCTURES_VPLINE_BIPARTITE_H
+#pragma once
 
 #include <pybind11/eigen.h>
 #include <pybind11/numpy.h>
@@ -47,5 +46,3 @@ std::map<int, VPLine_Bipartite2d> GetAllBipartites_VPLine2d(
 } // namespace structures
 
 } // namespace limap
-
-#endif
diff --git a/limap/triangulation/base_line_triangulator.h b/limap/triangulation/base_line_triangulator.h
index 6704ed3f..8d444a23 100644
--- a/limap/triangulation/base_line_triangulator.h
+++ b/limap/triangulation/base_line_triangulator.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_TRIANGULATION_BASE_LINE_TRIANGULATOR_H_
-#define LIMAP_TRIANGULATION_BASE_LINE_TRIANGULATOR_H_
+#pragma once
 
 #include "base/image_collection.h"
 #include "base/linebase.h"
@@ -147,5 +146,3 @@ class BaseLineTriangulator {
 } // namespace triangulation
 
 } // namespace limap
-
-#endif
diff --git a/limap/triangulation/functions.h b/limap/triangulation/functions.h
index 1c12b6a2..88eadbb0 100644
--- a/limap/triangulation/functions.h
+++ b/limap/triangulation/functions.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_TRIANGULATION_FUNCTIONS_H_
-#define LIMAP_TRIANGULATION_FUNCTIONS_H_
+#pragma once
 
 #include "base/camera_view.h"
 #include "base/infinite_line.h"
@@ -80,5 +79,3 @@ Line3d triangulate_line_with_direction(const Line2d &l1,
 } // namespace triangulation
 
 } // namespace limap
-
-#endif
diff --git a/limap/triangulation/global_line_triangulator.h b/limap/triangulation/global_line_triangulator.h
index ef9e55b3..223bfe79 100644
--- a/limap/triangulation/global_line_triangulator.h
+++ b/limap/triangulation/global_line_triangulator.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_TRIANGULATION_GLOBAL_LINE_TRIANGULATOR_H_
-#define LIMAP_TRIANGULATION_GLOBAL_LINE_TRIANGULATOR_H_
+#pragma once
 
 #include "base/graph.h"
 #include "base/line_linker.h"
@@ -107,5 +106,3 @@ class GlobalLineTriangulator : public BaseLineTriangulator {
 } // namespace triangulation
 
 } // namespace limap
-
-#endif
diff --git a/limap/undistortion/undistort.h b/limap/undistortion/undistort.h
index d8cb517c..75b71dc7 100644
--- a/limap/undistortion/undistort.h
+++ b/limap/undistortion/undistort.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_UNDISTORTION_UNDISTORT_H_
-#define LIMAP_UNDISTORTION_UNDISTORT_H_
+#pragma once
 
 #include "_limap/helpers.h"
 #include <pybind11/numpy.h>
@@ -30,5 +29,3 @@ std::vector<V2D> UndistortPoints(const std::vector<V2D> &points,
 } // namespace undistortion
 
 } // namespace limap
-
-#endif
diff --git a/limap/util/simple_logger.h b/limap/util/simple_logger.h
index b4a1a575..f150b6ec 100644
--- a/limap/util/simple_logger.h
+++ b/limap/util/simple_logger.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_UTIL_SIMPLE_LOGGER_H_
-#define LIMAP_UTIL_SIMPLE_LOGGER_H_
+#pragma once
 
 /*
  * File:   Log.h
@@ -200,5 +199,3 @@ class SyncLogProgressbar : public progressbar {
 };
 
 } // namespace limap
-
-#endif
diff --git a/limap/vplib/JLinkage/JLinkage.h b/limap/vplib/JLinkage/JLinkage.h
index 59504814..b85ce872 100644
--- a/limap/vplib/JLinkage/JLinkage.h
+++ b/limap/vplib/JLinkage/JLinkage.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_VPLIB_JLINKAGE_H_
-#define LIMAP_VPLIB_JLINKAGE_H_
+#pragma once
 
 #include <pybind11/eigen.h>
 #include <pybind11/pybind11.h>
@@ -44,5 +43,3 @@ class JLinkage : public BaseVPDetector {
 } // namespace vplib
 
 } // namespace limap
-
-#endif
diff --git a/limap/vplib/base_vp_detector.h b/limap/vplib/base_vp_detector.h
index 8febcaee..fe4b7162 100644
--- a/limap/vplib/base_vp_detector.h
+++ b/limap/vplib/base_vp_detector.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_VPLIB_BASE_VP_DETECTOR_H_
-#define LIMAP_VPLIB_BASE_VP_DETECTOR_H_
+#pragma once
 
 #include <pybind11/eigen.h>
 #include <pybind11/pybind11.h>
@@ -53,5 +52,3 @@ class BaseVPDetector {
 } // namespace vplib
 
 } // namespace limap
-
-#endif
diff --git a/limap/vplib/global_vptrack_constructor.h b/limap/vplib/global_vptrack_constructor.h
index 2146527a..01422bcd 100644
--- a/limap/vplib/global_vptrack_constructor.h
+++ b/limap/vplib/global_vptrack_constructor.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_VPLIB_GLOBAL_VPTRACK_CONSTRUCTOR_H_
-#define LIMAP_VPLIB_GLOBAL_VPTRACK_CONSTRUCTOR_H_
+#pragma once
 
 #include <pybind11/eigen.h>
 #include <pybind11/pybind11.h>
@@ -54,5 +53,3 @@ class GlobalVPTrackConstructor {
 } // namespace vplib
 
 } // namespace limap
-
-#endif
diff --git a/limap/vplib/vpbase.h b/limap/vplib/vpbase.h
index 63c59b9d..19ed9b3b 100644
--- a/limap/vplib/vpbase.h
+++ b/limap/vplib/vpbase.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_VPLIB_VPBASE_H_
-#define LIMAP_VPLIB_VPBASE_H_
+#pragma once
 
 #include <pybind11/eigen.h>
 #include <pybind11/pybind11.h>
@@ -48,5 +47,3 @@ class VPResult {
 } // namespace vplib
 
 } // namespace limap
-
-#endif
diff --git a/limap/vplib/vptrack.h b/limap/vplib/vptrack.h
index 47aac55b..3b40d532 100644
--- a/limap/vplib/vptrack.h
+++ b/limap/vplib/vptrack.h
@@ -1,5 +1,4 @@
-#ifndef LIMAP_VPLIB_VPTRACK_H_
-#define LIMAP_VPLIB_VPTRACK_H_
+#pragma once
 
 #include <pybind11/eigen.h>
 #include <pybind11/pybind11.h>
@@ -41,5 +40,3 @@ MergeVPTracksByDirection(const std::vector<VPTrack> &vptracks,
 } // namespace vplib
 
 } // namespace limap
-
-#endif