Changeset 39245

Timestamp:

Dec 11, 2015, 9:12:28 AM (11 years ago)

Author:

eugene

Message:

add api for arbitrary 3D rotations

Location:

trunk/Ohana/src

Files:

• libdvo/include/libdvo_astro.h (modified) (1 diff)
• libdvo/src/coord_systems.c (modified) (1 diff)
• opihi/cmd.astro/Makefile (modified) (1 diff)
• opihi/cmd.astro/crotation.c (added)
• opihi/cmd.astro/init.c (modified) (2 diffs)
• opihi/cmd.data/gaussdeviate.c (modified) (2 diffs)

trunk/Ohana/src/libdvo/include/libdvo_astro.h

@@ -138 +138 @@
 } Image;
 
+CoordTransform *AllocTransform (double phi, double Xo, double xo);
 CoordTransform *InitTransform (CoordTransformSystem input, CoordTransformSystem output);
 int ApplyTransform (double *x, double *y, double X, double Y, CoordTransform *transform);

trunk/Ohana/src/libdvo/src/coord_systems.c

@@ -1 +1 @@
 # include "dvo.h"
     
+CoordTransform *AllocTransform (double phi, double Xo, double xo) {
+
+  CoordTransform *transform;
+  ALLOCATE (transform, CoordTransform, 1);
+  transform->isIdentity = FALSE;
+  
+  transform->phi = phi*RAD_DEG;
+  transform->Xo  = Xo*RAD_DEG;
+  transform->xo  = xo;
+
+  // pre-calculated constants:
+  transform->sin_phi_cos_Xo = sin(transform->phi)*cos(transform->Xo);
+  transform->sin_phi_sin_Xo = sin(transform->phi)*sin(transform->Xo);
+  transform->cos_phi        = cos(transform->phi);
+  
+  transform->cos_phi_cos_Xo = cos(transform->phi)*cos(transform->Xo);
+  transform->cos_phi_sin_Xo = cos(transform->phi)*sin(transform->Xo);
+  transform->sin_phi        = sin(transform->phi);
+  transform->cos_Xo         = cos(transform->Xo);
+  transform->sin_Xo         = sin(transform->Xo);
+
+  return transform;
+}
+
 CoordTransform *InitTransform (CoordTransformSystem input, CoordTransformSystem output) {
 

trunk/Ohana/src/opihi/cmd.astro/Makefile

@@ -25 +25 @@
 $(SRC)/cneedles.$(ARCH).o                  \
 $(SRC)/cplot.$(ARCH).o             \
+$(SRC)/crotation.$(ARCH).o         \
 $(SRC)/csystem.$(ARCH).o           \
 $(SRC)/ctimes.$(ARCH).o    \

trunk/Ohana/src/opihi/cmd.astro/init.c

@@ -9 +9 @@
 int cneedles                PROTO((int, char **));
 int cplot                   PROTO((int, char **));
+int crotation               PROTO((int, char **));
 int csystem                 PROTO((int, char **));
 int ctimes                  PROTO((int, char **));
@@ -78 +79 @@
   {1, "cneedles",    cneedles,     "plot vectors in sky coordinates"},
   {1, "cplot",       cplot,        "plot vectors in sky coordinates"},
+  {1, "crotation",   crotation,    "rotate in 3D"},
   {1, "csystem",     csystem,      "convert between coordinate systems"},
   {1, "ctimes",      ctimes,       "convert between time formats"},

trunk/Ohana/src/opihi/cmd.data/gaussdeviate.c

@@ -31 +31 @@
 int gaussintegral (int argc, char **argv) {
   
+  gprint (GP_ERR, "fix gaussintegral\n");
+  return (FALSE);
+
+/*
   int i, Npts;
   Vector *vec;
@@ -53 +57 @@
   gprint (GP_ERR, "USAGE: gaussintegral Npts mean sigma\n");
   return (FALSE);
+*/
     
 }