Commit 100356fb authored by Stefan's avatar Stefan
Browse files

cosmetic changes

parent 0bb97466
......@@ -56,7 +56,7 @@ void ellipticBuildContinuous(elliptic_t* elliptic,
mesh_t *mesh = elliptic->mesh;
MPI_Barrier(mesh->comm);
const double tStart = MPI_Wtime();
if(mesh->rank == 0) printf("Building full FEM matrix ... ");
if(mesh->rank == 0) printf("building full FEM matrix ... ");
fflush(stdout);
switch(elliptic->elementType) {
......
......@@ -131,6 +131,11 @@ void ellipticBuildContinuousGalerkinHex3D(elliptic_t* elliptic,
mesh_t* mesh = elliptic->mesh;
setupAide options = elliptic->options;
MPI_Barrier(mesh->comm);
const double tStart = MPI_Wtime();
if(mesh->rank == 0) printf("building full FEM matrix using Galerkin projection ... ");
fflush(stdout);
int rank = mesh->rank;
//use the masked gs handle to define a global ordering
......@@ -175,9 +180,6 @@ void ellipticBuildContinuousGalerkinHex3D(elliptic_t* elliptic,
int* mask = (int*) calloc(mesh->Np * mesh->Nelements,sizeof(int));
for (dlong n = 0; n < elliptic->Nmasked; n++) mask[elliptic->maskIds[n]] = 1;
if(mesh->rank == 0) printf("Building full FEM matrix via Galerkin projection...");
fflush(stdout);
mesh_t* meshf = ellipticFine->mesh;
dfloat* b,* q,* Aq;
......@@ -308,7 +310,8 @@ void ellipticBuildContinuousGalerkinHex3D(elliptic_t* elliptic,
if (*nnz) cnt++;
*nnz = cnt;
if(mesh->rank == 0) printf("done.\n");
MPI_Barrier(mesh->comm);
if(mesh->rank == 0) printf("done (%gs)\n", MPI_Wtime() - tStart);
MPI_Barrier(mesh->comm);
MPI_Type_free(&MPI_NONZERO_T);
......
......@@ -80,7 +80,7 @@ void ellipticBuildJacobi(elliptic_t* elliptic, dfloat** invDiagA)
MPI_Barrier(mesh->comm);
const double tStart = MPI_Wtime();
if(mesh->rank == 0) printf("Building Jacobi ... ");
if(mesh->rank == 0) printf("building Jacobi ... ");
fflush(stdout);
// surface mass matrices MS = MM*LIFT
......
......@@ -10,7 +10,7 @@ void meshNekReaderHex3D(int N, mesh_t* mesh)
{
MPI_Barrier(mesh->comm);
const double tStart = MPI_Wtime();
if(mesh->rank == 0) printf("loading mesh from nek ...\n"); fflush(stdout);
if(mesh->rank == 0) printf("loading mesh from nek ... "); fflush(stdout);
mesh->dim = 3;
mesh->Nverts = 8;
......@@ -60,8 +60,7 @@ void meshNekReaderHex3D(int N, mesh_t* mesh)
if(mesh->rank == 0) {
int n = nekData.NboundaryIDt;
if(!mesh->cht) n = nekData.NboundaryID;
printf("NboundaryIDs: %d\n", n);
printf("NboundaryFaces: %d\n", mesh->NboundaryFaces);
printf("NboundaryIDs: %d, NboundaryFaces: %d ", n, mesh->NboundaryFaces);
}
int cnt = 0;
......@@ -121,5 +120,5 @@ void meshNekReaderHex3D(int N, mesh_t* mesh)
}
MPI_Barrier(mesh->comm);
if(mesh->rank == 0) printf("done (%gs)\n", MPI_Wtime() - tStart);
if(mesh->rank == 0) printf("done (%gs)\n", MPI_Wtime() - tStart); fflush(stdout);
}
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment