Remove wrong use of PROJ_CONST_DECL and replace it with PROJ_PURE_DECL (fixes #1224)

8 files changed, 101 insertions, 101 deletions

diff --git a/src/iso19111/common.cpp b/src/iso19111/common.cpp
index 32eb5a95..ad816275 100644
--- a/src/iso19111/common.cpp
+++ b/src/iso19111/common.cpp
@@ -329,7 +329,7 @@ Measure::~Measure() = default;
 
 /** \brief Return the unit of the Measure.
  */
-const UnitOfMeasure &Measure::unit() PROJ_CONST_DEFN { return d->unit_; }
+const UnitOfMeasure &Measure::unit() PROJ_PURE_DEFN { return d->unit_; }
 
 // ---------------------------------------------------------------------------
 
@@ -337,7 +337,7 @@ const UnitOfMeasure &Measure::unit() PROJ_CONST_DEFN { return d->unit_; }
  * corresponding
  * unit of the International System.
  */
-double Measure::getSIValue() PROJ_CONST_DEFN {
+double Measure::getSIValue() PROJ_PURE_DEFN {
     return d->value_ * d->unit_.conversionToSI();
 }
 
@@ -345,13 +345,13 @@ double Measure::getSIValue() PROJ_CONST_DEFN {
 
 /** \brief Return the value of the measure, expressed in the unit()
  */
-double Measure::value() PROJ_CONST_DEFN { return d->value_; }
+double Measure::value() PROJ_PURE_DEFN { return d->value_; }
 
 // ---------------------------------------------------------------------------
 
 /** \brief Return the value of this measure expressed into the provided unit.
  */
-double Measure::convertToUnit(const UnitOfMeasure &otherUnit) PROJ_CONST_DEFN {
+double Measure::convertToUnit(const UnitOfMeasure &otherUnit) PROJ_PURE_DEFN {
     return getSIValue() / otherUnit.conversionToSI();
 }
 
@@ -361,7 +361,7 @@ double Measure::convertToUnit(const UnitOfMeasure &otherUnit) PROJ_CONST_DEFN {
  *
  * The comparison is done both on the value and the unit.
  */
-bool Measure::operator==(const Measure &other) PROJ_CONST_DEFN {
+bool Measure::operator==(const Measure &other) PROJ_PURE_DEFN {
     return d->value_ == other.d->value_ && d->unit_ == other.d->unit_;
 }
 
@@ -567,7 +567,7 @@ IdentifiedObject::~IdentifiedObject() = default;
  * Generally, the only interesting field of the name will be
  * name()->description().
  */
-const IdentifierNNPtr &IdentifiedObject::name() PROJ_CONST_DEFN {
+const IdentifierNNPtr &IdentifiedObject::name() PROJ_PURE_DEFN {
     return d->name;
 }
 
@@ -577,7 +577,7 @@ const IdentifierNNPtr &IdentifiedObject::name() PROJ_CONST_DEFN {
  *
  * Return *(name()->description())
  */
-const std::string &IdentifiedObject::nameStr() PROJ_CONST_DEFN {
+const std::string &IdentifiedObject::nameStr() PROJ_PURE_DEFN {
     return *(d->name->description());
 }
 
@@ -589,7 +589,7 @@ const std::string &IdentifiedObject::nameStr() PROJ_CONST_DEFN {
  * filled.
  */
 const std::vector<IdentifierNNPtr> &
-IdentifiedObject::identifiers() PROJ_CONST_DEFN {
+IdentifiedObject::identifiers() PROJ_PURE_DEFN {
     return d->identifiers;
 }
 
@@ -598,7 +598,7 @@ IdentifiedObject::identifiers() PROJ_CONST_DEFN {
 /** \brief Return the alias(es) of the object.
  */
 const std::vector<GenericNameNNPtr> &
-IdentifiedObject::aliases() PROJ_CONST_DEFN {
+IdentifiedObject::aliases() PROJ_PURE_DEFN {
     return d->aliases;
 }
 
@@ -608,7 +608,7 @@ IdentifiedObject::aliases() PROJ_CONST_DEFN {
  *
  * Shortcut for aliases()[0]->toFullyQualifiedName()->toString()
  */
-std::string IdentifiedObject::alias() PROJ_CONST_DEFN {
+std::string IdentifiedObject::alias() PROJ_PURE_DEFN {
     if (d->aliases.empty())
         return std::string();
     return d->aliases[0]->toFullyQualifiedName()->toString();
@@ -619,7 +619,7 @@ std::string IdentifiedObject::alias() PROJ_CONST_DEFN {
 /** \brief Return the EPSG code.
  * @return code, or 0 if not found
  */
-int IdentifiedObject::getEPSGCode() PROJ_CONST_DEFN {
+int IdentifiedObject::getEPSGCode() PROJ_PURE_DEFN {
     for (const auto &id : identifiers()) {
         if (ci_equal(*(id->codeSpace()), metadata::Identifier::EPSG)) {
             return ::atoi(id->code().c_str());
@@ -632,7 +632,7 @@ int IdentifiedObject::getEPSGCode() PROJ_CONST_DEFN {
 
 /** \brief Return the remarks.
  */
-const std::string &IdentifiedObject::remarks() PROJ_CONST_DEFN {
+const std::string &IdentifiedObject::remarks() PROJ_PURE_DEFN {
     return d->remarks;
 }
 
@@ -642,7 +642,7 @@ const std::string &IdentifiedObject::remarks() PROJ_CONST_DEFN {
  *
  * \remark Extension of \ref ISO_19111_2018
  */
-bool IdentifiedObject::isDeprecated() PROJ_CONST_DEFN {
+bool IdentifiedObject::isDeprecated() PROJ_PURE_DEFN {
     return d->isDeprecated;
 }
 
@@ -829,7 +829,7 @@ bool IdentifiedObject::_isEquivalentTo(
 
 bool IdentifiedObject::_isEquivalentTo(const IdentifiedObject *otherIdObj,
                                        util::IComparable::Criterion criterion)
-    PROJ_CONST_DEFN {
+    PROJ_PURE_DEFN {
     if (criterion == util::IComparable::Criterion::STRICT) {
         if (!ci_equal(nameStr(), otherIdObj->nameStr())) {
             return false;
@@ -883,7 +883,7 @@ ObjectDomain::~ObjectDomain() = default;
  *
  * @return the scope, or empty.
  */
-const optional<std::string> &ObjectDomain::scope() PROJ_CONST_DEFN {
+const optional<std::string> &ObjectDomain::scope() PROJ_PURE_DEFN {
     return d->scope_;
 }
 
@@ -893,7 +893,7 @@ const optional<std::string> &ObjectDomain::scope() PROJ_CONST_DEFN {
  *
  * @return the domain of validity, or nullptr.
  */
-const ExtentPtr &ObjectDomain::domainOfValidity() PROJ_CONST_DEFN {
+const ExtentPtr &ObjectDomain::domainOfValidity() PROJ_PURE_DEFN {
     return d->domainOfValidity_;
 }
 
@@ -1013,7 +1013,7 @@ ObjectUsage::~ObjectUsage() = default;
 
 /** \brief Return the domains of the object.
  */
-const std::vector<ObjectDomainNNPtr> &ObjectUsage::domains() PROJ_CONST_DEFN {
+const std::vector<ObjectDomainNNPtr> &ObjectUsage::domains() PROJ_PURE_DEFN {
     return d->domains_;
 }
 
diff --git a/src/iso19111/coordinateoperation.cpp b/src/iso19111/coordinateoperation.cpp
index cc212314..4bdb3779 100644
--- a/src/iso19111/coordinateoperation.cpp
+++ b/src/iso19111/coordinateoperation.cpp
@@ -760,7 +760,7 @@ OperationMethod::~OperationMethod() = default;
  *
  * @return the formula, or empty.
  */
-const util::optional<std::string> &OperationMethod::formula() PROJ_CONST_DEFN {
+const util::optional<std::string> &OperationMethod::formula() PROJ_PURE_DEFN {
     return d->formula_;
 }
 
@@ -773,7 +773,7 @@ const util::optional<std::string> &OperationMethod::formula() PROJ_CONST_DEFN {
  * @return the formula citation, or empty.
  */
 const util::optional<metadata::Citation> &
-OperationMethod::formulaCitation() PROJ_CONST_DEFN {
+OperationMethod::formulaCitation() PROJ_PURE_DEFN {
     return d->formulaCitation_;
 }
 
@@ -784,7 +784,7 @@ OperationMethod::formulaCitation() PROJ_CONST_DEFN {
  * @return the parameters.
  */
 const std::vector<GeneralOperationParameterNNPtr> &
-OperationMethod::parameters() PROJ_CONST_DEFN {
+OperationMethod::parameters() PROJ_PURE_DEFN {
     return d->parameters_;
 }
 
@@ -835,7 +835,7 @@ OperationMethodNNPtr OperationMethod::create(
 /** \brief Return the EPSG code, either directly, or through the name
  * @return code, or 0 if not found
  */
-int OperationMethod::getEPSGCode() PROJ_CONST_DEFN {
+int OperationMethod::getEPSGCode() PROJ_PURE_DEFN {
     int epsg_code = IdentifiedObject::getEPSGCode();
     if (epsg_code == 0) {
         auto l_name = nameStr();
@@ -1012,7 +1012,7 @@ OperationParameterValue::~OperationParameterValue() = default;
  * @return the parameter (definition).
  */
 const OperationParameterNNPtr &
-OperationParameterValue::parameter() PROJ_CONST_DEFN {
+OperationParameterValue::parameter() PROJ_PURE_DEFN {
     return d->parameter;
 }
 
@@ -1023,7 +1023,7 @@ OperationParameterValue::parameter() PROJ_CONST_DEFN {
  * @return the parameter value.
  */
 const ParameterValueNNPtr &
-OperationParameterValue::parameterValue() PROJ_CONST_DEFN {
+OperationParameterValue::parameterValue() PROJ_PURE_DEFN {
     return d->parameterValue;
 }
 
@@ -1273,7 +1273,7 @@ const char *OperationParameter::getNameForEPSGCode(int epsg_code) noexcept {
 /** \brief Return the EPSG code, either directly, or through the name
  * @return code, or 0 if not found
  */
-int OperationParameter::getEPSGCode() PROJ_CONST_DEFN {
+int OperationParameter::getEPSGCode() PROJ_PURE_DEFN {
     int epsg_code = IdentifiedObject::getEPSGCode();
     if (epsg_code == 0) {
         const auto &l_name = nameStr();
@@ -1335,7 +1335,7 @@ SingleOperation::~SingleOperation() = default;
  * @return the parameter values.
  */
 const std::vector<GeneralParameterValueNNPtr> &
-SingleOperation::parameterValues() PROJ_CONST_DEFN {
+SingleOperation::parameterValues() PROJ_PURE_DEFN {
     return d->parameterValues_;
 }
 
@@ -1345,7 +1345,7 @@ SingleOperation::parameterValues() PROJ_CONST_DEFN {
  *
  * @return the operation method.
  */
-const OperationMethodNNPtr &SingleOperation::method() PROJ_CONST_DEFN {
+const OperationMethodNNPtr &SingleOperation::method() PROJ_PURE_DEFN {
     return d->method_;
 }
 
@@ -2072,7 +2072,7 @@ ParameterValueNNPtr ParameterValue::create(bool booleanValueIn) {
  *
  * @return the type.
  */
-const ParameterValue::Type &ParameterValue::type() PROJ_CONST_DEFN {
+const ParameterValue::Type &ParameterValue::type() PROJ_PURE_DEFN {
     return d->type_;
 }
 
@@ -2081,7 +2081,7 @@ const ParameterValue::Type &ParameterValue::type() PROJ_CONST_DEFN {
 /** \brief Returns the value as a Measure (assumes type() == Type::MEASURE)
  * @return the value as a Measure.
  */
-const common::Measure &ParameterValue::value() PROJ_CONST_DEFN {
+const common::Measure &ParameterValue::value() PROJ_PURE_DEFN {
     return *d->measure_;
 }
 
@@ -2090,7 +2090,7 @@ const common::Measure &ParameterValue::value() PROJ_CONST_DEFN {
 /** \brief Returns the value as a string (assumes type() == Type::STRING)
  * @return the value as a string.
  */
-const std::string &ParameterValue::stringValue() PROJ_CONST_DEFN {
+const std::string &ParameterValue::stringValue() PROJ_PURE_DEFN {
     return *d->stringValue_;
 }
 
@@ -2099,7 +2099,7 @@ const std::string &ParameterValue::stringValue() PROJ_CONST_DEFN {
 /** \brief Returns the value as a filename (assumes type() == Type::FILENAME)
  * @return the value as a filename.
  */
-const std::string &ParameterValue::valueFile() PROJ_CONST_DEFN {
+const std::string &ParameterValue::valueFile() PROJ_PURE_DEFN {
     return *d->stringValue_;
 }
 
@@ -2108,14 +2108,14 @@ const std::string &ParameterValue::valueFile() PROJ_CONST_DEFN {
 /** \brief Returns the value as a integer (assumes type() == Type::INTEGER)
  * @return the value as a integer.
  */
-int ParameterValue::integerValue() PROJ_CONST_DEFN { return d->integerValue_; }
+int ParameterValue::integerValue() PROJ_PURE_DEFN { return d->integerValue_; }
 
 // ---------------------------------------------------------------------------
 
 /** \brief Returns the value as a boolean (assumes type() == Type::BOOLEAN)
  * @return the value as a boolean.
  */
-bool ParameterValue::booleanValue() PROJ_CONST_DEFN { return d->booleanValue_; }
+bool ParameterValue::booleanValue() PROJ_PURE_DEFN { return d->booleanValue_; }
 
 // ---------------------------------------------------------------------------
 
@@ -5917,7 +5917,7 @@ Transformation::Transformation(const Transformation &other)
  *
  * @return the source CRS.
  */
-const crs::CRSNNPtr &Transformation::sourceCRS() PROJ_CONST_DEFN {
+const crs::CRSNNPtr &Transformation::sourceCRS() PROJ_PURE_DEFN {
     return CoordinateOperation::getPrivate()->strongRef_->sourceCRS_;
 }
 
@@ -5927,7 +5927,7 @@ const crs::CRSNNPtr &Transformation::sourceCRS() PROJ_CONST_DEFN {
  *
  * @return the target CRS.
  */
-const crs::CRSNNPtr &Transformation::targetCRS() PROJ_CONST_DEFN {
+const crs::CRSNNPtr &Transformation::targetCRS() PROJ_PURE_DEFN {
     return CoordinateOperation::getPrivate()->strongRef_->targetCRS_;
 }
 
diff --git a/src/iso19111/coordinatesystem.cpp b/src/iso19111/coordinatesystem.cpp
index 2281c78f..efba8c05 100644
--- a/src/iso19111/coordinatesystem.cpp
+++ b/src/iso19111/coordinatesystem.cpp
@@ -104,7 +104,7 @@ Meridian::~Meridian() = default;
  *
  * @return the longitude.
  */
-const common::Angle &Meridian::longitude() PROJ_CONST_DEFN {
+const common::Angle &Meridian::longitude() PROJ_PURE_DEFN {
     return d->longitude_;
 }
 
@@ -178,7 +178,7 @@ CoordinateSystemAxis::~CoordinateSystemAxis() = default;
  *
  * @return the abbreviation.
  */
-const std::string &CoordinateSystemAxis::abbreviation() PROJ_CONST_DEFN {
+const std::string &CoordinateSystemAxis::abbreviation() PROJ_PURE_DEFN {
     return d->abbreviation;
 }
 
@@ -199,7 +199,7 @@ const std::string &CoordinateSystemAxis::abbreviation() PROJ_CONST_DEFN {
  *
  * @return the direction.
  */
-const AxisDirection &CoordinateSystemAxis::direction() PROJ_CONST_DEFN {
+const AxisDirection &CoordinateSystemAxis::direction() PROJ_PURE_DEFN {
     return *(d->direction);
 }
 
@@ -213,7 +213,7 @@ const AxisDirection &CoordinateSystemAxis::direction() PROJ_CONST_DEFN {
  *
  * @return the axis unit.
  */
-const common::UnitOfMeasure &CoordinateSystemAxis::unit() PROJ_CONST_DEFN {
+const common::UnitOfMeasure &CoordinateSystemAxis::unit() PROJ_PURE_DEFN {
     return d->unit;
 }
 
@@ -225,7 +225,7 @@ const common::UnitOfMeasure &CoordinateSystemAxis::unit() PROJ_CONST_DEFN {
  * @return the minimum value, or empty.
  */
 const util::optional<double> &
-CoordinateSystemAxis::minimumValue() PROJ_CONST_DEFN {
+CoordinateSystemAxis::minimumValue() PROJ_PURE_DEFN {
     return d->minimumValue;
 }
 
@@ -237,7 +237,7 @@ CoordinateSystemAxis::minimumValue() PROJ_CONST_DEFN {
  * @return the maximum value, or empty.
  */
 const util::optional<double> &
-CoordinateSystemAxis::maximumValue() PROJ_CONST_DEFN {
+CoordinateSystemAxis::maximumValue() PROJ_PURE_DEFN {
     return d->maximumValue;
 }
 
@@ -249,7 +249,7 @@ CoordinateSystemAxis::maximumValue() PROJ_CONST_DEFN {
  *
  * @return the meridian, or null.
  */
-const MeridianPtr &CoordinateSystemAxis::meridian() PROJ_CONST_DEFN {
+const MeridianPtr &CoordinateSystemAxis::meridian() PROJ_PURE_DEFN {
     return d->meridian;
 }
 
@@ -466,7 +466,7 @@ CoordinateSystem::~CoordinateSystem() = default;
  * @return the axes.
  */
 const std::vector<CoordinateSystemAxisNNPtr> &
-CoordinateSystem::axisList() PROJ_CONST_DEFN {
+CoordinateSystem::axisList() PROJ_PURE_DEFN {
     return d->axisList;
 }
 
diff --git a/src/iso19111/crs.cpp b/src/iso19111/crs.cpp
index 81c70b7a..7bf40e0c 100644
--- a/src/iso19111/crs.cpp
+++ b/src/iso19111/crs.cpp
@@ -131,7 +131,7 @@ CRS::~CRS() = default;
  *
  * @return a BoundCRSPtr, that might be null.
  */
-const BoundCRSPtr &CRS::canonicalBoundCRS() PROJ_CONST_DEFN {
+const BoundCRSPtr &CRS::canonicalBoundCRS() PROJ_PURE_DEFN {
     return d->canonicalBoundCRS_;
 }
 
@@ -712,7 +712,7 @@ SingleCRS::~SingleCRS() = default;
  *
  * @return a Datum that might be null.
  */
-const datum::DatumPtr &SingleCRS::datum() PROJ_CONST_DEFN { return d->datum; }
+const datum::DatumPtr &SingleCRS::datum() PROJ_PURE_DEFN { return d->datum; }
 
 // ---------------------------------------------------------------------------
 
@@ -722,7 +722,7 @@ const datum::DatumPtr &SingleCRS::datum() PROJ_CONST_DEFN { return d->datum; }
  *
  * @return a DatumEnsemble that might be null.
  */
-const datum::DatumEnsemblePtr &SingleCRS::datumEnsemble() PROJ_CONST_DEFN {
+const datum::DatumEnsemblePtr &SingleCRS::datumEnsemble() PROJ_PURE_DEFN {
     return d->datumEnsemble;
 }
 
@@ -734,7 +734,7 @@ const datum::DatumEnsemblePtr &SingleCRS::datumEnsemble() PROJ_CONST_DEFN {
  *
  * @return a CoordinateSystem that might be null.
  */
-const cs::CoordinateSystemNNPtr &SingleCRS::coordinateSystem() PROJ_CONST_DEFN {
+const cs::CoordinateSystemNNPtr &SingleCRS::coordinateSystem() PROJ_PURE_DEFN {
     return d->coordinateSystem;
 }
 
@@ -873,7 +873,7 @@ CRSNNPtr GeodeticCRS::_shallowClone() const {
  * @return a GeodeticReferenceFrame or null (in which case datumEnsemble()
  * should return a non-null pointer.)
  */
-const datum::GeodeticReferenceFramePtr &GeodeticCRS::datum() PROJ_CONST_DEFN {
+const datum::GeodeticReferenceFramePtr &GeodeticCRS::datum() PROJ_PURE_DEFN {
     return d->datum_;
 }
 
@@ -895,7 +895,7 @@ static datum::GeodeticReferenceFrame *oneDatum(const GeodeticCRS *crs) {
  *
  * @return the PrimeMeridian.
  */
-const datum::PrimeMeridianNNPtr &GeodeticCRS::primeMeridian() PROJ_CONST_DEFN {
+const datum::PrimeMeridianNNPtr &GeodeticCRS::primeMeridian() PROJ_PURE_DEFN {
     if (d->datum_) {
         return d->datum_->primeMeridian();
     }
@@ -909,7 +909,7 @@ const datum::PrimeMeridianNNPtr &GeodeticCRS::primeMeridian() PROJ_CONST_DEFN {
  *
  * @return the PrimeMeridian.
  */
-const datum::EllipsoidNNPtr &GeodeticCRS::ellipsoid() PROJ_CONST_DEFN {
+const datum::EllipsoidNNPtr &GeodeticCRS::ellipsoid() PROJ_PURE_DEFN {
     if (d->datum_) {
         return d->datum_->ellipsoid();
     }
@@ -923,7 +923,7 @@ const datum::EllipsoidNNPtr &GeodeticCRS::ellipsoid() PROJ_CONST_DEFN {
  * @return a velocity model. might be null.
  */
 const std::vector<operation::PointMotionOperationNNPtr> &
-GeodeticCRS::velocityModel() PROJ_CONST_DEFN {
+GeodeticCRS::velocityModel() PROJ_PURE_DEFN {
     return d->velocityModel;
 }
 
@@ -938,7 +938,7 @@ GeodeticCRS::velocityModel() PROJ_CONST_DEFN {
  *
  * @return true if the CRS is a geocentric CRS.
  */
-bool GeodeticCRS::isGeocentric() PROJ_CONST_DEFN {
+bool GeodeticCRS::isGeocentric() PROJ_PURE_DEFN {
     const auto &cs = coordinateSystem();
     const auto &axisList = cs->axisList();
     return axisList.size() == 3 &&
@@ -1653,7 +1653,7 @@ CRSNNPtr GeographicCRS::_shallowClone() const {
  * @return a EllipsoidalCS.
  */
 const cs::EllipsoidalCSNNPtr &
-GeographicCRS::coordinateSystem() PROJ_CONST_DEFN {
+GeographicCRS::coordinateSystem() PROJ_PURE_DEFN {
     return d->coordinateSystem_;
 }
 
@@ -1715,7 +1715,7 @@ GeographicCRS::create(const util::PropertyMap &properties,
  * other 3D GeographicCRS.
  */
 bool GeographicCRS::is2DPartOf3D(util::nn<const GeographicCRS *> other)
-    PROJ_CONST_DEFN {
+    PROJ_PURE_DEFN {
     const auto &axis = d->coordinateSystem_->axisList();
     const auto &otherAxis = other->d->coordinateSystem_->axisList();
     if (!(axis.size() == 2 && otherAxis.size() == 3)) {
@@ -2009,7 +2009,7 @@ const datum::VerticalReferenceFramePtr VerticalCRS::datum() const {
  * @return a geoid model. might be null
  */
 const std::vector<operation::TransformationNNPtr> &
-VerticalCRS::geoidModel() PROJ_CONST_DEFN {
+VerticalCRS::geoidModel() PROJ_PURE_DEFN {
     return d->geoidModel;
 }
 
@@ -2020,7 +2020,7 @@ VerticalCRS::geoidModel() PROJ_CONST_DEFN {
  * @return a velocity model. might be null.
  */
 const std::vector<operation::PointMotionOperationNNPtr> &
-VerticalCRS::velocityModel() PROJ_CONST_DEFN {
+VerticalCRS::velocityModel() PROJ_PURE_DEFN {
     return d->velocityModel;
 }
 
@@ -2370,7 +2370,7 @@ DerivedCRS::~DerivedCRS() = default;
  *
  * @return the base CRS.
  */
-const SingleCRSNNPtr &DerivedCRS::baseCRS() PROJ_CONST_DEFN {
+const SingleCRSNNPtr &DerivedCRS::baseCRS() PROJ_PURE_DEFN {
     return d->baseCRS_;
 }
 
@@ -2388,7 +2388,7 @@ const operation::ConversionNNPtr DerivedCRS::derivingConversion() const {
 
 //! @cond Doxygen_Suppress
 const operation::ConversionNNPtr &
-DerivedCRS::derivingConversionRef() PROJ_CONST_DEFN {
+DerivedCRS::derivingConversionRef() PROJ_PURE_DEFN {
     return d->derivingConversion_;
 }
 //! @endcond
@@ -2495,7 +2495,7 @@ CRSNNPtr ProjectedCRS::_shallowClone() const {
  *
  * @return the base CRS.
  */
-const GeodeticCRSNNPtr &ProjectedCRS::baseCRS() PROJ_CONST_DEFN {
+const GeodeticCRSNNPtr &ProjectedCRS::baseCRS() PROJ_PURE_DEFN {
     return d->baseCRS();
 }
 
@@ -2505,7 +2505,7 @@ const GeodeticCRSNNPtr &ProjectedCRS::baseCRS() PROJ_CONST_DEFN {
  *
  * @return a CartesianCS
  */
-const cs::CartesianCSNNPtr &ProjectedCRS::coordinateSystem() PROJ_CONST_DEFN {
+const cs::CartesianCSNNPtr &ProjectedCRS::coordinateSystem() PROJ_PURE_DEFN {
     return d->coordinateSystem();
 }
 
@@ -3209,7 +3209,7 @@ CRSNNPtr CompoundCRS::_shallowClone() const {
  * @return the components.
  */
 const std::vector<CRSNNPtr> &
-CompoundCRS::componentReferenceSystems() PROJ_CONST_DEFN {
+CompoundCRS::componentReferenceSystems() PROJ_PURE_DEFN {
     return d->components_;
 }
 
@@ -3552,7 +3552,7 @@ CRSNNPtr BoundCRS::_shallowClone() const { return shallowCloneAsBoundCRS(); }
  *
  * @return the base CRS.
  */
-const CRSNNPtr &BoundCRS::baseCRS() PROJ_CONST_DEFN { return d->baseCRS_; }
+const CRSNNPtr &BoundCRS::baseCRS() PROJ_PURE_DEFN { return d->baseCRS_; }
 
 // ---------------------------------------------------------------------------
 
@@ -3589,7 +3589,7 @@ CRSNNPtr BoundCRS::baseCRSWithCanonicalBoundCRS() const {
  *
  * @return the hub CRS.
  */
-const CRSNNPtr &BoundCRS::hubCRS() PROJ_CONST_DEFN { return d->hubCRS_; }
+const CRSNNPtr &BoundCRS::hubCRS() PROJ_PURE_DEFN { return d->hubCRS_; }
 
 // ---------------------------------------------------------------------------
 
@@ -3598,7 +3598,7 @@ const CRSNNPtr &BoundCRS::hubCRS() PROJ_CONST_DEFN { return d->hubCRS_; }
  * @return transformation.
  */
 const operation::TransformationNNPtr &
-BoundCRS::transformation() PROJ_CONST_DEFN {
+BoundCRS::transformation() PROJ_PURE_DEFN {
     return d->transformation_;
 }
 
diff --git a/src/iso19111/datum.cpp b/src/iso19111/datum.cpp
index cbd15195..9d2e48a0 100644
--- a/src/iso19111/datum.cpp
+++ b/src/iso19111/datum.cpp
@@ -267,7 +267,7 @@ PrimeMeridian::~PrimeMeridian() = default;
  *
  * @return the longitude of the prime meridian.
  */
-const common::Angle &PrimeMeridian::longitude() PROJ_CONST_DEFN {
+const common::Angle &PrimeMeridian::longitude() PROJ_PURE_DEFN {
     return d->longitude_;
 }
 
@@ -477,7 +477,7 @@ Ellipsoid::Ellipsoid(const Ellipsoid &other)
  *
  * @return the semi-major axis.
  */
-const common::Length &Ellipsoid::semiMajorAxis() PROJ_CONST_DEFN {
+const common::Length &Ellipsoid::semiMajorAxis() PROJ_PURE_DEFN {
     return d->semiMajorAxis_;
 }
 
@@ -494,7 +494,7 @@ const common::Length &Ellipsoid::semiMajorAxis() PROJ_CONST_DEFN {
  * @return the inverse flattening value of the ellipsoid, or empty.
  */
 const util::optional<common::Scale> &
-Ellipsoid::inverseFlattening() PROJ_CONST_DEFN {
+Ellipsoid::inverseFlattening() PROJ_PURE_DEFN {
     return d->inverseFlattening_;
 }
 
@@ -511,7 +511,7 @@ Ellipsoid::inverseFlattening() PROJ_CONST_DEFN {
  * @return the semi-minor axis of the ellipsoid, or empty.
  */
 const util::optional<common::Length> &
-Ellipsoid::semiMinorAxis() PROJ_CONST_DEFN {
+Ellipsoid::semiMinorAxis() PROJ_PURE_DEFN {
     return d->semiMinorAxis_;
 }
 
@@ -526,7 +526,7 @@ Ellipsoid::semiMinorAxis() PROJ_CONST_DEFN {
  *
  * @return true if the ellipsoid is spherical.
  */
-bool Ellipsoid::isSphere() PROJ_CONST_DEFN {
+bool Ellipsoid::isSphere() PROJ_PURE_DEFN {
     if (d->inverseFlattening_.has_value()) {
         return d->inverseFlattening_->value() == 0;
     }
@@ -547,7 +547,7 @@ bool Ellipsoid::isSphere() PROJ_CONST_DEFN {
  * @return the semi-median axis of the ellipsoid, or empty.
  */
 const util::optional<common::Length> &
-Ellipsoid::semiMedianAxis() PROJ_CONST_DEFN {
+Ellipsoid::semiMedianAxis() PROJ_PURE_DEFN {
     return d->semiMedianAxis_;
 }
 
@@ -560,7 +560,7 @@ Ellipsoid::semiMedianAxis() PROJ_CONST_DEFN {
  *
  * @return the inverse flattening value of the ellipsoid, or 0 for a sphere.
  */
-double Ellipsoid::computedInverseFlattening() PROJ_CONST_DEFN {
+double Ellipsoid::computedInverseFlattening() PROJ_PURE_DEFN {
     if (d->inverseFlattening_.has_value()) {
         return d->inverseFlattening_->getSIValue();
     }
@@ -580,7 +580,7 @@ double Ellipsoid::computedInverseFlattening() PROJ_CONST_DEFN {
  *
  * @return the squared eccentricity, or a negative value if invalid.
  */
-double Ellipsoid::squaredEccentricity() PROJ_CONST_DEFN {
+double Ellipsoid::squaredEccentricity() PROJ_PURE_DEFN {
     const double rf = computedInverseFlattening();
     const double f = rf != 0.0 ? 1. / rf : 0.0;
     const double e2 = f * (2 - f);
@@ -616,7 +616,7 @@ common::Length Ellipsoid::computeSemiMinorAxis() const {
 /** \brief Return the name of the celestial body on which the ellipsoid refers
  * to.
  */
-const std::string &Ellipsoid::celestialBody() PROJ_CONST_DEFN {
+const std::string &Ellipsoid::celestialBody() PROJ_PURE_DEFN {
     return d->celestialBody_;
 }
 
@@ -998,7 +998,7 @@ GeodeticReferenceFrame::~GeodeticReferenceFrame() = default;
  * @return the PrimeMeridian.
  */
 const PrimeMeridianNNPtr &
-GeodeticReferenceFrame::primeMeridian() PROJ_CONST_DEFN {
+GeodeticReferenceFrame::primeMeridian() PROJ_PURE_DEFN {
     return d->primeMeridian_;
 }
 
@@ -1014,7 +1014,7 @@ GeodeticReferenceFrame::primeMeridian() PROJ_CONST_DEFN {
  *
  * @return the Ellipsoid.
  */
-const EllipsoidNNPtr &GeodeticReferenceFrame::ellipsoid() PROJ_CONST_DEFN {
+const EllipsoidNNPtr &GeodeticReferenceFrame::ellipsoid() PROJ_PURE_DEFN {
     return d->ellipsoid_;
 }
 // ---------------------------------------------------------------------------
diff --git a/src/iso19111/factory.cpp b/src/iso19111/factory.cpp
index 3b6563f3..fb75b7b6 100644
--- a/src/iso19111/factory.cpp
+++ b/src/iso19111/factory.cpp
@@ -1170,9 +1170,9 @@ struct AuthorityFactory::Private {
             const std::string &authorityName)
         : context_(contextIn), authority_(authorityName) {}
 
-    inline const std::string &authority() PROJ_CONST_DEFN { return authority_; }
+    inline const std::string &authority() PROJ_PURE_DEFN { return authority_; }
 
-    inline const DatabaseContextNNPtr &context() PROJ_CONST_DEFN {
+    inline const DatabaseContextNNPtr &context() PROJ_PURE_DEFN {
         return context_;
     }
 
@@ -3732,7 +3732,7 @@ AuthorityFactory::createFromCRSCodesWithIntermediates(
 /** \brief Returns the authority name associated to this factory.
  * @return name.
  */
-const std::string &AuthorityFactory::getAuthority() PROJ_CONST_DEFN {
+const std::string &AuthorityFactory::getAuthority() PROJ_PURE_DEFN {
     return d->authority();
 }
 
diff --git a/src/iso19111/io.cpp b/src/iso19111/io.cpp
index 5be02ffe..52c484b9 100644
--- a/src/iso19111/io.cpp
+++ b/src/iso19111/io.cpp
@@ -782,15 +782,15 @@ struct WKTNode::Private {
     explicit Private(const std::string &valueIn) : value_(valueIn) {}
 
     // cppcheck-suppress functionStatic
-    inline const std::string &value() PROJ_CONST_DEFN { return value_; }
+    inline const std::string &value() PROJ_PURE_DEFN { return value_; }
 
     // cppcheck-suppress functionStatic
-    inline const std::vector<WKTNodeNNPtr> &children() PROJ_CONST_DEFN {
+    inline const std::vector<WKTNodeNNPtr> &children() PROJ_PURE_DEFN {
         return children_;
     }
 
     // cppcheck-suppress functionStatic
-    inline size_t childrenSize() PROJ_CONST_DEFN { return children_.size(); }
+    inline size_t childrenSize() PROJ_PURE_DEFN { return children_.size(); }
 
     // cppcheck-suppress functionStatic
     const WKTNodeNNPtr &lookForChild(const std::string &childName,
diff --git a/src/iso19111/metadata.cpp b/src/iso19111/metadata.cpp
index ddda2987..6bf9d600 100644
--- a/src/iso19111/metadata.cpp
+++ b/src/iso19111/metadata.cpp
@@ -108,7 +108,7 @@ Citation &Citation::operator=(const Citation &other) {
 // ---------------------------------------------------------------------------
 
 /** \brief Returns the name by which the cited resource is known. */
-const optional<std::string> &Citation::title() PROJ_CONST_DEFN {
+const optional<std::string> &Citation::title() PROJ_PURE_DEFN {
     return d->title;
 }
 
@@ -169,7 +169,7 @@ GeographicBoundingBox::~GeographicBoundingBox() = default;
  * If eastBoundLongitude < westBoundLongitude(), then the bounding box crosses
  * the anti-meridian.
  */
-double GeographicBoundingBox::westBoundLongitude() PROJ_CONST_DEFN {
+double GeographicBoundingBox::westBoundLongitude() PROJ_PURE_DEFN {
     return d->west_;
 }
 
@@ -180,7 +180,7 @@ double GeographicBoundingBox::westBoundLongitude() PROJ_CONST_DEFN {
  *
  * The unit is degrees.
  */
-double GeographicBoundingBox::southBoundLatitude() PROJ_CONST_DEFN {
+double GeographicBoundingBox::southBoundLatitude() PROJ_PURE_DEFN {
     return d->south_;
 }
 
@@ -194,7 +194,7 @@ double GeographicBoundingBox::southBoundLatitude() PROJ_CONST_DEFN {
  * If eastBoundLongitude < westBoundLongitude(), then the bounding box crosses
  * the anti-meridian.
  */
-double GeographicBoundingBox::eastBoundLongitude() PROJ_CONST_DEFN {
+double GeographicBoundingBox::eastBoundLongitude() PROJ_PURE_DEFN {
     return d->east_;
 }
 
@@ -205,7 +205,7 @@ double GeographicBoundingBox::eastBoundLongitude() PROJ_CONST_DEFN {
  *
  * The unit is degrees.
  */
-double GeographicBoundingBox::northBoundLatitude() PROJ_CONST_DEFN {
+double GeographicBoundingBox::northBoundLatitude() PROJ_PURE_DEFN {
     return d->north_;
 }
 
@@ -466,19 +466,19 @@ VerticalExtent::~VerticalExtent() = default;
 
 /** \brief Returns the minimum of the vertical extent.
  */
-double VerticalExtent::minimumValue() PROJ_CONST_DEFN { return d->minimum_; }
+double VerticalExtent::minimumValue() PROJ_PURE_DEFN { return d->minimum_; }
 
 // ---------------------------------------------------------------------------
 
 /** \brief Returns the maximum of the vertical extent.
  */
-double VerticalExtent::maximumValue() PROJ_CONST_DEFN { return d->maximum_; }
+double VerticalExtent::maximumValue() PROJ_PURE_DEFN { return d->maximum_; }
 
 // ---------------------------------------------------------------------------
 
 /** \brief Returns the unit of the vertical extent.
  */
-common::UnitOfMeasureNNPtr &VerticalExtent::unit() PROJ_CONST_DEFN {
+common::UnitOfMeasureNNPtr &VerticalExtent::unit() PROJ_PURE_DEFN {
     return d->unit_;
 }
 
@@ -562,13 +562,13 @@ TemporalExtent::~TemporalExtent() = default;
 
 /** \brief Returns the start of the temporal extent.
  */
-const std::string &TemporalExtent::start() PROJ_CONST_DEFN { return d->start_; }
+const std::string &TemporalExtent::start() PROJ_PURE_DEFN { return d->start_; }
 
 // ---------------------------------------------------------------------------
 
 /** \brief Returns the end of the temporal extent.
  */
-const std::string &TemporalExtent::stop() PROJ_CONST_DEFN { return d->stop_; }
+const std::string &TemporalExtent::stop() PROJ_PURE_DEFN { return d->stop_; }
 
 // ---------------------------------------------------------------------------
 
@@ -643,7 +643,7 @@ Extent::~Extent() = default;
  *
  * @return the description, or empty.
  */
-const optional<std::string> &Extent::description() PROJ_CONST_DEFN {
+const optional<std::string> &Extent::description() PROJ_PURE_DEFN {
     return d->description_;
 }
 
@@ -654,7 +654,7 @@ const optional<std::string> &Extent::description() PROJ_CONST_DEFN {
  * @return the geographic element(s), or empty.
  */
 const std::vector<GeographicExtentNNPtr> &
-Extent::geographicElements() PROJ_CONST_DEFN {
+Extent::geographicElements() PROJ_PURE_DEFN {
     return d->geographicElements_;
 }
 
@@ -665,7 +665,7 @@ Extent::geographicElements() PROJ_CONST_DEFN {
  * @return the vertical element(s), or empty.
  */
 const std::vector<VerticalExtentNNPtr> &
-Extent::verticalElements() PROJ_CONST_DEFN {
+Extent::verticalElements() PROJ_PURE_DEFN {
     return d->verticalElements_;
 }
 
@@ -676,7 +676,7 @@ Extent::verticalElements() PROJ_CONST_DEFN {
  * @return the temporal element(s), or empty.
  */
 const std::vector<TemporalExtentNNPtr> &
-Extent::temporalElements() PROJ_CONST_DEFN {
+Extent::temporalElements() PROJ_PURE_DEFN {
     return d->temporalElements_;
 }
 
@@ -980,7 +980,7 @@ Identifier::createFromDescription(const std::string &descriptionIn) {
  *
  * @return the citation for the authority, or empty.
  */
-const optional<Citation> &Identifier::authority() PROJ_CONST_DEFN {
+const optional<Citation> &Identifier::authority() PROJ_PURE_DEFN {
     return d->authority_;
 }
 
@@ -993,7 +993,7 @@ const optional<Citation> &Identifier::authority() PROJ_CONST_DEFN {
  *
  * @return the code.
  */
-const std::string &Identifier::code() PROJ_CONST_DEFN { return d->code_; }
+const std::string &Identifier::code() PROJ_PURE_DEFN { return d->code_; }
 
 // ---------------------------------------------------------------------------
 
@@ -1004,7 +1004,7 @@ const std::string &Identifier::code() PROJ_CONST_DEFN { return d->code_; }
  *
  * @return the authority codespace, or empty.
  */
-const optional<std::string> &Identifier::codeSpace() PROJ_CONST_DEFN {
+const optional<std::string> &Identifier::codeSpace() PROJ_PURE_DEFN {
     return d->codeSpace_;
 }
 
@@ -1017,7 +1017,7 @@ const optional<std::string> &Identifier::codeSpace() PROJ_CONST_DEFN {
  *
  * @return the version or empty.
  */
-const optional<std::string> &Identifier::version() PROJ_CONST_DEFN {
+const optional<std::string> &Identifier::version() PROJ_PURE_DEFN {
     return d->version_;
 }
 
@@ -1028,7 +1028,7 @@ const optional<std::string> &Identifier::version() PROJ_CONST_DEFN {
  *
  * @return the description or empty.
  */
-const optional<std::string> &Identifier::description() PROJ_CONST_DEFN {
+const optional<std::string> &Identifier::description() PROJ_PURE_DEFN {
     return d->description_;
 }
 
@@ -1038,7 +1038,7 @@ const optional<std::string> &Identifier::description() PROJ_CONST_DEFN {
  *
  * @return the URI or empty.
  */
-const optional<std::string> &Identifier::uri() PROJ_CONST_DEFN {
+const optional<std::string> &Identifier::uri() PROJ_PURE_DEFN {
     return d->uri_;
 }
 
@@ -1266,7 +1266,7 @@ PositionalAccuracy::~PositionalAccuracy() = default;
 
 /** \brief Return the value of the positional accuracy.
  */
-const std::string &PositionalAccuracy::value() PROJ_CONST_DEFN {
+const std::string &PositionalAccuracy::value() PROJ_PURE_DEFN {
     return d->value_;
 }