Hi Hassan,
It appears that your BH_diagnostics.ah2.gp file is missing a single iteration (41440). My guess would be that this is because the apparent horizon finder failed to find a horizon at this iteration.
This causes trouble when trying to combine the two individual puncture trajectories into a single separation vector. Perhaps SimulationTools should somehow be more tolerant of this scenario, although my test shows to enabling resampling works around it. Alternatively, you can just work with the individual puncture trajectories with something like the following:
(* Read individual trajectories *) traj1 = ReadBinaryCoordinates["test", 1]; traj2 = ReadBinaryCoordinates["test", 2];
(* Cut out the bad bit at the end *)
traj1 = Slab[#, 0 ;; 121.312] & /@ traj1; traj2 = Slab[#, 0 ;; 121.312] & /@ traj2;
(* Compute coordinate separation *) sep = Norm[traj1 - traj2];
I hope this helps, Barry
On Mon, 9 Sep 2019 at 15:09, Hassan Khalvati hassan.kh92@gmail.com wrote:
Hello, I am using simulation tools in Mathematica to extract and visualize my simulation results and thanks to Ian, and Barry, it is really helpful. but I have just faced a problem in checking the binary separation. that is when I am asking to ReadBinarySeparation from my simulations, it is posting an error:
"Plus cannot operate on DataTables with different coordinates. To allow this, use the WithResampling function."
and I have tested for 3 different simulation results and all have the same error message. But if I am replacing the "BH_diagnostics.ah1.gp" file from the simulation tools' TestSimulation, I am getting output, however that output does not actually belong to my data, the point is that there is no error in this case.
And when I am using "WithResampling" function, I think some of my data is missing, because my data table is reducing to 990, from about 27000.
I have attached the diagnostics files from both, my results, and "Testsimulation", in case it is needed. I would be grateful if anyone could help me with this. With regards, Hassan
--
Hassan Khalvati Sharif University of Tehran
Users mailing list Users@einsteintoolkit.org http://lists.einsteintoolkit.org/mailman/listinfo/users