Changeset 14709 for trunk/FACT++/src/InterpreterV8.cc

Timestamp:

12/01/12 18:59:31 (12 years ago)

Author:

tbretz

Message:

added dist function to Sky class

File:

• trunk/FACT++/src/InterpreterV8.cc (modified) (4 diffs)

trunk/FACT++/src/InterpreterV8.cc

@@ -1250 +1250 @@
 }
 
-Handle<Value> InterpreterV8::LocalDist(const Arguments &args)
+Handle<Value> InterpreterV8::CalcDist(const Arguments &args, const bool local)
 {
     if (args.Length()!=2)
@@ -1266 +1266 @@
     };
 
-    const Handle<String> s_zd = String::New("zd");
-    const Handle<String> s_az = String::New("az");
-
-    const double zd0 = obj[0]->Get(s_zd)->NumberValue() * M_PI/180;
-    const double az0 = obj[0]->Get(s_az)->NumberValue() * M_PI/180;
-    const double zd1 = obj[1]->Get(s_zd)->NumberValue() * M_PI/180;
-    const double az1 = obj[1]->Get(s_az)->NumberValue() * M_PI/180;
-
-    if (!finite(zd0) || !finite(zd1) || !finite(az0) || !finite(az1))
+    const Handle<String> s_theta = String::New(local?"zd":"dec"); // was: zd
+    const Handle<String> s_phi   = String::New(local?"az":"ra");  // was: az
+
+    const double conv_t = M_PI/180;
+    const double conv_p = local ? -M_PI/180 : M_PI/12;
+    const double offset = local ? 0 : M_PI;
+
+    const double theta0 = offset - obj[0]->Get(s_theta)->NumberValue() * conv_t;
+    const double phi0   =          obj[0]->Get(s_phi  )->NumberValue() * conv_p;
+    const double theta1 = offset - obj[1]->Get(s_theta)->NumberValue() * conv_t;
+    const double phi1   =          obj[1]->Get(s_phi  )->NumberValue() * conv_p;
+
+    if (!finite(theta0) || !finite(theta1) || !finite(phi0) || !finite(phi1))
         return ThrowException(String::New("some values not valid or not finite."));
 
@@ -1291 +1295 @@
     // cos(az1-az0) = cos(az1)*cos(az0) + sin(az1)*sin(az0)
 
-    const double x = sin(zd0) * sin(zd1) * cos(az1-az0);
-    const double y = cos(zd0) * cos(zd1);
+    const double x = sin(theta0) * sin(theta1) * cos(phi1-phi0);
+    const double y = cos(theta0) * cos(theta1);
 
     const double res = acos(x + y) * 180/M_PI;
 
     return handle_scope.Close(Number::New(res));
+}
+
+Handle<Value> InterpreterV8::LocalDist(const Arguments &args)
+{
+    return CalcDist(args, true);
+}
+
+Handle<Value> InterpreterV8::SkyDist(const Arguments &args)
+{
+    return CalcDist(args, false);
 }
 
@@ -1970 +1984 @@
     Handle<FunctionTemplate> sky = FunctionTemplate::New(ConstructorSky);
     sky->SetClassName(String::New("Sky"));
+    sky->Set(String::New("dist"),  FunctionTemplate::New(SkyDist), ReadOnly);
     global->Set(String::New("Sky"), sky, ReadOnly);