add getAngle judge

GWmodel-Lab · Jan 8, 2025 · fde6bb4 · fde6bb4
1 parent 7ec996e
commit fde6bb4
Show file tree

Hide file tree

Showing 2 changed files with 6 additions and 1 deletion.
diff --git a/include/gwmodelpp/GTWR.h b/include/gwmodelpp/GTWR.h
@@ -86,7 +86,11 @@ class GTWR :  public GWRBase, public IBandwidthSelectable, public IParallelizabl
     }
 
     const double getAngle(){
-        return mStdistance->angle();
+        if (mStdistance != nullptr) {
+            return mStdistance->angle();
+        } else {
+            throw std::runtime_error("mStdistance is not initialized");
+        }
     }
 
 private:

diff --git a/test/testGTWR.cpp b/test/testGTWR.cpp
@@ -80,6 +80,7 @@ TEST_CASE("GTWR: londonhp100")
         REQUIRE(algorithm.hasIntercept() == true);
         ////lambda=1
         REQUIRE_THAT(algorithm.getLambda(), Catch::Matchers::WithinAbs(1.0, 1e-12));
+        REQUIRE_THAT(algorithm.getAngle(), Catch::Matchers::WithinAbs(arma::datum::pi / 2, 1e-12));
         REQUIRE_THAT(diagnostic.AIC, Catch::Matchers::WithinAbs(2436.6044573089, 1e-8));
         REQUIRE_THAT(diagnostic.AICc, Catch::Matchers::WithinAbs(2448.2720652516, 1e-8));
         REQUIRE_THAT(diagnostic.RSquare, Catch::Matchers::WithinAbs(0.7080106320292, 1e-8));