paso/profiling/Test.c

#include <stdio.h>
#include "paso/Common.h"
#include "paso/Solver.h"
#include "paso/SystemMatrix.h"
#include "Paso_tests.h"

int main (int argc, char *argv[]) {
    Paso_SystemMatrix *A = NULL;
    double *b,*x;
    dim_t i,n,level=1;
   
    if (argc<2) {
        fprintf(stderr,"Please enter the filename\n");
        return -1;    
    }

    if (argc==3) {
      level=atoi(argv[2]);
    }
    
    A=MEMALLOC(1,Paso_SystemMatrix);
   
    A=Paso_SystemMatrix_loadMM_toCSR(argv[1]);
    if (A==NULL) {
      fprintf(stderr,"CSR Matrix not Loaded\n");
      return 0;
    }
    n=Paso_SystemMatrix_getTotalNumRows(A);
    b=MEMALLOC(n,double);
    x=MEMALLOC(n,double);
    for(i=0;i<n;i++)
    {
     b[i]=1;
    }
    
    Paso_test_run(A,b,level);
    MEMFREE(b);
    MEMFREE(x);
    Paso_SystemMatrix_free(A);
return 1;
}                                                                                                                                                                                                     
1	artak	2160	#include <stdio.h>
2			#include "paso/Common.h"
3			#include "paso/Solver.h"
4			#include "paso/SystemMatrix.h"
5	artak	2163	#include "Paso_tests.h"
6	artak	2160
7			int main (int argc, char *argv[]) {
8			Paso_SystemMatrix *A = NULL;
9	artak	2162	double b,x;
10	artak	2165	dim_t i,n,level=1;
11	artak	2160
12			if (argc<2) {
13			fprintf(stderr,"Please enter the filename\n");
14			return -1;
15			}
16	artak	2165
17			if (argc==3) {
18			level=atoi(argv[2]);
19			}
20	artak	2160
21			A=MEMALLOC(1,Paso_SystemMatrix);
22
23			A=Paso_SystemMatrix_loadMM_toCSR(argv[1]);
24			if (A==NULL) {
25	artak	2165	fprintf(stderr,"CSR Matrix not Loaded\n");
26	artak	2160	return 0;
27			}
28			n=Paso_SystemMatrix_getTotalNumRows(A);
29			b=MEMALLOC(n,double);
30			x=MEMALLOC(n,double);
31			for(i=0;i<n;i++)
32			{
33			b[i]=1;
34			}
35
36	artak	2165	Paso_test_run(A,b,level);
37	artak	2160	MEMFREE(b);
38			MEMFREE(x);
39			Paso_SystemMatrix_free(A);
40			return 1;
41			}