On 7 Nov 2017, at 17:37, Nisa Amir <nisaa...@math.qau.edu.pk> wrote:
> I am just trying to add this spacetime, I dont want to use it in some other
> thorns of the toolkit. and I am not familiar with xAct package.
Hi,
The below patch modifies EinsteinExact so that it is able to generate thorns
from metrics which are in spherical polar coordinates. It achieves this by
mapping the polar coordinate r,th,ph to the CartGrid3D x,y,z required by
ADMBase. You can use this to generate the existing Schwarzschild metric from
the Metrics package. Note that this means that your x, y, z coordinates in the
simulation are really r, th and ph, and that much of the rest of the toolkit
will not know how to deal with this, but if you really don't want to use any
other thorns, then that shouldn't be a problem.
diff --git a/m/EinsteinExact.m b/m/EinsteinExact.m
index 80eaff2..93b7fff 100644
--- a/m/EinsteinExact.m
+++ b/m/EinsteinExact.m
@@ -247,11 +247,12 @@ idThorn[spacetime_, thorn_] :=
Print["Generating thorn ", thorn, " for ", spacetime, " spacetime."];
(* Load the spacetime: coordinates, metric, inverse metric *)
- coordRule = {t -> T, x -> X, y -> Y, z -> Z, None -> {}};
+ coordRule = {t -> T, x -> X, y -> Y, z -> Z, None -> {},
+ r -> X, \[Theta] -> Y, \[Phi] -> Z};
coords = MetricProperty[spacetime, "Coordinates"] /. coordRule;
If[coords =!= {T, X, Y, Z},
- Throw["Error, only metrics in Cartesian coordinates are supported"];
+ Throw["Error, only metrics in Cartesian or spherical polar coordinates are
supported"];
];
spatialCoords = coords[[2;;]];
@@ -456,7 +457,8 @@ g_{ab} =
StringForm[doctext]
]
-spacetimes = {"GaugeWave", "KerrSchild", "Minkowski", "ShiftedGaugeWave",
"Vaidya", "ModifiedSchwarzschildBL"};
+(* spacetimes = {"GaugeWave", "KerrSchild", "Minkowski", "ShiftedGaugeWave",
"Vaidya", "ModifiedSchwarzschildBL"}; *)
+spacetimes = {"Schwarzschild"};
thornNameRules = {"Vaidya" -> "Vaidya2"};
thorns = spacetimes /. thornNameRules;
--
Ian Hinder
http://members.aei.mpg.de/ianhin
_______________________________________________
Users mailing list
Users@einsteintoolkit.org
http://lists.einsteintoolkit.org/mailman/listinfo/users
