Hi, Philip,
It looks the performance of MatPtAP is pretty bad. There are a lot of
issues with PtAP, which I am going to address.
MatPtAPNumeric 181 1.0 nan nan 0.00e+00 0.0 3.3e+03 1.8e+04 0.0e+00
56 0 4 21 0 56 0 4 21 0 -nan -nan 0 0.00e+000 0.00e+00
0
Thanks.
The following is the log_view output for the ported case using 4 MPI tasks.
***WIDEN YOUR WINDOW TO 160
Hi, Philip,
Could you add -log_view and see what functions are used in the solve? Since
it is CPU-only, perhaps with -log_view of different runs, we can easily see
which functions slowed down.
--Junchao Zhang
From: Fackler, Philip
Sent: Tuesday, January 17,
In Xolotl's feature-petsc-kokkos branch I have ported the code to use petsc's
COO interface for creating the Jacobian matrix (and the Kokkos interface for
interacting with Vec entries). As the attached plots show for one case, while
the code for computing the RHSFunction and RHSJacobian perform
