Removed drivers from ParaSails
This commit is contained in:
Rob Falgout
parent
120dbe0f9f
commit
040df7cdd3
@ -1 +0,0 @@
|
||||
insure++.file_ignore Mem.c
|
||||
@ -1,250 +0,0 @@
|
||||
/*BHEADER**********************************************************************
|
||||
* Copyright (c) 2008, Lawrence Livermore National Security, LLC.
|
||||
* Produced at the Lawrence Livermore National Laboratory.
|
||||
* This file is part of HYPRE. See file COPYRIGHT for details.
|
||||
*
|
||||
* HYPRE is free software; you can redistribute it and/or modify it under the
|
||||
* terms of the GNU Lesser General Public License (as published by the Free
|
||||
* Software Foundation) version 2.1 dated February 1999.
|
||||
*
|
||||
* $Revision$
|
||||
***********************************************************************EHEADER*/
|
||||
|
||||
|
||||
|
||||
|
||||
#include <stdlib.h>
|
||||
#include <assert.h>
|
||||
#include "Common.h"
|
||||
#include "Matrix.h"
|
||||
#include "ParaSails.h"
|
||||
#include "ConjGrad.h"
|
||||
|
||||
/*
|
||||
* Usage: driver symmetric num_runs matrixfile [rhsfile]
|
||||
*
|
||||
* If num_runs == 1, then hard-coded parameters will be used; else the
|
||||
* user will be prompted for parameters (except on the final run).
|
||||
*
|
||||
* To simulate diagonal preconditioning, use a large value of thresh,
|
||||
* e.g., thresh > 10.
|
||||
*/
|
||||
|
||||
HYPRE_Int main(HYPRE_Int argc, char *argv[])
|
||||
{
|
||||
HYPRE_Int mype, npes;
|
||||
HYPRE_Int symmetric;
|
||||
HYPRE_Int num_runs;
|
||||
Matrix *A;
|
||||
ParaSails *ps;
|
||||
FILE *file;
|
||||
HYPRE_Int n, beg_row, end_row;
|
||||
HYPRE_Real time0, time1;
|
||||
HYPRE_Real setup_time, solve_time;
|
||||
HYPRE_Real max_setup_time, max_solve_time;
|
||||
HYPRE_Real cost;
|
||||
|
||||
HYPRE_Real *x, *b;
|
||||
HYPRE_Int i, niter;
|
||||
HYPRE_Real thresh;
|
||||
HYPRE_Real threshg;
|
||||
HYPRE_Int nlevels;
|
||||
HYPRE_Real filter;
|
||||
HYPRE_Real loadbal;
|
||||
|
||||
hypre_MPI_Init(&argc, &argv);
|
||||
hypre_MPI_Comm_rank(hypre_MPI_COMM_WORLD, &mype);
|
||||
hypre_MPI_Comm_size(hypre_MPI_COMM_WORLD, &npes);
|
||||
|
||||
/* Read number of rows in matrix */
|
||||
symmetric = atoi(argv[1]);
|
||||
num_runs = atoi(argv[2]);
|
||||
|
||||
file = fopen(argv[3], "r");
|
||||
assert(file != NULL);
|
||||
#ifdef EMSOLVE
|
||||
hypre_fscanf(file, "%*d %d\n", &n);
|
||||
#else
|
||||
hypre_fscanf(file, "%d\n", &n);
|
||||
#endif
|
||||
fclose(file);
|
||||
assert(n >= npes);
|
||||
|
||||
beg_row = (HYPRE_Int) ((HYPRE_Real)(mype*n) / npes) + 1; /* assumes 1-based */
|
||||
end_row = (HYPRE_Int) ((HYPRE_Real)((mype+1)* n) / npes);
|
||||
|
||||
if (mype == 0)
|
||||
assert(beg_row == 1);
|
||||
if (mype == npes-1)
|
||||
assert(end_row == n);
|
||||
|
||||
#ifdef EMSOLVE
|
||||
beg_row--;
|
||||
end_row--;
|
||||
#endif
|
||||
|
||||
x = hypre_TAlloc(HYPRE_Real, (end_row-beg_row+1) , HYPRE_MEMORY_HOST);
|
||||
b = hypre_TAlloc(HYPRE_Real, (end_row-beg_row+1) , HYPRE_MEMORY_HOST);
|
||||
|
||||
A = MatrixCreate(hypre_MPI_COMM_WORLD, beg_row, end_row);
|
||||
|
||||
MatrixRead(A, argv[3]);
|
||||
if (mype == 0)
|
||||
hypre_printf("%s\n", argv[3]);
|
||||
|
||||
/* MatrixPrint(A, "A"); */
|
||||
|
||||
/* Right-hand side */
|
||||
if (argc > 4)
|
||||
{
|
||||
RhsRead(b, A, argv[4]);
|
||||
if (mype == 0)
|
||||
hypre_printf("Using rhs from %s\n", argv[4]);
|
||||
}
|
||||
else
|
||||
{
|
||||
for (i=0; i<end_row-beg_row+1; i++)
|
||||
b[i] = (HYPRE_Real) (2*rand()) / (HYPRE_Real) RAND_MAX - 1.0;
|
||||
}
|
||||
|
||||
while (num_runs && num_runs >= -1)
|
||||
{
|
||||
/* Initial guess */
|
||||
for (i=0; i<end_row-beg_row+1; i++)
|
||||
x[i] = 0.0;
|
||||
|
||||
if (num_runs == -1)
|
||||
{
|
||||
thresh = 0.0;
|
||||
nlevels = 0;
|
||||
filter = 0.0;
|
||||
loadbal = 0.0;
|
||||
}
|
||||
else
|
||||
{
|
||||
if (mype == 0)
|
||||
{
|
||||
#if PARASAILS_EXT_PATTERN
|
||||
hypre_printf("Enter parameters threshg, thresh, nlevels, "
|
||||
"filter, beta:\n");
|
||||
fflush(stdout);
|
||||
hypre_scanf("%lf %lf %d %lf %lf", &threshg, &thresh, &nlevels,
|
||||
&filter, &loadbal);
|
||||
#else
|
||||
hypre_printf("Enter parameters thresh, nlevels, "
|
||||
"filter, beta:\n");
|
||||
fflush(stdout);
|
||||
hypre_scanf("%lf %d %lf %lf", &thresh, &nlevels,
|
||||
&filter, &loadbal);
|
||||
#endif
|
||||
}
|
||||
|
||||
hypre_MPI_Bcast(&threshg, 1, hypre_MPI_REAL, 0, hypre_MPI_COMM_WORLD);
|
||||
hypre_MPI_Bcast(&thresh, 1, hypre_MPI_REAL, 0, hypre_MPI_COMM_WORLD);
|
||||
hypre_MPI_Bcast(&nlevels, 1, HYPRE_MPI_INT, 0, hypre_MPI_COMM_WORLD);
|
||||
hypre_MPI_Bcast(&filter, 1, hypre_MPI_REAL, 0, hypre_MPI_COMM_WORLD);
|
||||
hypre_MPI_Bcast(&loadbal, 1, hypre_MPI_REAL, 0, hypre_MPI_COMM_WORLD);
|
||||
|
||||
if (nlevels < 0)
|
||||
break;
|
||||
}
|
||||
|
||||
/**************
|
||||
* Setup phase
|
||||
**************/
|
||||
|
||||
hypre_MPI_Barrier(hypre_MPI_COMM_WORLD);
|
||||
time0 = hypre_MPI_Wtime();
|
||||
|
||||
ps = ParaSailsCreate(hypre_MPI_COMM_WORLD, beg_row, end_row, symmetric);
|
||||
|
||||
ps->loadbal_beta = loadbal;
|
||||
|
||||
#if PARASAILS_EXT_PATTERN
|
||||
ParaSailsSetupPatternExt(ps, A, threshg, thresh, nlevels);
|
||||
#else
|
||||
ParaSailsSetupPattern(ps, A, thresh, nlevels);
|
||||
#endif
|
||||
|
||||
time1 = hypre_MPI_Wtime();
|
||||
setup_time = time1-time0;
|
||||
|
||||
cost = ParaSailsStatsPattern(ps, A);
|
||||
if (cost > 5.e11)
|
||||
{
|
||||
hypre_printf("Aborting setup and solve due to high cost.\n");
|
||||
goto cleanup;
|
||||
}
|
||||
|
||||
hypre_MPI_Barrier(hypre_MPI_COMM_WORLD);
|
||||
time0 = hypre_MPI_Wtime();
|
||||
|
||||
err = ParaSailsSetupValues(ps, A, filter);
|
||||
if (err != 0)
|
||||
{
|
||||
hypre_printf("ParaSailsSetupValues returned error.\n");
|
||||
goto cleanup;
|
||||
}
|
||||
|
||||
time1 = hypre_MPI_Wtime();
|
||||
setup_time += (time1-time0);
|
||||
|
||||
ParaSailsStatsValues(ps, A);
|
||||
|
||||
if (!strncmp(argv[3], "testpsmat", 8))
|
||||
MatrixPrint(ps->M, "M");
|
||||
|
||||
#if 0
|
||||
if (mype == 0)
|
||||
hypre_printf("SETTING UP VALUES AGAIN WITH FILTERED PATTERN\n");
|
||||
ps->loadbal_beta = 0;
|
||||
ParaSailsSetupValues(ps, A, 0.0);
|
||||
#endif
|
||||
|
||||
/*****************
|
||||
* Solution phase
|
||||
*****************/
|
||||
|
||||
niter = 3000;
|
||||
if (MatrixNnz(ps->M) == n) /* if diagonal preconditioner */
|
||||
niter = 5000;
|
||||
|
||||
hypre_MPI_Barrier(hypre_MPI_COMM_WORLD);
|
||||
time0 = hypre_MPI_Wtime();
|
||||
|
||||
if (symmetric == 1)
|
||||
PCG_ParaSails(A, ps, b, x, 1.e-8, niter);
|
||||
else
|
||||
FGMRES_ParaSails(A, ps, b, x, 50, 1.e-8, niter);
|
||||
|
||||
time1 = hypre_MPI_Wtime();
|
||||
solve_time = time1-time0;
|
||||
|
||||
hypre_MPI_Reduce(&setup_time, &max_setup_time, 1, hypre_MPI_REAL, hypre_MPI_MAX, 0,
|
||||
hypre_MPI_COMM_WORLD);
|
||||
hypre_MPI_Reduce(&solve_time, &max_solve_time, 1, hypre_MPI_REAL, hypre_MPI_MAX, 0,
|
||||
hypre_MPI_COMM_WORLD);
|
||||
|
||||
if (mype == 0)
|
||||
{
|
||||
hypre_printf("**********************************************\n");
|
||||
hypre_printf("*** Setup Solve Total\n");
|
||||
hypre_printf("III %8.1f %8.1f %8.1f\n", max_setup_time, max_solve_time,
|
||||
max_setup_time+max_solve_time);
|
||||
hypre_printf("**********************************************\n");
|
||||
}
|
||||
|
||||
cleanup:
|
||||
ParaSailsDestroy(ps);
|
||||
|
||||
num_runs--;
|
||||
}
|
||||
|
||||
free(x);
|
||||
free(b);
|
||||
|
||||
MatrixDestroy(A);
|
||||
hypre_MPI_Finalize();
|
||||
|
||||
return 0;
|
||||
}
|
||||
@ -1,336 +0,0 @@
|
||||
/*BHEADER**********************************************************************
|
||||
* Copyright (c) 2008, Lawrence Livermore National Security, LLC.
|
||||
* Produced at the Lawrence Livermore National Laboratory.
|
||||
* This file is part of HYPRE. See file COPYRIGHT for details.
|
||||
*
|
||||
* HYPRE is free software; you can redistribute it and/or modify it under the
|
||||
* terms of the GNU Lesser General Public License (as published by the Free
|
||||
* Software Foundation) version 2.1 dated February 1999.
|
||||
*
|
||||
* $Revision$
|
||||
***********************************************************************EHEADER*/
|
||||
|
||||
|
||||
|
||||
|
||||
#include <stdlib.h>
|
||||
#include <assert.h>
|
||||
#include "Common.h"
|
||||
#include "Matrix.h"
|
||||
#include "ParaSails.h"
|
||||
#include "ConjGrad.h"
|
||||
|
||||
/* functions for instrumenting code */
|
||||
/* compile using
|
||||
cc -bI:/usr/lpp/pmtoolkit/lib/pmsvcs.exp -o counter counter.c -lpmapi */
|
||||
|
||||
#include <stdio.h>
|
||||
#if 0
|
||||
#include "pmapi.h"
|
||||
|
||||
#define ERROR_CODE -1
|
||||
#define OK_CODE 0
|
||||
|
||||
void hcinit();
|
||||
void hcreset(); /* will accumulate if not reset */
|
||||
void hcstart(); /* use start and stop to surround sections to instrument */
|
||||
void hcstop(); /* each stop will print a message */
|
||||
void hcfinish();
|
||||
|
||||
void hcinit()
|
||||
{
|
||||
pm_info_t myinfo;
|
||||
pm_prog_t prog;
|
||||
HYPRE_Int rc;
|
||||
HYPRE_Int filter = 0;
|
||||
|
||||
prog.mode.w = 0;
|
||||
|
||||
prog.events[0] = 15; /* floating point ops */
|
||||
prog.events[1] = 6; /* data cache misses */
|
||||
prog.events[2] = 1; /* cycles */
|
||||
prog.events[3] = 2; /* instr completed per cycle */
|
||||
|
||||
/* prog.mode.w |= PM_KERNEL; *//* count in kernel mode */
|
||||
prog.mode.w |= PM_USER; /* count in user mode */
|
||||
|
||||
prog.mode.w |= PM_PROCESS;
|
||||
|
||||
if ( (rc = pm_init(filter, &myinfo)) != OK_CODE) {
|
||||
pm_error("pm_init", rc);
|
||||
exit(ERROR_CODE);
|
||||
}
|
||||
|
||||
prog.mode.b.count = 1;
|
||||
|
||||
if ( (rc = pm_set_program_mygroup(&prog)) != OK_CODE) {
|
||||
pm_error("pm_set_program_mygroup", rc);
|
||||
exit(ERROR_CODE);
|
||||
}
|
||||
}
|
||||
|
||||
void hcreset()
|
||||
{
|
||||
HYPRE_Int rc;
|
||||
|
||||
if ( (rc = pm_reset_data_mygroup()) != OK_CODE)
|
||||
pm_error("pm_reset_data_mygroup", rc);
|
||||
}
|
||||
|
||||
void hcstart()
|
||||
{
|
||||
HYPRE_Int rc;
|
||||
|
||||
if ( (rc = pm_start_mygroup()) != OK_CODE)
|
||||
pm_error("pm_start_mygroup", rc);
|
||||
}
|
||||
|
||||
void hcstop()
|
||||
{
|
||||
pm_prog_t getprog;
|
||||
pm_data_t mydata;
|
||||
HYPRE_Int rc;
|
||||
|
||||
if ( (rc = pm_stop_mygroup()) != OK_CODE)
|
||||
pm_error("pm_stop_mygroup", rc);
|
||||
|
||||
if ( (rc = pm_get_program_mygroup(&getprog)) != OK_CODE)
|
||||
pm_error("pm_get_program_mygroup", rc);
|
||||
|
||||
if ( (rc = pm_get_data_mygroup(&mydata)) != OK_CODE)
|
||||
pm_error("pm_get_data_mygroup", rc);
|
||||
|
||||
hypre_printf("--------------\n");
|
||||
hypre_printf("flops: %-8lld\n", mydata.accu[0]);
|
||||
hypre_printf("dcmiss: %-8lld\n", mydata.accu[1]);
|
||||
hypre_printf("cycles: %-8lld\n", mydata.accu[2]);
|
||||
hypre_printf("instr: %-8lld\n", mydata.accu[3]);
|
||||
}
|
||||
|
||||
void hcfinish()
|
||||
{
|
||||
HYPRE_Int rc;
|
||||
|
||||
if ( (rc = pm_delete_program_mygroup()) != OK_CODE)
|
||||
pm_error("pm_delete_program_mygroup", rc);
|
||||
}
|
||||
#endif
|
||||
|
||||
HYPRE_Int rownum(const HYPRE_Int x, const HYPRE_Int y, const HYPRE_Int z,
|
||||
const HYPRE_Int nx, const HYPRE_Int ny, const HYPRE_Int nz, HYPRE_Int P, HYPRE_Int Q)
|
||||
{
|
||||
HYPRE_Int p, q, r;
|
||||
HYPRE_Int lowerx, lowery, lowerz;
|
||||
HYPRE_Int id, startrow;
|
||||
|
||||
p = (x-1) / nx;
|
||||
q = (y-1) / ny;
|
||||
r = (z-1) / nz;
|
||||
id = r*P*Q+q*P+p;
|
||||
startrow = id*(nx*ny*nz) + 1;
|
||||
lowerx = nx*p + 1;
|
||||
lowery = ny*q + 1;
|
||||
lowerz = nz*r + 1;
|
||||
|
||||
return startrow + nx*ny*(z-lowerz) + nx*(y-lowery) + (x-lowerx);
|
||||
}
|
||||
|
||||
HYPRE_Int main(HYPRE_Int argc, char *argv[])
|
||||
{
|
||||
HYPRE_Int npes, mype;
|
||||
|
||||
HYPRE_Int nx, ny, nz;
|
||||
HYPRE_Int P, Q, R;
|
||||
HYPRE_Real dx, dy, dz;
|
||||
HYPRE_Int p, q, r;
|
||||
HYPRE_Int lowerx, lowery, lowerz;
|
||||
HYPRE_Int upperx, uppery, upperz;
|
||||
HYPRE_Int x, y, z;
|
||||
HYPRE_Int num_rows;
|
||||
HYPRE_Int row;
|
||||
HYPRE_Int inds[100], *inds_p;
|
||||
HYPRE_Real coefs[100], *coefs_p;
|
||||
HYPRE_Int beg_row, end_row;
|
||||
HYPRE_Real time0, time1;
|
||||
HYPRE_Real setup_time, solve_time;
|
||||
HYPRE_Real max_setup_time, max_solve_time;
|
||||
|
||||
HYPRE_Real *x0, *b;
|
||||
HYPRE_Int i;
|
||||
Matrix *A;
|
||||
ParaSails *ps;
|
||||
|
||||
hypre_MPI_Init(&argc, &argv);
|
||||
|
||||
hypre_MPI_Comm_rank(hypre_MPI_COMM_WORLD, &mype);
|
||||
hypre_MPI_Comm_size(hypre_MPI_COMM_WORLD, &npes);
|
||||
|
||||
if (argc > 9)
|
||||
{
|
||||
nx = atoi(argv[1]);
|
||||
ny = atoi(argv[2]);
|
||||
nz = atoi(argv[3]);
|
||||
P = atoi(argv[4]);
|
||||
Q = atoi(argv[5]);
|
||||
R = atoi(argv[6]);
|
||||
dx = atof(argv[7]);
|
||||
dy = atof(argv[8]);
|
||||
dz = atof(argv[9]);
|
||||
}
|
||||
else
|
||||
{
|
||||
hypre_printf("Usage: mpirun -np %d %s <nx,ny,nz,P,Q,R,dx,dy,dz> ,\n\n",
|
||||
npes, argv[0]);
|
||||
hypre_printf(" where nx X ny X nz is the problem size per processor;\n");
|
||||
hypre_printf(" P X Q X R is the processor topology;\n");
|
||||
hypre_printf(" dx, dy, dz are the diffusion coefficients.\n");
|
||||
|
||||
exit(1);
|
||||
}
|
||||
|
||||
assert(npes == P*Q*R);
|
||||
|
||||
hypre_printf("XX side: %4d npes: %4d n: %10d\n", nx, npes, npes*nx*ny*nz);
|
||||
|
||||
/* compute p,q,r from P,Q,R and mype */
|
||||
p = mype % P;
|
||||
q = (( mype - p)/P) % Q;
|
||||
r = ( mype - p - P*q)/( P*Q );
|
||||
|
||||
/* compute ilower and iupper from p,q,r and nx,ny,nz */
|
||||
/* 1-based */
|
||||
|
||||
lowerx = nx*p + 1;
|
||||
upperx = nx*(p+1);
|
||||
lowery = ny*q + 1;
|
||||
uppery = ny*(q+1);
|
||||
lowerz = nz*r + 1;
|
||||
upperz = nz*(r+1);
|
||||
|
||||
num_rows = P*Q*R*nx*ny*nz;
|
||||
beg_row = mype*(nx*ny*nz)+1;
|
||||
end_row = (mype+1)*(nx*ny*nz);
|
||||
|
||||
x0 = hypre_TAlloc(HYPRE_Real, (end_row-beg_row+1) , HYPRE_MEMORY_HOST);
|
||||
b = hypre_TAlloc(HYPRE_Real, (end_row-beg_row+1) , HYPRE_MEMORY_HOST);
|
||||
|
||||
for (i=0; i<end_row-beg_row+1; i++)
|
||||
{
|
||||
b[i] = (HYPRE_Real) (2*rand()) / (HYPRE_Real) RAND_MAX - 1.0;
|
||||
x0[i] = 0.0;
|
||||
}
|
||||
|
||||
|
||||
A = MatrixCreate(hypre_MPI_COMM_WORLD, beg_row, end_row);
|
||||
|
||||
|
||||
for (z=lowerz; z<=upperz; z++)
|
||||
for (y=lowery; y<=uppery; y++)
|
||||
for (x=lowerx; x<=upperx; x++)
|
||||
{
|
||||
HYPRE_Int temp;
|
||||
|
||||
coefs_p = coefs;
|
||||
inds_p = inds;
|
||||
row = rownum(x,y,z,nx,ny,nz,P,Q);
|
||||
|
||||
*coefs_p++ = 2.0*(dx+dy+dz);
|
||||
*inds_p++ = row;
|
||||
if (x != 1)
|
||||
{*coefs_p++ = -dx; *inds_p++ = rownum(x-1,y,z,nx,ny,nz,P,Q);}
|
||||
if (x != P*nx)
|
||||
{*coefs_p++ = -dx; *inds_p++ = rownum(x+1,y,z,nx,ny,nz,P,Q);}
|
||||
if (y != 1)
|
||||
{*coefs_p++ = -dy; *inds_p++ = rownum(x,y-1,z,nx,ny,nz,P,Q);}
|
||||
if (y != Q*ny)
|
||||
{*coefs_p++ = -dy; *inds_p++ = rownum(x,y+1,z,nx,ny,nz,P,Q);}
|
||||
if (z != 1)
|
||||
{*coefs_p++ = -dz; *inds_p++ = rownum(x,y,z-1,nx,ny,nz,P,Q);}
|
||||
if (z != R*nz)
|
||||
{*coefs_p++ = -dz; *inds_p++ = rownum(x,y,z+1,nx,ny,nz,P,Q);}
|
||||
|
||||
temp = inds_p-inds;
|
||||
MatrixSetRow(A, row, temp, inds, coefs);
|
||||
}
|
||||
|
||||
MatrixComplete(A);
|
||||
|
||||
/**************
|
||||
* Setup phase
|
||||
**************/
|
||||
|
||||
/*
|
||||
hcinit();
|
||||
hcreset();
|
||||
hcstart();
|
||||
*/
|
||||
hypre_MPI_Barrier(hypre_MPI_COMM_WORLD);
|
||||
time0 = hypre_MPI_Wtime();
|
||||
|
||||
ps = ParaSailsCreate(hypre_MPI_COMM_WORLD, beg_row, end_row, 1);
|
||||
|
||||
ParaSailsSetupPattern(ps, A, .1, 3);
|
||||
ParaSailsSetupValues(ps, A, 0.00);
|
||||
|
||||
time1 = hypre_MPI_Wtime();
|
||||
/*
|
||||
hcstop();
|
||||
*/
|
||||
setup_time = time1-time0;
|
||||
fflush(stdout);
|
||||
hypre_MPI_Barrier(hypre_MPI_COMM_WORLD);
|
||||
|
||||
/*****************
|
||||
* Solution phase
|
||||
*****************/
|
||||
|
||||
/*
|
||||
hcreset();
|
||||
hcstart();
|
||||
*/
|
||||
hypre_MPI_Barrier(hypre_MPI_COMM_WORLD);
|
||||
time0 = hypre_MPI_Wtime();
|
||||
|
||||
/*
|
||||
PCG_ParaSails(A, ps, b, x0, 1.e-8, 1500);
|
||||
*/
|
||||
|
||||
time1 = hypre_MPI_Wtime();
|
||||
/*
|
||||
hcstop();
|
||||
*/
|
||||
solve_time = time1-time0;
|
||||
|
||||
ParaSailsStatsPattern(ps, A);
|
||||
ParaSailsStatsValues(ps, A);
|
||||
|
||||
hypre_MPI_Reduce(&setup_time, &max_setup_time, 1, hypre_MPI_REAL, hypre_MPI_MAX, 0,
|
||||
hypre_MPI_COMM_WORLD);
|
||||
hypre_MPI_Reduce(&solve_time, &max_solve_time, 1, hypre_MPI_REAL, hypre_MPI_MAX, 0,
|
||||
hypre_MPI_COMM_WORLD);
|
||||
|
||||
if (mype == 0)
|
||||
{
|
||||
hypre_printf("**********************************************\n");
|
||||
hypre_printf("*** Setup Solve Total\n");
|
||||
hypre_printf("III %8.1f %8.1f %8.1f\n", max_setup_time, max_solve_time,
|
||||
max_setup_time+max_solve_time);
|
||||
hypre_printf("**********************************************\n");
|
||||
}
|
||||
|
||||
/*
|
||||
hcfinish();
|
||||
*/
|
||||
|
||||
ParaSailsDestroy(ps);
|
||||
/*
|
||||
if (mype == 6)
|
||||
HashPrint(A->numb->hash);
|
||||
*/
|
||||
MatrixDestroy(A);
|
||||
|
||||
free(x0);
|
||||
free(b);
|
||||
|
||||
hypre_MPI_Finalize();
|
||||
}
|
||||
Loading…
Reference in New Issue
Block a user