Changeset 26891 for trunk/Ohana/src/opihi/cmd.data

Timestamp:

Feb 10, 2010, 7:24:46 PM (16 years ago)

Author:

eugene

Message:

updates from eam_branches/20091201

Location:

trunk/Ohana/src/opihi/cmd.data

Files:

• Makefile (modified) (2 diffs)
• hermitian1d.c (copied) (copied from branches/eam_branches/20091201/Ohana/src/opihi/cmd.data/hermitian1d.c )
• hermitian2d.c (copied) (copied from branches/eam_branches/20091201/Ohana/src/opihi/cmd.data/hermitian2d.c )
• init.c (modified) (4 diffs)
• matrix.c (copied) (copied from branches/eam_branches/20091201/Ohana/src/opihi/cmd.data/matrix.c )
• svd.c (modified) (1 diff)
• test/matrix.sh (copied) (copied from branches/eam_branches/20091201/Ohana/src/opihi/cmd.data/test/matrix.sh )
• test/svd.sh (copied) (copied from branches/eam_branches/20091201/Ohana/src/opihi/cmd.data/test/svd.sh )

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

@@ -53 +53 @@
 $(SRC)/ungridify.$(ARCH).o     \
 $(SRC)/histogram.$(ARCH).o      \
+$(SRC)/hermitian1d.$(ARCH).o    \
+$(SRC)/hermitian2d.$(ARCH).o    \
 $(SRC)/imcut.$(ARCH).o          \
 $(SRC)/imhist.$(ARCH).o \
@@ -69 +71 @@
 $(SRC)/load.$(ARCH).o           \
 $(SRC)/lookup.$(ARCH).o \
+$(SRC)/matrix.$(ARCH).o \
 $(SRC)/mkrgb.$(ARCH).o  \
 $(SRC)/mcreate.$(ARCH).o        \

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

@@ -42 +42 @@
 int ungridify        PROTO((int, char **));
 int histogram        PROTO((int, char **));
+int hermitian1d      PROTO((int, char **));
+int hermitian2d      PROTO((int, char **));
 int imcut            PROTO((int, char **));
 int imhist           PROTO((int, char **));
@@ -58 +60 @@
 int load             PROTO((int, char **));
 int lookup           PROTO((int, char **));
+int matrix           PROTO((int, char **));
 int mkrgb            PROTO((int, char **));
 int mcreate          PROTO((int, char **));
@@ -172 +175 @@
   {1, "header",       header,           "print image header"},
   {1, "histogram",    histogram,        "generate histogram from vector"},
+  {1, "hermitian1d",  hermitian1d,      "generate 1-D Hermitian Polynomial"},
+  {1, "hermitian2d",  hermitian2d,      "generate 2-D Hermitian Polynomial"},
   {1, "imbin",        rebin,            "rebin image data by factor of N"},
   {1, "imclip",       imclip,           "clip values in an image to be within a range"},
@@ -195 +200 @@
   {1, "minterp",      minterp,          "interpolate image pixels"},
   {1, "iminterp",     minterp,          "interpolate image pixels"},
+  {1, "matrix",       matrix,           "matrix math operations"},
   {1, "mkrgb",        mkrgb,            "convert 3 images to rgb jpeg (use Kapa for better control)"},
   {1, "mset",         mset,             "insert a vector in an image"},

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

@@ -57 +57 @@
   /* use a bcopy instead? */
 
+  // try C.R. Bond's version -- requires matrices in the form A[row][col] not A[row*Ncol + col]
+  if (1) { 
+      int j;
+      double **a, **u, **v, *q;
+      ALLOCATE (a, double *, Ny);
+      ALLOCATE (u, double *, Ny);
+      ALLOCATE (v, double *, Ny);
+      for (i = 0; i < Ny; i++) {
+          ALLOCATE (a[i], double, Nx);
+          ALLOCATE (u[i], double, Nx);
+          ALLOCATE (v[i], double, Nx);
+      }   
+      ALLOCATE (q, double, Nx);
+
+      for (j = 0; j < Ny; j++) {
+          for (i = 0; i < Nx; i++) {
+              a[j][i] = A[j*Nx + i];
+              u[j][i] = 0;
+              v[j][i] = 0;
+          }
+      }
+      for (i = 0; i < Nx; i++) {
+          q[i] = 0;
+      }
+
+      status = svdcmp_bond_raw (Ny, Nx, 1, 1, FLT_EPSILON, 1e-6, a, q, u, v);
+      fprintf (stderr, "status: %d\n", status);
+
+      // copy u q v back to U W V:
+      for (j = 0; j < Ny; j++) {
+          for (i = 0; i < Nx; i++) {
+              U[j*Nx + i] = u[j][i];
+              V[j*Nx + i] = v[j][i];
+          }
+      }
+      for (i = 0; i < Nx; i++) {
+          W[i] = q[i];
+      }
+
+      for (j = 0; j < Ny; j++) {
+          free(a[j]);
+          free(u[j]);
+          free(v[j]);
+      }
+      free (a);
+      free (u);
+      free (v);
+      free (q);
+
+      return TRUE;
+  }
+
   status = svdcmp (U, W, V, Nx, Ny);
   if (!status) {