Atilde and Ricci check

Author: at0m741

Tests/Matrix/Matrix.cpp

@@ -294,7 +294,7 @@ int matrix_tests() {
 
 	tensorium::Vector<double> X(dim);
 	X(0) = 0.0; 
-	X(1) = 10.0;
+	X(1) = 3.0;
 	X(2) = M_PI / 2.0;
 	X(3) = 0.0;
 
@@ -303,7 +303,7 @@ int matrix_tests() {
 
 	std::cout << "X = " << X(0) << " " << X(1) << " " << X(2) << " " << X(3) << "\n";
 
-	tensorium_RG::Metric<double> metric("kerr_schild", 1.0, 0.9);
+	tensorium_RG::Metric<double> metric("kerr_schild", 1.0, 0.935);
 	metric(X, g);
 
 	std::cout << "Metric tensor g at X = (t=0, r=10, θ=π/2, φ=0):\n";
@@ -313,7 +313,7 @@ int matrix_tests() {
 	std::cout << "Christoffel symbols Γ^λ_{μν} at X = (t=0, r=10, θ=π/2, φ=0):\n";
 	auto gamma = tensorium::compute_christoffel(X, 1e-5, g, g_inv, metric);
 	gamma.print();
-	auto R = tensorium::compute_riemann_tensor<double>(X, 1e-5, tensorium_RG::Metric<double>("kerr_schild", 1.0, 0.9));
+	auto R = tensorium::compute_riemann_tensor<double>(X, 1e-5, tensorium_RG::Metric<double>("kerr_schild", 1.0, 0.935));
 	tensorium::print_riemann_tensor(R);
 	tensorium::contract_tensor<0, 1>(R);
 	std::cout << "Riemann tensor contracted:\n";
@@ -324,7 +324,7 @@ int matrix_tests() {
 	tensorium::Vector<double> beta(3);
 	tensorium::Tensor<double, 2> gammaj({3, 3});
 	metric.BSSN(X, alpha, beta, gammaj);
-
+ 
 	tensorium::Tensor<double, 2> gammaj_inv = tensorium::inv_mat_tensor(gammaj);
 
 	std::cout << "\n--- BSSN 3+1 Decomposition ---\n";

includes/Tensorium/DiffGeometry/BSSN/BSSNSetup.hpp

@@ -147,7 +147,7 @@ template <typename T> class BSSN {
         tensorium_RG::BSSNAtildeTensor<T> Aij;
         auto AtildeTensor = Aij.compute_Atilde_tensor(Kij, gamma_ij_inv, gamma_ij, chi);
 
-        const size_t NX = 1, NY = 1, NZ = 1;
+        const size_t NX = 32, NY = 32, NZ = 32;
 
         tensorium::Tensor<T, 5> Atilde_full({NX, NY, NZ, 3, 3});
         tensorium::Tensor<T, 5> gtilde_inv_full({NX, NY, NZ, 3, 3});
@@ -167,9 +167,9 @@ template <typename T> class BSSN {
         }
 
         auto psi_full = ConstraintSolver<T>::solveLichnerowicz(Atilde_full, gtilde_inv_full, dx, dy,
-                                                               dz, 2000, T(1e-8));
+                                                               dz, 4000, T(1e-4));
         {
-            T psi_000 = std::fmax(psi_full(std::array<size_t, 3>{0, 0, 0}), T(1e-8));
+            T psi_000 = std::fmax(psi_full(std::array<size_t, 3>{0, 0, 0}), T(1e-4));
             T new_chi = T(1) / std::pow(psi_000, T(4));
             grid.chi = {new_chi};
         }
@@ -188,7 +188,7 @@ template <typename T> class BSSN {
         grid.beta = {beta};
         grid.gamma_ij = {gamma_ij};
         grid.gamma_ij_inv = {gamma_ij_inv};
-        grid.chi = {chi};
+        // grid.chi = {new_chi};
         grid.gamma_tilde = {gamma_tilde};
         grid.gamma_tilde_inv = {gamma_tilde_inv};
         grid.dgamma_tilde = {dgamma_tilde};

shell.nix

@@ -12,7 +12,7 @@ pkgs.mkShell {
     gcc
     openblas
     openmpi
-    valgrind
+    # valgrind
     cloc
     tree
 	doxygen 
@@ -30,7 +30,7 @@ pkgs.mkShell {
       pybind11
     ]))
 
-  ] ++ (with llvmPackages_18; [
+  ] ++ (with llvmPackages_19; [
     mlir
     clang
     llvm