bpsraw.c (cb32e2e7f026784d97a57f1901677e9727def907) bpsraw.c (dc7d240c1c3ecbc40e893ea38e37c8e19d48593c)
1// Copyright (c) 2017, Lawrence Livermore National Security, LLC. Produced at
2// the Lawrence Livermore National Laboratory. LLNL-CODE-734707. All Rights
3// reserved. See files LICENSE and NOTICE for details.
4//
5// This file is part of CEED, a collection of benchmarks, miniapps, software
6// libraries and APIs for efficient high-order finite element and spectral
7// element discretizations for exascale applications. For more information and
8// source code availability see http://github.com/ceed.

--- 467 unchanged lines hidden (view full) ---

476 ierr = VecSetSizes(X, mnodes[0]*mnodes[1]*mnodes[2]*ncompu, PETSC_DECIDE);
477 CHKERRQ(ierr);
478 ierr = VecSetUp(X); CHKERRQ(ierr);
479
480 // Set up libCEED
481 CeedInit(ceedresource, &ceed);
482
483 // Print summary
1// Copyright (c) 2017, Lawrence Livermore National Security, LLC. Produced at
2// the Lawrence Livermore National Laboratory. LLNL-CODE-734707. All Rights
3// reserved. See files LICENSE and NOTICE for details.
4//
5// This file is part of CEED, a collection of benchmarks, miniapps, software
6// libraries and APIs for efficient high-order finite element and spectral
7// element discretizations for exascale applications. For more information and
8// source code availability see http://github.com/ceed.

--- 467 unchanged lines hidden (view full) ---

476 ierr = VecSetSizes(X, mnodes[0]*mnodes[1]*mnodes[2]*ncompu, PETSC_DECIDE);
477 CHKERRQ(ierr);
478 ierr = VecSetUp(X); CHKERRQ(ierr);
479
480 // Set up libCEED
481 CeedInit(ceedresource, &ceed);
482
483 // Print summary
484 CeedInt gsize;
485 ierr = VecGetSize(X, &gsize); CHKERRQ(ierr);
484 if (!test_mode) {
486 if (!test_mode) {
485 CeedInt gsize;
486 ierr = VecGetSize(X, &gsize); CHKERRQ(ierr);
487 const char *usedresource;
488 CeedGetResource(ceed, &usedresource);
489 ierr = PetscPrintf(comm,
490 "\n-- CEED Benchmark Problem %d -- libCEED + PETSc --\n"
491 " libCEED:\n"
492 " libCEED Backend : %s\n"
493 " Mesh:\n"
494 " Number of 1D Basis Nodes (p) : %d\n"

--- 309 unchanged lines hidden (view full) ---

804 KSPConvergedReason reason;
805 PetscReal rnorm;
806 PetscInt its;
807 ierr = KSPGetType(ksp, &ksptype); CHKERRQ(ierr);
808 ierr = KSPGetConvergedReason(ksp, &reason); CHKERRQ(ierr);
809 ierr = KSPGetIterationNumber(ksp, &its); CHKERRQ(ierr);
810 ierr = KSPGetResidualNorm(ksp, &rnorm); CHKERRQ(ierr);
811 if (!test_mode || reason < 0 || rnorm > 1e-8) {
487 const char *usedresource;
488 CeedGetResource(ceed, &usedresource);
489 ierr = PetscPrintf(comm,
490 "\n-- CEED Benchmark Problem %d -- libCEED + PETSc --\n"
491 " libCEED:\n"
492 " libCEED Backend : %s\n"
493 " Mesh:\n"
494 " Number of 1D Basis Nodes (p) : %d\n"

--- 309 unchanged lines hidden (view full) ---

804 KSPConvergedReason reason;
805 PetscReal rnorm;
806 PetscInt its;
807 ierr = KSPGetType(ksp, &ksptype); CHKERRQ(ierr);
808 ierr = KSPGetConvergedReason(ksp, &reason); CHKERRQ(ierr);
809 ierr = KSPGetIterationNumber(ksp, &its); CHKERRQ(ierr);
810 ierr = KSPGetResidualNorm(ksp, &rnorm); CHKERRQ(ierr);
811 if (!test_mode || reason < 0 || rnorm > 1e-8) {
812 ierr = MPI_Allreduce(&my_rt, &rt_min, 1, MPI_DOUBLE, MPI_MIN, comm);
813 CHKERRQ(ierr);
814 ierr = MPI_Allreduce(&my_rt, &rt_max, 1, MPI_DOUBLE, MPI_MAX, comm);
815 CHKERRQ(ierr);
812 ierr = PetscPrintf(comm,
813 " KSP:\n"
814 " KSP Type : %s\n"
815 " KSP Convergence : %s\n"
816 " Total KSP Iterations : %D\n"
817 " Final rnorm : %e\n",
818 ksptype, KSPConvergedReasons[reason], its,
819 (double)rnorm); CHKERRQ(ierr);
816 ierr = PetscPrintf(comm,
817 " KSP:\n"
818 " KSP Type : %s\n"
819 " KSP Convergence : %s\n"
820 " Total KSP Iterations : %D\n"
821 " Final rnorm : %e\n",
822 ksptype, KSPConvergedReasons[reason], its,
823 (double)rnorm); CHKERRQ(ierr);
824 ierr = PetscPrintf(comm,
825 " Performance:\n"
826 " CG Solve Time : %g (%g) sec\n"
827 " DoFs/Sec in CG : %g (%g) million\n",
828 rt_max, rt_min, 1e-6*gsize*its/rt_max,
829 1e-6*gsize*its/rt_min); CHKERRQ(ierr);
820 }
821 if (benchmark_mode && (!test_mode)) {
830 }
831 if (benchmark_mode && (!test_mode)) {
822 CeedInt gsize;
823 ierr = VecGetSize(X, &gsize); CHKERRQ(ierr);
824 ierr = MPI_Allreduce(&my_rt, &rt_min, 1, MPI_DOUBLE, MPI_MIN, comm);
825 CHKERRQ(ierr);
826 ierr = MPI_Allreduce(&my_rt, &rt_max, 1, MPI_DOUBLE, MPI_MAX, comm);
827 CHKERRQ(ierr);
828 ierr = PetscPrintf(comm,
829 " Performance:\n"
830 " CG Solve Time : %g (%g) sec\n"
831 " DoFs/Sec in CG : %g (%g) million\n",

--- 60 unchanged lines hidden --- 		832 ierr = MPI_Allreduce(&my_rt, &rt_min, 1, MPI_DOUBLE, MPI_MIN, comm);
833 CHKERRQ(ierr);
834 ierr = MPI_Allreduce(&my_rt, &rt_max, 1, MPI_DOUBLE, MPI_MAX, comm);
835 CHKERRQ(ierr);
836 ierr = PetscPrintf(comm,
837 " Performance:\n"
838 " CG Solve Time : %g (%g) sec\n"
839 " DoFs/Sec in CG : %g (%g) million\n",

--- 60 unchanged lines hidden ---