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