[Debian-astro-commits] [gyoto] 12/221: Correct also the generic implementation of circularVelocity()

[Thibaut Jean-Claude Paumard]

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