permrv.c

http://www.math.tu-berlin.de/~bolle/SPARSPAK/index-e.html

#include <sparspak.h>

/* #include <exec/types.h>  */
/* #include <exec/memory.h> */
#include <malloc.h>

#define MC_EXTERNAL       0
#define Malloc(A,B)       malloc(A)
#define Free(A)           free(A)
#define Realloc(A,B,C)    realloc(A,B)


/*****************************************************************************/
/* 'permrv' maps rhs[1],...,rhs[n] onto rhs[perm[1]],...,rhs[perm[n]].       */
/* In C, 'rhs, perm' have SUBSCRIPTS 0,1,2,...,*n-1, while the components of */
/* 'perm' are between 1 and *n.                                              */
/*****************************************************************************/
void permrv(int *n, double *rhs, int *perm)
{
   int    i;
   double *dummy;
   
   dummy=(double *)Malloc(*n*sizeof(double), MC_EXTERNAL);

   for (i=0; i<*n; i++)
       dummy[i]=rhs[perm[i]-1];
   for (i=0; i<*n; i++)
       rhs[i]=dummy[i];

   Free(dummy);
} /* end permrv */