Index: /branches/eam_branches/ipp-20120405/Ohana/src/libohana/include/ohana_sort.h
===================================================================
--- /branches/eam_branches/ipp-20120405/Ohana/src/libohana/include/ohana_sort.h	(revision 33887)
+++ /branches/eam_branches/ipp-20120405/Ohana/src/libohana/include/ohana_sort.h	(revision 33888)
@@ -73,4 +73,6 @@
 void fsortthree (float *X, float *Y, float *Z, int N);
 void dsortthree (double *X, double *Y, double *Z, int N);
+void fsortfour (float *X, float *Y, float *Z, float *W, int N);
+void dsortfour (double *X, double *Y, double *Z, double *W, int N);
 
 #endif
Index: /branches/eam_branches/ipp-20120405/Ohana/src/libohana/src/sorts.c
===================================================================
--- /branches/eam_branches/ipp-20120405/Ohana/src/libohana/src/sorts.c	(revision 33887)
+++ /branches/eam_branches/ipp-20120405/Ohana/src/libohana/src/sorts.c	(revision 33888)
@@ -120,2 +120,36 @@
 
 }
+
+void fsortfour (float *X, float *Y, float *Z, float *W, int N) {
+
+# define SWAPFUNC(A,B){ float tmp; \
+  tmp = X[A]; X[A] = X[B]; X[B] = tmp; \
+  tmp = Y[A]; Y[A] = Y[B]; Y[B] = tmp; \
+  tmp = Z[A]; Z[A] = Z[B]; Z[B] = tmp; \
+  tmp = W[A]; W[A] = W[B]; W[B] = tmp; \
+}
+# define COMPARE(A,B)(X[A] < X[B])
+
+  OHANA_SORT (N, COMPARE, SWAPFUNC);
+
+# undef SWAPFUNC
+# undef COMPARE
+
+}
+
+void dsortfour (double *X, double *Y, double *Z, double *W, int N) {
+
+# define SWAPFUNC(A,B){ double tmp; \
+  tmp = X[A]; X[A] = X[B]; X[B] = tmp; \
+  tmp = Y[A]; Y[A] = Y[B]; Y[B] = tmp; \
+  tmp = Z[A]; Z[A] = Z[B]; Z[B] = tmp; \
+  tmp = W[A]; W[A] = W[B]; W[B] = tmp; \
+}
+# define COMPARE(A,B)(X[A] < X[B])
+
+  OHANA_SORT (N, COMPARE, SWAPFUNC);
+
+# undef SWAPFUNC
+# undef COMPARE
+
+}
Index: /branches/eam_branches/ipp-20120405/Ohana/src/opihi/cmd.data/vstats.c
===================================================================
--- /branches/eam_branches/ipp-20120405/Ohana/src/opihi/cmd.data/vstats.c	(revision 33887)
+++ /branches/eam_branches/ipp-20120405/Ohana/src/opihi/cmd.data/vstats.c	(revision 33888)
@@ -144,5 +144,5 @@
       }      
     }
-    stdev = sqrt (var / N);
+    stdev = sqrt (var / (N - 1));
 
     Nmode = 0;