casa-mp-base
Public
Source
+static int std_drop(int lfil, int i, double tolL, double tolU, double toldiag) ;--- SRC/MatOps.c.orig 2013-10-08 23:09:05.000000000 -0400+++ SRC/MatOps.c 2013-10-08 23:09:17.000000000 -0400-#include <stdio.h>-#include <stdlib.h>-#include <string.h>-#include <math.h> -#include "globheads.h" #include "protos.h"- int diag_scal( vbsptr vbmat ){ /*---------------------------------------------------------------------------- * Diagonal scaling:--- SRC/PQ.c.orig 2013-10-08 23:09:05.000000000 -0400+++ SRC/PQ.c 2013-10-08 23:09:17.000000000 -0400-#include <stdio.h>-#include <stdlib.h>-#include <string.h>-#include <math.h>-#include "globheads.h" #include "protos.h" #define ALPHA 0.00001 |-----end-of-indsetPQ-------------------------------------------------- |--------------------------------------------------------------------*/-int add2is(int *last, int nod, int *iord, int *riord)-{-/*-----------------------------------------------------------------------| adds element nod to independent set-|---------------------------------------------------------------------*/- (*last)++;- iord[nod] = *last;- riord[*last] = nod;- return 0;-}-/*----------------------------------------------------------------------|---- end of add2is ----------------------------------------------------|--------------------------------------------------------------------*/-int add2com(int *nback, int nod, int *iord, int *riord) -{-/*-----------------------------------------------------------------------| adds element nod to independent set-|---------------------------------------------------------------------*/- iord[nod] = *nback;- riord[*nback] = nod;- (*nback)--;- return 0;-}-/*----------------------------------------------------------------------|---- end of add2com ---------------------------------------------------|--------------------------------------------------------------------*/-int indsetC(csptr mat, int bsize, int *iord, int *nnod, double tol) -{-/*--------------------------------------------------------------------- -| greedy algorithm for independent set ordering -- -|-----------------------------------------------------------------------| Input parameters:-| ------------------| (mat) = matrix in SpaFmt format