zsolve/vector.c

   1 /*
   2 4ti2 -- A software package for algebraic, geometric and combinatorial
   3 problems on linear spaces.
   4
   5 Copyright (C) 2006 4ti2 team.
   6 Main author(s): Matthias Walter.
   7
   8 This program is free software; you can redistribute it and/or
   9 modify it under the terms of the GNU General Public License
  10 as published by the Free Software Foundation; either version 2
  11 of the License, or (at your option) any later version.
  12
  13 This program is distributed in the hope that it will be useful,
  14 but WITHOUT ANY WARRANTY; without even the implied warranty of
  15 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
  16 GNU General Public License for more details.
  17
  18 You should have received a copy of the GNU General Public License
  19 along with this program; if not, write to the Free Software
  20 Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
  21 */
  22
  23 #include "vector.h"
  24
  25 #include <assert.h>
  26 #include <stdlib.h>
  27 #include "defs.h"
  28
  29 //                                                                            //
  30
  31 Vector createVector(int size)
  32 {
  33         Vector vector = (Vector)malloc(size*sizeof(vector_t));
  34
  35         if (vector==NULL)
  36         {
  37                 fprintf(stderr, "Fatal Error (%s/%d): Could not allocate memory for Vector!\n", __FILE__, __LINE__);
  38                 exit(1);
  39         }
  40
  41         return vector;
  42 }
  43 //                                                                            //
  44
  45 Vector createZeroVector(int size)
  46 {
  47         Vector vector = (Vector)calloc(size, sizeof(vector_t));
  48
  49         if (vector==NULL)
  50         {
  51                 fprintf(stderr, "Fatal Error (%s/%d): Could not allocate memory for Vector!\n", __FILE__, __LINE__);
  52                 exit(1);
  53         }
  54
  55         return vector;
  56 }
  57
  58 //                                                                            //
  59
  60 inline void deleteVector(Vector vector)
  61 {
  62         free(vector);
  63 }
  64
  65 //                                                                            //
  66
  67 Vector copyVector(Vector old, int size)
  68 {
  69         Vector new = (Vector)malloc(size*sizeof(vector_t));
  70
  71         if (new==NULL)
  72         {
  73                 fprintf(stderr, "Fatal Error (%s/%d): Could not allocate memory for Vector!\n", __FILE__, __LINE__);
  74                 exit(1);
  75         }
  76
  77         while (size--)
  78                 new[size] = old[size];
  79
  80         return new;
  81 }
  82
  83 //                                                                            //
  84
  85 void fprintVector(FILE *stream, Vector vector, int size)
  86 {
  87         int i;
  88
  89         assert(stream);
  90         assert(vector);
  91
  92         fprintf(stream, "%2d", vector[0]);
  93         for (i=1; i<size; i++)
  94                 fprintf(stream, " %2d", vector[i]);
  95 }
  96
  97 //                                                                            //
  98
  99 inline void printVector(Vector vector, int size)
 100 {
 101         assert(vector);
 102
 103         fprintVector(stdout, vector, size);
 104 }
 105
 106 //                                                                            //
 107
 108 Vector readVector(FILE *stream, int size)
 109 {
 110         int i;
 111         Vector vector;
 112
 113         assert(stream);
 114
 115         vector = createVector(size);
 116         assert(vector);
 117
 118         for (i=0; i<size; i++)
 119         {
 120                 if (fscanf(stream, "%d", &(vector[i]))<1)
 121                 {
 122                         deleteVector(vector);
 123                         return NULL;
 124                 }
 125         }
 126
 127         return vector;
 128 }
 129
 130 //                                                                            //
 131
 132 int normVector(Vector vector, int size)
 133 {
 134         int norm = 0;
 135
 136         assert(vector);
 137
 138         if (size>=0)
 139                 while (size--)
 140                         norm += abs(vector[size]);
 141
 142         return norm;
 143 }
 144
 145 //                                                                            //
 146
 147 void swapVector(Vector vector, int a, int b)
 148 {
 149         int tmp;
 150
 151         assert(vector);
 152
 153         tmp = vector[a];
 154         vector[a] = vector[b];
 155         vector[b] = tmp;
 156 }
 157
 158 //                                                                            //
 159
 160 int gcdVector(Vector vector, int size)
 161 {
 162         int i,result;
 163
 164         if (size==0)
 165                 return 0;
 166
 167         result = vector[0];
 168         for (i=1; i<size; i++)
 169                 result = gcd(result, vector[i]);
 170
 171         return result;
 172 }
 173
 174 //                                                                            //
 175
 176 int lexCompareInverseVector(Vector vector, int size) {
 177         int i = 0;
 178         while (i<size && vector[i]==0)
 179                 i++;
 180         return i==size ? 0 : vector[i];
 181 }
 182
 183 //                                                                            //