lib/linalgdata.c

   1 /* linalgdata - allocation support for basic linear algebra routines.  */
   2 /* Copyright (c) 1990, by Luke Tierney                                 */
   3
   4 #include "linalg.h"
   5
   6 #ifdef INTPTR
   7 typedef int PTR;
   8 #else
   9 typedef char *PTR;
  10 #endif
  11
  12 extern PTR la_allocate(int, int);
  13 extern void la_free_alloc(PTR);
  14 extern void xlfail(char *);
  15
  16
  17 /************************************************************************/
  18 /**                                                                    **/
  19 /**                 Storage Allocation Functions                       **/
  20 /**                                                                    **/
  21 /************************************************************************/
  22
  23 static char
  24 *allocate(unsigned n, unsigned m)
  25 {
  26   char *p = (char *) la_allocate(n, m);
  27   if (p == nil) xlfail("allocation failed");
  28   return(p);
  29 }
  30
  31 static void
  32 free_alloc(char *p)
  33 {
  34   if (p != nil) la_free_alloc((PTR) p);
  35 }
  36
  37 IVector
  38 ivector(unsigned n)
  39 {
  40   return((IVector) allocate(n, sizeof(int)));
  41 }
  42
  43 double
  44 *rvector(unsigned n)
  45 {
  46   return((double *) allocate(n, sizeof(double)));
  47 }
  48
  49 CVector cvector(n)
  50         unsigned n;
  51 {
  52   return((CVector) allocate(n, sizeof(Complex)));
  53 }
  54
  55 void
  56 free_vector(double *v)
  57 {
  58   free_alloc((char *)v);
  59 }
  60
  61 IMatrix imatrix(unsigned n, unsigned m)
  62 {
  63   int i;
  64   IMatrix mat = (IMatrix) allocate(n, sizeof(IVector));
  65   for (i = 0; i < n; i++) mat[i] = (IVector) allocate(m, sizeof(int));
  66   return(mat);
  67 }
  68
  69 RMatrix
  70 rmatrix(unsigned n, unsigned m)
  71 {
  72   int i;
  73   RMatrix mat = (RMatrix) allocate(n, sizeof(RVector));
  74   for (i = 0; i < n; i++) mat[i] = (RVector) allocate(m, sizeof(double));
  75   return(mat);
  76 }
  77
  78 CMatrix
  79 cmatrix(unsigned n, unsigned m)
  80 {
  81   int i;
  82   CMatrix mat = (CMatrix) allocate(n, sizeof(CVector));
  83   for (i = 0; i < n; i++) mat[i] = (CVector) allocate(m, sizeof(Complex));
  84   return(mat);
  85 }
  86
  87 void
  88 free_matrix(double **mat, int n) /* Matrix?? Not RMatrix?? */
  89 {
  90   int i;
  91
  92   if (mat != nil) for (i = 0; i < n; i++) free_alloc((char *)mat[i]);
  93   free_alloc((char *)mat);
  94 }