void apbsfinal_() {
unsigned long int bytesTotal, highWater;
int i;
VASSERT(nosh != VNULL);
/* Kill the saved potential Vgrids */
for (i=0; i<2; i++){
Vgrid_dtor(&permU[i]);
Vgrid_dtor(&indU[i]);
Vgrid_dtor(&nlIndU[i]);
}
Valist_dtor(&alist[0]);
/* Saving the kappa and dielectric maps is under consideration.
killKappaMaps(nosh, kappaMap);
killDielMaps(nosh, dielXMap, dielYMap, dielZMap);
*/
NOsh_dtor(&nosh);
/* Clean up MALOC structures */
bytesTotal = Vmem_bytesTotal();
highWater = Vmem_highWaterTotal();
/*
printf(" Final APBS memory usage: %4.3f MB total, %4.3f MB high water\n\n",
(double)(bytesTotal)/(1024.*1024.),(double)(highWater)/(1024.*1024.));
*/
Vmem_dtor(&mem);
Vnm_tstop(APBS_TIMER_WALL_CLOCK, "APBS WALL CLOCK");
Vcom_finalize();
Vcom_dtor(&com);
}
/*
* ***************************************************************************
* Routine: Vmp_finalize
*
* Purpose: The Vmp finalizer.
*
* Author: Michael Holst
* ***************************************************************************
*/
VPUBLIC int Vmp_finalize(void)
{
# if defined(USE_VCOM)
return Vcom_finalize();
# else
return Vmpi_finalize();
# endif
}
int main(int argc, char *argv[]) {
int i, testi, testj;
int myrank;
int nproc;
char string[100];
Vcom *vcom;
VASSERT( Vcom_init(&argc,&argv) );
vcom = Vcom_ctor(1);
myrank = Vcom_rank(vcom);
nproc = Vcom_size(vcom);
printf("PE %d: Starting test program with %d total PEs.\n",myrank,nproc);
/* Character send/recv test */
if (myrank == 0) {
printf("\n\nPE %d: STARTING SEND/RECV TEST.\n",myrank);
printf("PE %d: Sending non-blocked messages.\n",myrank);
for (i=1; i<nproc; i++) {
sprintf(&(string[0]), "non-blocked message (PE %d --> %d)",
myrank, i);
printf( "PE %d: Sent non-blocked message to PE %d\n", myrank, i);
assert(Vcom_send(vcom,i,string,100,3,0));
}
printf("PE %d: Sending blocked messages.\n",myrank);
for (i=1; i<nproc; i++) {
sprintf(&(string[0]),"blocked message (PE %d --> %d)",i,nproc);
printf( "PE %d: Sent blocked message to PE %d\n", myrank, i);
assert(Vcom_send(vcom,i,string,100,3,1));
}
} else {
printf("PE %d: Recving blocked message from PE 0.\n",myrank);
/* Get blocked sent messages */
assert(Vcom_recv(vcom,0,string,100,3,1));
printf( "PE %d: Blocked message is: '%s'\n",myrank,string);
/* Get blocked sent messages */
printf("PE %d: Recving non-blocked message from PE 0.\n",myrank);
assert(Vcom_recv(vcom,0,string,100,3,1));
printf( "PE %d: Non-blocked message is: '%s'\n",myrank,string);
}
fflush(stdout);
Vcom_barr(vcom);
/* Barrier test */
if (myrank == 0) printf("\n\nPE %d: STARTING BARRIER TEST.\n",myrank);
printf("PE %d: Hit barrier\n", myrank);
Vcom_barr(vcom);
printf("PE %d: Passed barrier\n", myrank);
fflush(stdout);
Vcom_barr(vcom);
/* getCount test */
if (myrank == 0) {
printf("\n\nPE %d: STARTING PROBE TEST.\n",myrank);
printf("PE %d: Sending string of length 100\n", myrank);
for (i=1; i<nproc; i++) assert(Vcom_send(vcom,i,string,100,3,0));
} else {
Vcom_getCount(vcom, 0, &testi, 3);
printf("PE %d: Probed for string of length %d\n", myrank, testi);
}
fflush(stdout);
Vcom_barr(vcom);
/* Reduction test */
testi = 4;
if (myrank == 0) {
printf("\n\nPE %d: STARTING REDUCTION TEST.\n",myrank);
printf("PE %d: All %d PEs have the number %d.\n", myrank, nproc, testi);
printf("PE %d: Performing global sum reduction.\n", myrank);
}
Vcom_reduce(vcom, &testi, &testj, 1, 1, 0);
printf("PE %d: Global sum = %d\n", myrank, testj);
fflush(stdout);
Vcom_barr(vcom);
testi = 4;
if (myrank == 0) {
testj = 8;
printf("PE %d: PE %d has the number %d; all others have %d.\n", myrank,
myrank, testj, testi);
testi = testj;
printf("PE %d: Performing global max reduction.\n", myrank);
}
Vcom_reduce(vcom, &testi, &testj, 1, 1, 3);
printf("PE %d: Global max = %d\n", myrank, testj);
fflush(stdout);
Vcom_barr(vcom);
/* Resize the communicator */
printf("PE %d: Resizing communicator from %d to %d.\n", myrank, nproc,
(int)(nproc/2));
Vcom_resize(vcom, (int)(nproc/2));
myrank = Vcom_rank(vcom);
if (Vcom_rank(vcom) != -1) printf("PE %d: Hello world from resized communicator\n", myrank);
Vcom_barr(vcom);
/* Finish up */
printf("PE %d: Exiting test program.\n",myrank);
fflush(stdout);
Vcom_dtor(&vcom);
VASSERT( Vcom_finalize() );
return 0;
}
