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(size_t, size_t);
  13 extern void la_free_alloc(PTR);
  14 extern void xlfail(char *);
  15
  16 /************************************************************************/
  17 /**                                                                    **/
  18 /**                 Storage Allocation Functions                       **/
  19 /**                                                                    **/
  20 /************************************************************************/
  21
  22 static char
  23 *allocate(size_t n, size_t m)
  24 {
  25   char *p = (char *) la_allocate(n, m);
  26   if (p == nil) xlfail("allocation failed");
  27   return(p);
  28 }
  29
  30 static void
  31 free_alloc(char *p)
  32 {
  33   if (p != nil) la_free_alloc((PTR) p);
  34 }
  35
  36 IVector
  37 ivector(size_t n)
  38 {
  39   return((IVector) allocate(n, sizeof(int)));
  40 }
  41
  42 double
  43 *rvector(size_t n)
  44 {
  45   return((double *) allocate(n, sizeof(double)));
  46 }
  47
  48 CVector cvector(size_t n)
  49 {
  50   return((CVector) allocate(n, sizeof(Complex)));
  51 }
  52
  53 void
  54 free_vector(double *v)
  55 {
  56   free_alloc((char *)v);
  57 }
  58
  59 IMatrix imatrix(size_t n, size_t m)
  60 {
  61   size_t i;
  62   IMatrix mat = (IMatrix) allocate(n, sizeof(IVector));
  63   for (i = 0; i < n; i++) mat[i] = (IVector) allocate(m, sizeof(int));
  64   return(mat);
  65 }
  66
  67 RMatrix
  68 rmatrix(size_t n, size_t m)
  69 {
  70   size_t i;
  71   RMatrix mat = (RMatrix) allocate(n, sizeof(RVector));
  72   for (i = 0; i < n; i++) mat[i] = (RVector) allocate(m, sizeof(double));
  73   return(mat);
  74 }
  75
  76 CMatrix
  77 cmatrix(size_t n, size_t m)
  78 {
  79   size_t i;
  80   CMatrix mat = (CMatrix) allocate(n, sizeof(CVector));
  81   for (i = 0; i < n; i++) mat[i] = (CVector) allocate(m, sizeof(Complex));
  82   return(mat);
  83 }
  84
  85 void
  86 free_matrix(double **mat, int n) /* Matrix?? Not RMatrix?? */
  87 {
  88   size_t i;
  89
  90   if (mat != nil) for (i = 0; i < n; i++) free_alloc((char *)mat[i]);
  91   free_alloc((char *)mat);
  92 }