WSTester updated to work plus hopefully all the other changes that need to go into...

[jabaws.git] / binaries / src / ViennaRNA / RNAforester / g2-0.70 / src / g2_util.c

diff --git a/binaries/src/ViennaRNA/RNAforester/g2-0.70/src/g2_util.c b/binaries/src/ViennaRNA/RNAforester/g2-0.70/src/g2_util.c
new file mode 100644 (file)
index 0000000..00f04e6
--- /dev/null
+++ b/binaries/src/ViennaRNA/RNAforester/g2-0.70/src/g2_util.c
@@ -0,0 +1,224 @@
+/*****************************************************************************
+**  Copyright (C) 1998-2001  Ljubomir Milanovic & Horst Wagner
+**  This file is part of the g2 library
+**
+**  This library is free software; you can redistribute it and/or
+**  modify it under the terms of the GNU Lesser General Public
+**  License as published by the Free Software Foundation; either
+**  version 2.1 of the License, or (at your option) any later version.
+**
+**  This library is distributed in the hope that it will be useful,
+**  but WITHOUT ANY WARRANTY; without even the implied warranty of
+**  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+**  Lesser General Public License for more details.
+**
+**  You should have received a copy of the GNU Lesser General Public
+**  License along with this library; if not, write to the Free Software
+**  Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
+******************************************************************************/
+#include <stdio.h>
+#include <stdlib.h>
+#include <math.h>
+#include <stdarg.h>
+
+#include "g2_util.h"
+#include "g2_physical_device.h"
+#include "g2_config.h"
+
+
+/*
+ *
+ * Double to integer
+ *
+ */
+int dtoi(double x)
+{
+    return (int)(x+0.5);
+}
+
+
+
+/*
+ *
+ * Transform user coord. in physical device coord (int)
+ *
+ */
+void g2_uc2pdc_int(g2_physical_device *pd, double x, double y,
+                  int *ix, int *iy)
+{
+    double vx, vy;
+    
+    vx=pd->x_mul*x+pd->x_origin;
+    vy=pd->y_mul*y+pd->y_origin;
+    
+    *ix = dtoi(pd->a11*vx + pd->b1);
+    *iy = dtoi(pd->a22*vy + pd->b2);
+}
+
+
+/*
+ *
+ * Transform user coord. in physical device coord (double)
+ *
+ */
+void g2_uc2pdc_double(g2_physical_device *pd, double x, double y,
+                     double *dx, double *dy)
+{
+    double vx, vy;
+    
+    vx=pd->x_mul*x+pd->x_origin;
+    vy=pd->y_mul*y+pd->y_origin;
+    
+    *dx = pd->a11*vx + pd->b1;
+    *dy = pd->a22*vy + pd->b2;
+}
+
+
+/*
+ *
+ * Transform user size in physical device size (int)
+ *
+ */
+void g2_us2pds_int(g2_physical_device *pd, double x, double y,
+                  int *ix, int *iy)
+{
+    if(ix!=NULL)
+       *ix=dtoi(x*fabs(pd->x_mul*pd->a11));
+    if(iy!=NULL)
+       *iy=dtoi(y*fabs(pd->y_mul*pd->a22));
+}
+
+
+/*
+ *
+ * Transform user size in physical device size (double)
+ *
+ */
+void g2_us2pds_double(g2_physical_device *pd, double x, double y,
+                     double *dx, double *dy)
+{
+    if(dx!=NULL)
+       *dx=x*fabs(pd->x_mul*pd->a11);
+    if(dy!=NULL)
+       *dy=y*fabs(pd->y_mul*pd->a22);
+}
+
+
+/*
+ *
+ * Transform physical device coord in user coord
+ *
+ */
+void g2_pdc2uc(g2_physical_device *pd, double ix, double iy,
+                  double *x, double *y)
+{
+    double pcx, pcy;
+    pcx=(ix-pd->b1)/pd->a11;
+    pcy=(iy-pd->b2)/pd->a22;
+
+    *x=(pcx-pd->x_origin)/pd->x_mul;
+    *y=(pcy-pd->y_origin)/pd->y_mul;
+}
+
+
+/*
+ *  return a < b
+ */
+void g2_sort2_i(int *a, int *b)
+{
+    if(*a>*b) {
+       int t=*a;
+       *a=*b; *b=t;
+    }
+}
+
+void g2_sort2_d(double *a, double *b)
+{
+    if(*a>*b) {
+       double t=*a;
+       *a=*b; *b=t;
+    }
+}
+
+
+
+/*
+ *
+ * g2 malloc (with error message)
+ *
+ */
+void *g2_malloc(size_t size)
+{
+    void *rv;
+
+    if((rv=malloc(size))==NULL) {
+       fprintf(stderr, "g2_malloc: Can not allocate memory\n");
+       exit(-1);
+    }
+
+    return rv;
+}
+
+
+/*
+ *
+ * g2 realloc (with error message)
+ *
+ */
+void *g2_realloc(void *p, size_t size)
+{
+    void *rv;
+    if((rv=realloc(p, size))==NULL) {
+       fprintf(stderr, "g2_realloc: Can not allocate memory\n");
+       exit(-1);
+    }
+    return rv;
+}
+
+
+
+/*
+ *
+ * g2 free
+ *
+ */
+void g2_free(void *p)
+{
+    if(p!=NULL)
+       free(p);
+}
+
+
+/*
+ *
+ * transform float* to double* for N elements
+ *
+ * Note: don't forget to free d
+ *
+ */
+double *g2_floatp2doublep(float *f, int N)
+{
+    int i;
+    double *d;
+    d=(double *)g2_malloc(N*sizeof(double));
+    for(i=0;i<N;++i)
+       d[i]=f[i];
+    return d;
+}
+
+
+/*
+ *
+ * log messages to stderr
+ *
+ */
+void g2_log(enum g2_log_level log_level, const char *format, ...)
+{
+    va_list arg;
+    if(log_level > g2_LogLevel) {
+       return;
+    }
+    va_start(arg, format);
+    vfprintf(stderr, format, arg);
+    va_end(arg);    
+}