[Debian-astro-commits] [gyoto] 12/221: Correct also the generic implementation of circularVelocity()
Thibaut Jean-Claude Paumard
thibaut at moszumanska.debian.org
Fri May 22 20:52:29 UTC 2015
This is an automated email from the git hooks/post-receive script.
thibaut pushed a commit to branch master
in repository gyoto.
commit f9f33402a22908cd25d1d0b71ac6a40f2317489d
Author: Thibaut Paumard <paumard at users.sourceforge.net>
Date: Tue Sep 9 11:22:24 2014 +0200
Correct also the generic implementation of circularVelocity()
---
lib/Metric.C | 2 +-
lib/PolishDoughnut.C | 4 ++--
2 files changed, 3 insertions(+), 3 deletions(-)
diff --git a/lib/Metric.C b/lib/Metric.C
index 70fbd05..80d5404 100644
--- a/lib/Metric.C
+++ b/lib/Metric.C
@@ -329,7 +329,7 @@ void Metric::Generic::circularVelocity(double const * coor, double* vel,
double coord[4] = {coor[0], coor[1]*sinth, M_PI*0.5, coor[3]};
vel[1] = vel[2] = 0.;
- vel[3] = 1./(dir*pow(coord[1], 1.5)*sinth);
+ vel[3] = 1./(dir*pow(coord[1], 1.5));
vel[0] = SysPrimeToTdot(coor, vel+1);
vel[3] *= vel[0];
diff --git a/lib/PolishDoughnut.C b/lib/PolishDoughnut.C
index ee00b20..5ce6e37 100644
--- a/lib/PolishDoughnut.C
+++ b/lib/PolishDoughnut.C
@@ -333,8 +333,8 @@ double PolishDoughnut::operator()(double const coord[4]) {
void PolishDoughnut::getVelocity(double const pos[4], double vel[4])
{
if (adaf_) {
- // This will return the Keplerian velocity at the
- // radius projected on the equat plane
+ // This will return the circular velocity at the
+ // radius projected on the equat plane, or it's Keplerian approximation
return gg_->circularVelocity(pos,vel,1);
}
double gtt=gg_->gmunu(pos,0,0);
--
Alioth's /usr/local/bin/git-commit-notice on /srv/git.debian.org/git/debian-astro/packages/gyoto.git
More information about the Debian-astro-commits
mailing list