dglprecsol.c

a software code for computing selected eigenvalues of large sparse symmetric matrices

#include <string.h>
#include <ilupack.h>
#include <ilupackmacros.h>

#define MAX(A,B)   (((A)>=(B))?(A):(B))
#define MIN(A,B)   (((A)<=(B))?(A):(B))

void  dglprecsol(size_t    *IPdiag, 
		 integer   *n,
		 FLOAT     *target,
		 FLOAT     *source) {

  FLOAT *diag;
  integer i;

  memcpy(&diag, IPdiag, sizeof(size_t));

  for (i=0; i<*n; i++) 
      target[i]=diag[i]*source[i];
}