Pull-out test case elastodamage#
Setting#
For parallel testing#
Create a server controlling n>1 MPI engines (consider OpenMPI as MPI backend)
#nopy
import ipyparallel as ipp
nw=8
cluster = ipp.Cluster(engines="mpi", n=nw)
rc = cluster.start_and_connect_sync()
view=rc[:]
Import TwoScale library#
%%px
from twoscale import core
from twoscale import linear
from twoscale import util
Test case construction#
Use dirty constants used to generate meshes representing some dimensions (stub, …)
Load pre-refined mesh
Scale jump construction#
Covering function (thick cone)#
Covering analysis based on the following geometry. A conic envelope is defined by:
axis is y that passes by 0,0,0
2 cone equations:
origin 0,ho1,0 and angle a
origin 0,ho2,0 and angle a
2 xz planes passing by 0,h1,0 (delimit inferior part) and 0,h2,0 (upper part) where ho1, ho2, a, h1, h2 are constants pre-computed
%%px
grow=lc*0.79
a1=35.*np.pi/180.
#a1=55.*np.pi/180.
tga=np.tan(a1)
sina=np.sin(a1)
cosa=np.cos(a1)
cosa2=cosa*cosa
hconnect=-450.-0.8*grow/sina+3.
ho1= hconnect - 60. / tga - grow / sina
ho2= ho1 + 2. * grow / sina
h1= -469.
#
def coneCovering(coords,level):
r=np.sqrt(coords[0]*coords[0]+coords[2]*coords[2])
ho1l=ho1*(1+0/(level+1))
ho2l=ho2*(1-0/(level+1))
return np.logical_and(
np.logical_not(np.logical_or(coords[1]>h2,coords[1]<h1)),
np.logical_and(
(tga*(coords[1]-ho1l)-r)>0.,
(r-tga*(coords[1]-ho2l))>0.)
)
Choose a mesh refinement localization criterion#
Use the coneCovering function for that
%%px
def crit(coords,level):
return coneCovering(coords,level)
Choose macro node to enrich#
All nodes with support in the cone envelope
%%px
def enriched(coords):
return coneCovering(coords,0)
Compute scale jump object#
%%px
j=core.topDown(domain,enriched,crit,ref,False,False)
%%px
fdomain=j.getFineMesh
cdomain=j.getCoarseMesh
Fine scale discretisation#
Space#
%%px
space = fem.functionspace(fdomain, fdomain.ufl_domain().ufl_coordinate_element())
Lamè#
%%px
E=26400.00
nu=0.193
mu=fem.Constant(fdomain,E / (2.0 * (1.0 + nu)))
lmbda = fem.Constant(fdomain,E * nu / ((1.0 + nu) * (1.0 - 2.0 * nu)))
Damage#
%%px
sp=0.8/(grow*grow)
lp=0.2
def profile(d):
return (1.-d/grow+sp*d*(grow-d)*np.exp(-d*lp))
#nopy
grow=view['grow'][0]
sp=view['sp'][0]
lp=view['lp'][0]
def profile(d):
return 1.-d/grow+sp*d*(grow-d)*np.exp(-d*lp)
x = np.linspace(0, grow, 200)
prof=pv.Chart2D()
prof.line(x,profile(x), color='r')
prof.line(x,1.-x/grow)
prof.show()
%%px
def dam(coords):
r=np.sqrt(coords[0]*coords[0]+coords[2]*coords[2])
ho=(ho1+ho2)/2
test=np.logical_and(
np.logical_not(np.logical_or(coords[1]>h2,coords[1]<h1)),
np.logical_and(
(tga*(coords[1]-ho1)-r)>0.,
(r-tga*(coords[1]-ho2))>0.))
yc=(coords[1]-ho+r*tga)*cosa2+ho
return np.where(test,profile(np.sqrt(np.power(coords[1]-yc,2)+np.power(r-(yc-ho)*tga,2))),0)
%%px
element_scal = basix.ufl.element("Lagrange", "tetrahedron", 1)
space_scal = fem.functionspace(fdomain, element_scal)
space_scalc = fem.functionspace(cdomain, element_scal)
d=fem.Function(space_scal,name="d")
d.interpolate(dam)
dc=fem.Function(space_scalc,name="dc")
dc.interpolate(dam)
Formulation#
Bilinear form#
%%px
u = ufl.TrialFunction(space)
v = ufl.TestFunction(space)
%%px
def eps(v):
return 0.5*(ufl.grad(v) + ufl.grad(v).T)
%%px
def sigma(strain):
return (1-d)*(2.0*mu*strain + lmbda*ufl.tr(strain)*ufl.Identity(dim))
%%px
a_ufl=ufl.inner(sigma(eps(u)), eps(v)) * ufl.dx
Surfacic load#
Define ds to integrate surface loading
%%px
r2=(dh/2)**2
small=0.00001
reducer2=1.5
if deri_load:
r2=r2/reducer2
load_faces=mesh.locate_entities(fdomain, 2, lambda x: np.logical_and(np.logical_and(x[1]-px+hef<=small,x[1]-px+hef>=-small),
x[0]*x[0]+x[2]*x[2]-r2<=small))
%%px
facet_markers = np.full_like(load_faces, 100,dtype=np.int32)
sorted_facets = np.argsort(load_faces)
nfacets_tag = mesh.meshtags(
fdomain, 2, load_faces[sorted_facets], facet_markers
)
ds = ufl.Measure("ds", domain=fdomain, subdomain_data=nfacets_tag)
Define load in y direction on surface under stub head
%%px
if deri_load:
tload=0.
else:
tload=load
stuby=fem.Constant(fdomain,[0.,tload,0.])
b_ufl=ufl.inner(stuby,v)*ds(100)
%%px
if deri_load:
load_dofs=fem.locate_dofs_topological(space.sub(1),2,load_faces)
bc_load=fem.dirichletbc(imp,load_dofs,space.sub(1))
Dirichelet#
Support clamped in y direction
%%px
def locate_support(teta,dom):
cs = np.cos(teta)
si = np.sin(teta)
lower = dhef-a/2
up = dhef+a/2
c2 = c/2.
support_face=mesh.locate_entities(dom, 2, lambda x: np.logical_and(x[1]==py,
np.logical_and(np.logical_and(si*x[0]-cs*x[2]+c2>=-small,si*x[0]-cs*x[2]-c2<=small),
np.logical_and(cs*x[0]+si*x[2]-lower>=-small,cs*x[0]+si*x[2]-up<=small))))
return support_face
delta = np.pi/4.;
support_faces=np.empty(0,dtype=np.int32)
for k in range(8):
teta = k*delta+delta
support_faces=np.concatenate((support_faces,locate_support(teta,fdomain)))
%%px
support_dofs=fem.locate_dofs_topological(space.sub(1),2,support_faces)
bc_support=fem.dirichletbc(0.,support_dofs,space.sub(1))
stub cylindar clamped in x,z direction
%%px
r2=(d0/2)**2
#dd=hef+h2
dd=0
if deri_load:
r2=r2/reducer2
stub_faces=mesh.locate_entities(fdomain, 2, lambda x: np.logical_and(np.logical_and(x[1]<=small,x[1]-px+hef-dd>=-small),
x[0]*x[0]+x[2]*x[2]-r2<=small))
%%px
stub1_x_dofs=fem.locate_dofs_topological(space.sub(0),2,stub_faces)
bc_stub1_x=fem.dirichletbc(0.,stub1_x_dofs,space.sub(0))
stub1_z_dofs=fem.locate_dofs_topological(space.sub(2),2,stub_faces)
bc_stub1_z=fem.dirichletbc(0.,stub1_z_dofs,space.sub(2))
#stub2_x_dofs=fem.locate_dofs_topological(space.sub(0),2,facets_tag.find(400))
#bc_stub2_x=fem.dirichletbc(0.,stub1_x_dofs,space.sub(0))
#stub2_z_dofs=fem.locate_dofs_topological(space.sub(2),2,facets_tag.find(400))
#bc_stub2_z=fem.dirichletbc(0.,stub1_z_dofs,space.sub(2))
Dirichlet is all condition above
%%px
if deri_load:
bcs=[bc_support,bc_stub1_x,bc_stub1_z,bc_load]
else:
bcs=[bc_support,bc_stub1_x,bc_stub1_z]
#bcs=[bc_support,bc_stub1_x,bc_stub1_z,bc_stub2_x,bc_stub2_z]
mpc#
%%px
mpc=core.generateMPC(j,space)
Solve with direct resolution#
%%px
#problem = petsc.LinearProblem(a_ufl, b_ufl,bcs=bcs,petsc_options=petsc_options)
problem = LinearProblem_mpc(a_ufl, b_ufl,mpc,bcs=bcs,petsc_options=petsc_options)
%%px
stageF = petsc4py.PETSc.Log.Stage("solveFine")
stageF.push()
dt=common.Timer("Fine solve")
dt.start()
disp = problem.solve()
dt.stop()
stageF.pop()
del problem
%%px
disp.name="displacements"
print(plt.camera_position)
[(-33.0238996678173, -354.60971081475634, -386.6175118462205),
(49.764908472958275, -465.79163270141413, 466.73312098224534),
(0.002783776868024144, 0.9916497049175196, 0.12893065315643742)]
Two-scale approach#
Use nested strategy
Enriched space#
If the nested strategy is used, the enriched space is just the standard space. Otherwise, the enriched space is a mixed space.
%%px
if nested_strategy:
std=fem.functionspace(cdomain, cdomain.ufl_domain().ufl_coordinate_element())
else:
el_mixed = basix.ufl.mixed_element([cdomain.ufl_domain().ufl_coordinate_element(), cdomain.ufl_domain().ufl_coordinate_element()])
coarse_enriched_space = fem.functionspace(cdomain, el_mixed)
std_=coarse_enriched_space.sub(0)
std, std_to_mix=std_.collapse()
enr_=coarse_enriched_space.sub(1)
enr, enr_to_mix=enr_.collapse()
std_to_mix_np=np.array(std_to_mix[0],dtype=int)
disp_ce=fem.Function(coarse_enriched_space,name='disp_coarse_enriched')
Solve coarse non-enriched problem (standard part) to initialize two-scale loop#
Formulation#
%%px
muc=fem.Constant(cdomain,E / (2.0 * (1.0 + nu)))
lmbdac = fem.Constant(cdomain,E * nu / ((1.0 + nu) * (1.0 - 2.0 * nu)))
uc = ufl.TrialFunction(std)
vc = ufl.TestFunction(std)
def sigmac(strain):
return (1-dc)*(2.0*muc*strain + lmbdac*ufl.tr(strain)*ufl.Identity(dim))
ac_ufl=ufl.inner(sigmac(eps(uc)), eps(vc)) * ufl.dx
%%px
r2=(dh/2)**2
load_facesc=mesh.locate_entities(cdomain, 2, lambda x: np.logical_and(np.logical_and(x[1]-px+hef<=small,x[1]-px+hef>=-small),x[0]*x[0]+x[2]*x[2]-r2<=small))
facet_markersc = np.full_like(load_facesc, 100,dtype=np.int32)
sorted_facetsc = np.argsort(load_facesc)
nfacets_tagc = mesh.meshtags(
cdomain, 2, load_facesc[sorted_facetsc], facet_markersc
)
dsc = ufl.Measure("ds", domain=cdomain, subdomain_data=nfacets_tagc)
%%px
stubyc=fem.Constant(cdomain,[0.,tload,0.])
dsx = ufl.Measure("ds", domain=cdomain, subdomain_data=nfacets_tagc)
bc_ufl=ufl.inner(stubyc,vc)*dsc(100)
Dirichlet#
%%px
support_facesc=np.empty(0,dtype=np.int32)
for k in range(8):
teta = k*delta+delta
support_facesc=np.concatenate((support_facesc,locate_support(teta,cdomain)))
r2=(d0/2)**2
stub_facesc=mesh.locate_entities(cdomain, 2, lambda x: np.logical_and(np.logical_and(x[1]<=small,x[1]-px+hef-dd>=-small),
x[0]*x[0]+x[2]*x[2]-r2<=small))
%%px
support_dofsc=fem.locate_dofs_topological(std.sub(1),2,support_facesc)
bc_supportc=fem.dirichletbc(0.,support_dofsc,std.sub(1))
stub1c_x_dofsc=fem.locate_dofs_topological(std.sub(0),2,stub_facesc)
bc_stub1c_x=fem.dirichletbc(0.,stub1c_x_dofsc,std.sub(0))
stub1c_z_dofs=fem.locate_dofs_topological(std.sub(2),2,stub_facesc)
bc_stub1c_z=fem.dirichletbc(0.,stub1c_z_dofs,std.sub(2))
if deri_load:
load_dofsc=fem.locate_dofs_topological(std.sub(1),2,load_facesc)
bc_loadc=fem.dirichletbc(imp,load_dofsc,std.sub(1))
bcsc=[bc_supportc,bc_stub1c_x,bc_stub1c_z,bc_loadc]
else:
bcsc=[bc_supportc,bc_stub1c_x,bc_stub1c_z]
Resolution#
%%px
problemc = petsc.LinearProblem(ac_ufl, bc_ufl,bcs=bcsc,petsc_options=petsc_options,petsc_options_prefix="c_")
dispc = problemc.solve()
dispc.name="disp_coarse"
#nopy
plt=pv.Plotter()
for i in range(nw):
plt.add_mesh(view['clippedc'][i],**pvoptpw)
plt.add_title("Coarse non enriched",font_size=15,color='white')
plt.camera_position = [(-33.0238996678173, -354.60971081475634, -386.6175118462205),
(49.764908472958275, -465.79163270141413, 466.73312098224534),
(0.002783776868024144, 0.9916497049175196, 0.12893065315643742)]
plt.add_axes()
plt.show()
Store coarse standard in coarse enriched field#
%%px
if nested_strategy:
disp_ce=dispc
else:
disp_ce.x.array[std_to_mix_np]=dispc.x.array
# print(disp_ce.x.array)
Twoscale loop#
Choose an enriched function#
%%px
enriched_shift=core.generateEnrichedShiftFunction(disp)
Resolution with linearBasicLoop (i.e. scale loop resolution)#
Dirichlet boundary condition must be defined on mixed space and merged into one Dirichlet BC object for non nested strategy
%%px
if nested_strategy:
BC_ce=bcsc
else:
support_dofsc=std_to_mix_np[support_dofsc]
stub1c_x_dofsc=std_to_mix_np[stub1c_x_dofsc]
stub1c_z_dofs=std_to_mix_np[stub1c_z_dofs]
if deri_load:
load_dofsc=std_to_mix_np[load_dofsc]
all_BC_values=fem.Function(coarse_enriched_space)
if deri_load:
all_BC_dof=np.concat((support_dofsc,stub1c_x_dofsc,stub1c_z_dofs,load_dofsc))
all_BC_values.x.array[load_dofsc]=imp
else:
all_BC_dof=np.concat((support_dofsc,stub1c_x_dofsc,stub1c_z_dofs))
all_BC_dof.sort()
BC_ce=[fem.dirichletbc(all_BC_values,all_BC_dof)]
Systems are created and TS resolution is done considering Dirichlet BC imposed at patch level:
%%px
[A,AD,BND,BD]=util.createFineScaleSytems(a_ufl,b_ufl,bcs=bcs,MPC=mpc)
%%px
dt=common.Timer("TS solve")
dt.start()
[dispf, ra,nm, it,hrb,hrr]=linear.linearBasicLoop(j,space,A,AD,BND,BD,enriched_shift,disp_ce,mpc,BC_ce,itmax,epsr)
dt.stop()
%%px
if nested_strategy:
dt=common.Timer("TS solve Eimp")
dt.start()
[dispf3, ra3,nm3, it3,hrb3,hrr3]=linear.linearBasicLoopEimp(j,space,A,AD,BND,BD,enriched_shift,disp_ce,mpc,BC_ce,10,itmax,epsr)
dt.stop()
Results#
Displacement#
Solutions with direct solver and TS solver are plotted below with their differences
convergence#
Text(0, 0.5, 'Relative residual')
Residual#
%%px
Rf=fem.Function(space)
R=Rf.x.petsc_vec
AD.mult(dispf.x.petsc_vec,R)
with io.VTXWriter(MPI.COMM_WORLD, "PO_Ax.bp", [Rf]) as ofile:
ofile.write(0.)
R.aypx(-1.,BD)
with io.VTXWriter(MPI.COMM_WORLD, "PO_Axb.bp", [Rf]) as ofile:
ofile.write(0.)
Video#
%%px
#nopy
[dispfv, r, it,pm,cm,hrb,hrr,hds]=linear.linearBasicLoopVideo("PO.gif",scale,j,disp.function_space,A,AD,BND,BD,enriched_shift,disp_ce,mpc,BC_ce,itmax,epsr,None,patch=False)
[stdout:0] At twoscale iteration 0 residual is 0.5910222721886301(b) 1.(r) 1.(ds)
At twoscale iteration 1 residual is 0.2816299962544677(b) 0.47651333884856195(r) 0.07553055014976795(ds)
At twoscale iteration 2 residual is 0.06251403193363851(b) 0.10577271767126668(r) 0.014900066431402605(ds)
At twoscale iteration 3 residual is 0.015027994449228584(b) 0.025427120358050167(r) 0.002057706656086269(ds)
At twoscale iteration 4 residual is 0.005304651070938736(b) 0.008975382689547288(r) 0.0005455883303826798(ds)
At twoscale iteration 5 residual is 0.0021662969076530645(b) 0.003665338870616489(r) 0.00019807656883422814(ds)
At twoscale iteration 6 residual is 0.0009732303564592294(b) 0.0016466898156904212(r) 8.034360186751716e-05(ds)
At twoscale iteration 7 residual is 0.00045980338475946875(b) 0.0007779797926341401(r) 3.4641116599365894e-05(ds)
At twoscale iteration 8 residual is 0.00022696498744742186(b) 0.0003840210396927037(r) 1.5544114380123173e-05(ds)
At twoscale iteration 9 residual is 0.00011577881294876418(b) 0.00019589585434745226(r) 7.19180626585551e-06(ds)
At twoscale iteration 10 residual is 6.0516717536985315e-05(b) 0.00010239329444029962(r) 3.4207021463252684e-06(ds)
At twoscale iteration 11 residual is 3.239058561635266e-05(b) 5.4804340107871444e-05(r) 1.6659285607131173e-06(ds)
At twoscale iteration 12 residual is 1.758275006658677e-05(b) 2.9749725000168992e-05(r) 8.315915616252038e-07(ds)
At twoscale iteration 13 residual is 9.705198674696362e-06(b) 1.6421037127343416e-05(r) 4.240381258163602e-07(ds)
At twoscale iteration 14 residual is 5.404177448500229e-06(b) 9.143779689533318e-06(r) 2.2109466603186254e-07(ds)
At twoscale iteration 15 residual is 3.0444184205130474e-06(b) 5.1511060813988975e-06(r) 1.1741807126188597e-07(ds)
At twoscale iteration 16 residual is 1.7256159073299897e-06(b) 2.9197138391076465e-06(r) 6.350920197494203e-08(ds)
At twoscale iteration 17 residual is 9.862214439272467e-07(b) 1.6686705228131997e-06(r) 3.48454940342552e-08(ds)
At twoscale iteration 18 residual is 5.663620765587203e-07(b) 9.582753530783366e-07(r) 1.9379492343828394e-08(ds)
At twoscale iteration 19 residual is 3.272422717941777e-07(b) 5.536885616549749e-07(r) 1.0888660023343596e-08(ds)
At twoscale iteration 20 residual is 1.8984584597984198e-07(b) 3.212160605670897e-07(r) 6.175875686633799e-09(ds)
At twoscale iteration 21 residual is 1.1065795710842277e-07(b) 1.8723145017639078e-07(r) 3.5275917919809405e-09(ds)
At twoscale iteration 22 residual is 6.472704987592293e-08(b) 1.0951710776690443e-07(r) 2.027947978786385e-09(ds)
Statistics#
[stdout:0] ****************************************************************************************************************************************************************
*** WIDEN YOUR WINDOW TO 160 CHARACTERS. Use 'enscript -r -fCourier9' to print this document ***
****************************************************************************************************************************************************************
------------------------------------------------------------------ PETSc Performance Summary: ------------------------------------------------------------------
Unknown Name on a named xxxxxxxxxxxxxxxxxxxxx with 8 processes, by salzman on xxxxxxxxxxxxxxxxxxx 2026
Using Petsc Release Version 3.22.2, unknown
Max Max/Min Avg Total
Time (sec): 2.345e+02 1.000 2.345e+02
Objects: 0.000e+00 0.000 0.000e+00
Flops: 8.379e+10 1.347 7.175e+10 5.740e+11
Flops/sec: 3.574e+08 1.347 3.060e+08 2.448e+09
MPI Msg Count: 1.702e+05 1.419 1.478e+05 1.182e+06
MPI Msg Len (bytes): 8.458e+08 2.040 3.793e+03 4.483e+09
MPI Reductions: 6.211e+05 1.014
Flop counting convention: 1 flop = 1 real number operation of type (multiply/divide/add/subtract)
e.g., VecAXPY() for real vectors of length N --> 2N flops
and VecAXPY() for complex vectors of length N --> 8N flops
Summary of Stages: ----- Time ------ ----- Flop ------ --- Messages --- -- Message Lengths -- -- Reductions --
Avg %Total Avg %Total Count %Total Avg %Total Count %Total
0: Main Stage: 3.2765e+01 14.0% 3.4765e+09 0.6% 1.128e+04 1.0% 8.899e+03 2.2% 1.601e+04 2.6%
1: generateCM: 2.8753e-01 0.1% 0.0000e+00 0.0% 4.480e+02 0.0% 1.212e+04 0.1% 3.800e+01 0.0%
2: generatePM: 6.9840e-02 0.0% 0.0000e+00 0.0% 0.000e+00 0.0% 0.000e+00 0.0% 0.000e+00 0.0%
3: genprbPM: 9.1757e-01 0.4% 0.0000e+00 0.0% 7.494e+03 0.6% 7.382e+01 0.0% 1.571e+04 2.5%
4: genstdCM: 4.5659e-01 0.2% 7.8036e+08 0.1% 1.340e+03 0.1% 4.062e+04 1.2% 9.300e+01 0.0%
5: solvePM: 4.6181e+00 2.0% 4.5445e+09 0.8% 1.200e+05 10.2% 2.132e+02 0.6% 1.315e+05 21.3%
6: resetCM: 1.8317e-04 0.0% 0.0000e+00 0.0% 0.000e+00 0.0% 0.000e+00 0.0% 0.000e+00 0.0%
7: PEupdateCM: 4.5934e-01 0.2% 0.0000e+00 0.0% 1.186e+04 1.0% 8.694e+02 0.2% 5.823e+03 0.9%
8: updateCM: 2.0762e+00 0.9% 9.3579e+09 1.6% 4.004e+03 0.3% 1.023e+05 9.1% 1.100e+02 0.0%
9: solveCM: 3.3462e+01 14.3% 2.0898e+11 36.4% 2.945e+05 24.9% 4.082e+03 26.8% 5.513e+04 8.9%
10: residual: 1.7179e-01 0.1% 8.6653e+08 0.2% 8.360e+02 0.1% 1.272e+04 0.2% 2.200e+01 0.0%
11: generateCMI: 2.8509e-01 0.1% 0.0000e+00 0.0% 4.480e+02 0.0% 1.212e+04 0.1% 3.800e+01 0.0%
12: generatePMI: 7.5345e-02 0.0% 0.0000e+00 0.0% 0.000e+00 0.0% 0.000e+00 0.0% 0.000e+00 0.0%
13: genprbPMI: 9.1640e-01 0.4% 0.0000e+00 0.0% 7.494e+03 0.6% 7.382e+01 0.0% 1.571e+04 2.5%
14: genstdCMI: 4.3080e-01 0.2% 7.8036e+08 0.1% 1.340e+03 0.1% 4.062e+04 1.2% 9.300e+01 0.0%
15: solvePMI: 1.6459e+01 7.0% 5.1555e+09 0.9% 1.354e+05 11.5% 2.135e+02 0.6% 1.493e+05 24.1%
16: resetCMI: 2.0428e-04 0.0% 0.0000e+00 0.0% 0.000e+00 0.0% 0.000e+00 0.0% 0.000e+00 0.0%
17: PEupdateCMI: 5.5849e-01 0.2% 0.0000e+00 0.0% 1.334e+04 1.1% 8.758e+02 0.3% 5.826e+03 0.9%
18: updateCMI: 1.0752e+00 0.5% 4.6790e+09 0.8% 2.002e+03 0.2% 1.023e+05 4.6% 5.500e+01 0.0%
19: solveCMI: 3.5713e+01 15.2% 1.0257e+11 17.9% 1.399e+05 11.8% 4.469e+03 13.9% 2.746e+04 4.4%
20: residualI: 2.7531e-01 0.1% 9.8470e+08 0.2% 9.500e+02 0.1% 1.272e+04 0.3% 2.500e+01 0.0%
21: updateEImpCMI: 6.2156e-02 0.0% 3.1121e+08 0.1% 8.120e+02 0.1% 8.666e+03 0.2% 0.000e+00 0.0%
22: solveEImpCMI: 1.3294e+00 0.6% 1.7000e+09 0.3% 1.275e+03 0.1% 4.372e+03 0.1% 2.100e+01 0.0%
23: generateCMV: 3.1594e-01 0.1% 0.0000e+00 0.0% 4.480e+02 0.0% 1.212e+04 0.1% 3.800e+01 0.0%
24: generatePMV: 7.3499e-02 0.0% 0.0000e+00 0.0% 0.000e+00 0.0% 0.000e+00 0.0% 0.000e+00 0.0%
25: solvePMV: 2.4190e+01 10.3% 4.5445e+09 0.8% 1.200e+05 10.2% 2.133e+02 0.6% 1.315e+05 21.3%
26: resetCMV: 1.8846e-04 0.0% 0.0000e+00 0.0% 0.000e+00 0.0% 0.000e+00 0.0% 0.000e+00 0.0%
27: PEupdateCMV: 5.1888e-01 0.2% 0.0000e+00 0.0% 1.186e+04 1.0% 8.694e+02 0.2% 5.823e+03 0.9%
28: updateCMV: 2.3290e+00 1.0% 9.3579e+09 1.6% 4.004e+03 0.3% 1.023e+05 9.1% 1.100e+02 0.0%
29: solveCMV: 7.3927e+01 31.5% 2.1507e+11 37.5% 2.900e+05 24.5% 4.297e+03 27.8% 5.529e+04 8.9%
30: residualV: 1.8178e-01 0.1% 8.6653e+08 0.2% 8.360e+02 0.1% 1.272e+04 0.2% 2.200e+01 0.0%
------------------------------------------------------------------------------------------------------------------------
See the 'Profiling' chapter of the users' manual for details on interpreting output.
Phase summary info:
Count: number of times phase was executed
Time and Flop: Max - maximum over all processors
Ratio - ratio of maximum to minimum over all processors
Mess: number of messages sent
AvgLen: average message length (bytes)
Reduct: number of global reductions
Global: entire computation
Stage: stages of a computation. Set stages with PetscLogStagePush() and PetscLogStagePop().
%T - percent time in this phase %F - percent flop in this phase
%M - percent messages in this phase %L - percent message lengths in this phase
%R - percent reductions in this phase
Total Mflop/s: 10e-6 * (sum of flop over all processors)/(max time over all processors)
------------------------------------------------------------------------------------------------------------------------
Event Count Time (sec) Flop --- Global --- --- Stage ---- Total
Max Ratio Max Ratio Max Ratio Mess AvgLen Reduct %T %F %M %L %R %T %F %M %L %R Mflop/s
------------------------------------------------------------------------------------------------------------------------
--- Event Stage 0: Main Stage
BuildTwoSided 679 1.3 1.3571e+00 15.4 0.00e+00 0.0 3.2e+03 4.0e+00 5.8e+02 0 0 0 0 0 3 0 29 0 4 0
BuildTwoSidedF 338 1.3 1.3524e+00 15.9 0.00e+00 0.0 5.4e+02 7.6e+04 2.9e+02 0 0 0 1 0 3 0 5 41 2 0
SFSetGraph 341 1.3 3.5148e-03 2.7 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFSetUp 341 1.3 2.2707e-02 2.4 0.00e+00 0.0 5.9e+03 6.4e+02 2.9e+02 0 0 0 0 0 0 0 52 4 2 0
SFReduceBegin 2 1.0 1.6854e-05 1.2 0.00e+00 0.0 4.2e+01 4.9e+03 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFReduceEnd 2 1.0 1.5167e-05 1.7 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFPack 16 1.0 3.1649e-04 1.3 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFUnpack 16 1.0 2.2984e-04 2.9 8.48e+03 1.4 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 250
VecDot 23 1.0 6.6318e-03 4.6 3.77e+06 2.6 0.0e+00 0.0e+00 2.3e+01 0 0 0 0 0 0 1 0 0 0 2890
VecNorm 4 1.0 9.1270e-03 22.2 6.56e+05 2.6 0.0e+00 0.0e+00 4.0e+00 0 0 0 0 0 0 0 0 0 0 365
VecSet 7 1.0 1.9741e-04 2.4 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
VecAYPX 4 1.0 2.5946e-04 1.5 3.28e+05 2.6 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 6424
VecScatterBegin 14 1.0 1.0520e-02 16.2 0.00e+00 0.0 5.3e+02 1.1e+04 0.0e+00 0 0 0 0 0 0 0 5 6 0 0
VecScatterEnd 14 1.0 2.7855e-02 3.5 8.48e+03 1.4 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 2
MatMult 4 1.0 3.6199e-02 1.2 2.99e+07 2.5 1.5e+02 1.3e+04 0.0e+00 0 0 0 0 0 0 4 1 2 0 4214
MatMultAdd 3 1.0 4.7966e-03 2.2 3.44e+06 2.6 1.1e+02 7.4e+02 0.0e+00 0 0 0 0 0 0 0 1 0 0 3578
MatMultTranspose 1 1.0 1.4609e-03 1.0 1.15e+06 2.6 3.8e+01 7.4e+02 0.0e+00 0 0 0 0 0 0 0 0 0 0 3918
MatSolve 2 1.0 1.5231e+00 1.0 1.55e+08 1.3 5.0e+02 1.8e+04 1.0e+01 1 0 0 0 0 5 31 4 9 0 717
MatLUFactorSym 2 1.0 1.3210e+00 1.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 1 0 0 0 0 4 0 0 0 0 0
MatLUFactorNum 2 1.0 7.7517e+00 1.0 1.43e+08 1.2 0.0e+00 0.0e+00 0.0e+00 3 0 0 0 0 24 30 0 0 0 134
MatAssemblyBegin 2599 1.1 1.4275e+00 8.3 0.00e+00 0.0 5.4e+02 7.6e+04 2.9e+02 0 0 0 1 0 3 0 5 41 2 0
MatAssemblyEnd 2599 1.1 9.7575e-02 1.1 1.90e+05 1.4 4.8e+03 3.2e+02 1.4e+03 0 0 0 0 0 0 0 42 2 9 13
MatCreateSubMats 1934 1.0 8.6224e-01 1.1 0.00e+00 0.0 4.1e+03 1.2e+02 1.1e+04 0 0 0 0 2 2 0 36 1 67 0
MatZeroEntries 9 1.0 4.9212e-02 4.3 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
MatTranspose 6 1.0 3.3496e-02 1.6 0.00e+00 0.0 2.9e+02 1.4e+04 1.8e+01 0 0 0 0 0 0 0 3 4 0 0
MatMatMultSym 2 1.0 9.7290e-02 1.0 0.00e+00 0.0 2.9e+02 4.7e+04 2.0e+01 0 0 0 0 0 0 0 3 14 0 0
MatMatMultNum 2 1.0 5.8588e-02 1.0 6.80e+07 2.5 5.8e+01 1.5e+05 2.0e+00 0 0 0 0 0 0 10 1 9 0 5804
MatRARtSym 1 1.0 1.2670e-01 1.0 0.00e+00 0.0 4.0e+02 4.1e+04 2.9e+01 0 0 0 0 0 0 0 4 16 0 0
MatRARtNum 1 1.0 6.3536e-02 1.0 6.64e+07 2.5 5.8e+01 1.5e+05 2.0e+00 0 0 0 0 0 0 9 1 9 0 5185
MatTrnMatMultSym 1 1.0 7.4665e-02 1.0 0.00e+00 0.0 2.7e+02 1.2e+04 1.3e+01 0 0 0 0 0 0 0 2 3 0 0
MatTrnMatMultNum 1 1.0 3.5466e-02 1.0 2.05e+07 2.6 7.6e+01 3.3e+04 1.0e+00 0 0 0 0 0 0 3 1 2 0 2965
MatGetLocalMat 10 1.0 5.5383e-02 1.9 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
MatGetBrAoCol 8 1.0 8.9766e-02 6.3 0.00e+00 0.0 4.6e+02 9.6e+04 0.0e+00 0 0 0 1 0 0 0 4 44 0 0
PCSetUp 2 1.0 9.0730e+00 1.0 1.43e+08 1.2 0.0e+00 0.0e+00 4.0e+00 4 0 0 0 0 28 30 0 0 0 114
PCApply 2 1.0 1.5231e+00 1.0 1.55e+08 1.3 5.0e+02 1.8e+04 1.0e+01 1 0 0 0 0 5 31 4 9 0 717
KSPSetUp 2 1.0 2.1540e-06 1.8 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
KSPSolve 2 1.0 1.5232e+00 1.0 1.55e+08 1.3 5.0e+02 1.8e+04 1.2e+01 1 0 0 0 0 5 31 4 9 0 717
--- Event Stage 1: generateCM
BuildTwoSided 8 1.0 1.0456e-01 79.1 0.00e+00 0.0 2.1e+02 4.0e+00 8.0e+00 0 0 0 0 0 22 0 48 0 21 0
BuildTwoSidedF 3 1.0 1.0421e-01 9.8 0.00e+00 0.0 1.6e+02 3.3e+04 3.0e+00 0 0 0 0 0 23 0 35 96 8 0
SFSetGraph 5 1.0 6.0782e-05 1.7 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFSetUp 5 1.0 1.1121e-02 27.9 0.00e+00 0.0 2.7e+02 7.9e+02 5.0e+00 0 0 0 0 0 3 0 61 4 13 0
SFReduceBegin 1 1.0 1.0408e-05 1.7 0.00e+00 0.0 2.0e+01 3.0e+02 0.0e+00 0 0 0 0 0 0 0 4 0 0 0
SFReduceEnd 1 1.0 3.6990e-06 2.2 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFPack 1 1.0 1.6300e-07 1.6 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFUnpack 1 1.0 6.0600e-07 3.4 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
VecNorm 1 1.0 1.4735e-04 7.6 0.00e+00 0.0 0.0e+00 0.0e+00 1.0e+00 0 0 0 0 0 0 0 0 0 3 0
MatAssemblyBegin 3 1.0 1.0779e-01 7.2 0.00e+00 0.0 1.6e+02 3.3e+04 3.0e+00 0 0 0 0 0 24 0 35 96 8 0
MatAssemblyEnd 3 1.0 3.0040e-02 1.5 0.00e+00 0.0 1.6e+02 1.3e+03 1.5e+01 0 0 0 0 0 9 0 35 4 39 0
MatTranspose 1 1.0 1.9881e-02 1.0 0.00e+00 0.0 1.4e+02 2.1e+04 9.0e+00 0 0 0 0 0 7 0 31 54 24 0
--- Event Stage 2: generatePM
--- Event Stage 3: genprbPM
BuildTwoSided 628 1.3 8.0761e-03 1.5 0.00e+00 0.0 2.0e+03 4.0e+00 5.3e+02 0 0 0 0 0 1 0 27 1 3 0
BuildTwoSidedF 314 1.3 7.3701e-03 1.5 0.00e+00 0.0 0.0e+00 0.0e+00 2.6e+02 0 0 0 0 0 1 0 0 0 2 0
SFSetGraph 314 1.3 2.1723e-04 1.4 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFSetUp 314 1.3 1.9876e-03 1.3 0.00e+00 0.0 4.1e+03 1.2e+02 2.6e+02 0 0 0 0 0 0 0 55 92 2 0
MatAssemblyBegin 2562 1.1 8.2574e-03 1.4 0.00e+00 0.0 0.0e+00 0.0e+00 2.6e+02 0 0 0 0 0 1 0 0 0 2 0
MatAssemblyEnd 2562 1.1 1.9463e-02 1.4 0.00e+00 0.0 4.1e+03 1.2e+02 1.3e+03 0 0 0 0 0 2 0 55 92 8 0
MatCreateSubMats 1934 1.0 6.8838e-01 1.2 0.00e+00 0.0 4.1e+03 1.2e+02 1.1e+04 0 0 0 0 2 68 0 55 92 68 0
--- Event Stage 4: genstdCM
BuildTwoSided 22 1.0 1.0176e-01 2.6 0.00e+00 0.0 3.4e+02 4.0e+00 2.2e+01 0 0 0 0 0 15 0 25 0 24 0
BuildTwoSidedF 13 1.0 1.1128e-01 2.4 0.00e+00 0.0 2.4e+02 3.6e+04 1.3e+01 0 0 0 0 0 16 0 18 16 14 0
SFSetGraph 9 1.0 6.0517e-05 1.5 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFSetUp 9 1.0 1.4757e-03 3.6 0.00e+00 0.0 4.5e+02 9.5e+02 9.0e+00 0 0 0 0 0 0 0 33 1 10 0
SFReduceBegin 2 1.0 4.5243e-04 20.0 0.00e+00 0.0 4.2e+01 4.9e+03 0.0e+00 0 0 0 0 0 0 0 3 0 0 0
SFReduceEnd 2 1.0 6.8833e-04 54.3 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFPack 3 1.0 5.9400e-07 1.9 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFUnpack 3 1.0 1.4585e-05 1.5 5.43e+02 1.6 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 240
VecScatterBegin 1 1.0 5.3973e-05 2.4 0.00e+00 0.0 3.8e+01 7.4e+02 0.0e+00 0 0 0 0 0 0 0 3 0 0 0
VecScatterEnd 1 1.0 5.4804e-04 49.2 5.43e+02 1.6 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 6
MatMultTranspose 1 1.0 1.4703e-03 1.0 1.15e+06 2.6 3.8e+01 7.4e+02 0.0e+00 0 0 0 0 0 0 1 3 0 0 3893
MatAssemblyBegin 26 1.0 1.1411e-01 2.3 0.00e+00 0.0 2.4e+02 3.6e+04 1.3e+01 0 0 0 0 0 17 0 18 16 14 0
MatAssemblyEnd 26 1.0 3.3938e-02 1.5 2.58e+04 1.7 3.6e+02 6.1e+02 3.7e+01 0 0 0 0 0 6 0 27 0 40 5
MatZeroEntries 2 1.0 4.0810e-03 2.6 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 1 0 0 0 0 0
MatTranspose 6 1.0 2.7172e-02 1.4 0.00e+00 0.0 2.9e+02 1.4e+04 1.8e+01 0 0 0 0 0 5 0 22 7 19 0
MatMatMultSym 2 1.0 9.8302e-02 1.0 0.00e+00 0.0 2.9e+02 4.7e+04 2.0e+01 0 0 0 0 0 22 0 22 26 22 0
MatMatMultNum 2 1.0 7.4073e-02 1.0 6.80e+07 2.5 5.8e+01 1.5e+05 2.0e+00 0 0 0 0 0 16 44 4 16 2 4591
MatRARtSym 1 1.0 9.7966e-02 1.0 0.00e+00 0.0 4.0e+02 4.1e+04 2.9e+01 0 0 0 0 0 21 0 30 30 31 0
MatRARtNum 1 1.0 6.3192e-02 1.0 6.64e+07 2.5 5.8e+01 1.5e+05 2.0e+00 0 0 0 0 0 14 42 4 16 2 5213
MatTrnMatMultSym 1 1.0 8.3580e-02 1.0 0.00e+00 0.0 2.7e+02 1.2e+04 1.3e+01 0 0 0 0 0 18 0 20 6 14 0
MatTrnMatMultNum 1 1.0 3.5227e-02 1.0 2.05e+07 2.6 7.6e+01 3.3e+04 1.0e+00 0 0 0 0 0 8 13 6 5 1 2985
MatGetLocalMat 10 1.0 5.5547e-02 2.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 9 0 0 0 0 0
MatGetBrAoCol 8 1.0 1.0320e-01 3.7 0.00e+00 0.0 4.6e+02 9.6e+04 0.0e+00 0 0 0 1 0 13 0 35 82 0 0
--- Event Stage 5: solvePM
BuildTwoSided 21431 1.0 1.8730e-01 1.6 0.00e+00 0.0 2.4e+04 4.0e+00 2.1e+04 0 0 2 0 3 4 0 20 0 16 0
SFSetGraph 21431 1.0 1.7056e-02 2.1 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFSetUp 21431 1.0 2.2302e-01 1.3 0.00e+00 0.0 4.8e+04 1.9e+01 2.1e+04 0 0 4 0 3 4 0 40 3 16 0
SFPack 28182 1.1 4.0220e-03 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFUnpack 28182 1.1 2.0267e-03 1.1 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
VecCopy 18348 2.9 6.3931e-03 2.5 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
VecAYPX 21274 2.3 9.6721e-03 1.9 4.43e+06 2.5 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 2265
VecScatterBegin 28182 1.1 1.0568e-01 1.4 0.00e+00 0.0 6.7e+04 3.6e+02 0.0e+00 0 0 6 1 0 2 0 56 93 0 0
VecScatterEnd 28182 1.1 4.0724e-02 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 1 0 0 0 0 0
MatMult 21274 2.3 1.3168e-01 1.7 1.21e+08 2.4 2.3e+04 4.1e+02 0.0e+00 0 0 2 0 0 2 13 19 37 0 4578
MatSolve 21274 2.3 1.5132e+00 1.8 7.76e+08 2.5 2.8e+04 4.3e+02 6.6e+02 0 1 2 0 0 24 85 24 48 0 2548
MatLUFactorSym 967 2.3 2.6574e-01 1.8 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 4 0 0 0 0 0
MatLUFactorNum 967 2.3 4.2842e-01 2.1 1.28e+07 2.5 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 6 1 0 0 0 148
PCSetUp 967 2.3 7.0531e-01 2.0 1.28e+07 2.5 0.0e+00 0.0e+00 2.6e+02 0 0 0 0 0 10 1 0 0 0 90
PCApply 21274 2.3 1.5247e+00 1.8 7.76e+08 2.5 2.8e+04 4.3e+02 6.6e+02 0 1 2 0 0 24 85 24 48 0 2529
KSPSetUp 967 2.3 3.6157e-04 2.4 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
KSPSolve 21274 2.3 1.5552e+00 1.8 7.76e+08 2.5 2.8e+04 4.3e+02 3.5e+03 0 1 2 0 1 24 85 24 48 3 2479
--- Event Stage 6: resetCM
--- Event Stage 7: PEupdateCM
BuildTwoSided 1955 1.0 1.3886e-02 1.3 0.00e+00 0.0 4.9e+02 4.0e+00 2.0e+03 0 0 0 0 0 3 0 4 0 34 0
BuildTwoSidedF 988 1.0 4.0223e-03 1.0 0.00e+00 0.0 0.0e+00 0.0e+00 9.9e+02 0 0 0 0 0 1 0 0 0 17 0
SFSetGraph 967 1.0 1.1905e-03 2.3 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFSetUp 967 1.0 1.2045e-02 1.2 0.00e+00 0.0 9.9e+02 2.3e+02 9.7e+02 0 0 0 0 0 2 0 8 2 17 0
SFPack 21274 1.0 1.4891e-03 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFUnpack 21274 1.0 9.7865e-04 1.1 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
VecScatterBegin 21274 1.0 2.3143e-02 1.4 0.00e+00 0.0 1.1e+04 9.3e+02 0.0e+00 0 0 1 0 0 4 0 92 98 0 0
VecScatterEnd 21274 1.0 7.6172e-03 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 2 0 0 0 0 0
MatAssemblyBegin 988 1.0 5.6782e-03 1.0 0.00e+00 0.0 0.0e+00 0.0e+00 9.9e+02 0 0 0 0 0 1 0 0 0 17 0
MatAssemblyEnd 988 1.0 6.6553e-04 1.1 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
--- Event Stage 8: updateCM
BuildTwoSided 110 1.0 1.1280e+00 16.0 0.00e+00 0.0 9.2e+02 4.0e+00 1.1e+02 0 0 0 0 0 35 0 23 0 100 0
BuildTwoSidedF 110 1.0 1.1285e+00 15.9 0.00e+00 0.0 1.8e+03 4.3e+04 1.1e+02 0 0 0 2 0 35 0 46 19 100 0
SFPack 22 1.0 1.1559e-03 1.5 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFUnpack 22 1.0 4.5730e-06 1.6 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
VecScatterBegin 22 1.0 1.0914e-02 9.8 0.00e+00 0.0 4.4e+02 1.8e+04 0.0e+00 0 0 0 0 0 0 0 11 2 0 0
VecScatterEnd 22 1.0 3.1935e-02 61.8 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 1 0 0 0 0 0
MatMult 22 1.0 5.4390e-02 1.6 2.43e+07 2.5 4.4e+02 1.8e+04 0.0e+00 0 0 0 0 0 2 1 11 2 0 2211
MatAssemblyBegin 110 1.0 1.1562e+00 11.9 0.00e+00 0.0 1.8e+03 4.3e+04 1.1e+02 0 0 0 2 0 36 0 46 19 100 0
MatAssemblyEnd 110 1.0 6.0092e-02 1.9 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 2 0 0 0 0 0
MatTranspose 44 1.0 1.7828e-01 1.1 0.00e+00 0.0 1.8e+03 4.3e+04 4.4e+01 0 0 0 2 0 8 0 46 19 40 0
MatMatMultNum 22 1.0 5.1473e-01 1.0 4.02e+08 2.5 4.4e+02 3.1e+05 2.2e+01 0 0 0 3 0 25 21 11 33 20 3867
MatRARtNum 22 1.0 1.4827e+00 1.0 1.46e+09 2.5 2.2e+03 1.2e+05 6.6e+01 1 1 0 6 0 71 77 54 61 60 4888
MatGetLocalMat 66 1.0 2.1289e-01 2.7 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 6 0 0 0 0 0
MatGetBrAoCol 66 1.0 6.3642e-01 10.3 0.00e+00 0.0 1.7e+03 1.9e+05 0.0e+00 0 0 0 7 0 19 0 43 79 0 0
--- Event Stage 9: solveCM
BuildTwoSided 25 1.0 3.6922e-03 2.6 0.00e+00 0.0 6.6e+02 4.0e+00 2.5e+01 0 0 0 0 0 0 0 0 0 0 0
SFSetGraph 25 1.0 2.8774e-03 2.8 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFSetUp 25 1.0 4.2901e-03 1.7 0.00e+00 0.0 1.3e+03 9.2e+02 2.5e+01 0 0 0 0 0 0 0 0 0 0 0
SFPack 16636 1.0 6.6476e-02 2.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFUnpack 16636 1.0 3.5471e-02 4.4 1.52e+05 1.5 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 28
VecDot 24163 1.0 8.1334e-01 1.6 1.17e+09 7.1 0.0e+00 0.0e+00 4.8e+04 0 1 0 0 8 2 2 0 0 88 3953
VecMDot 1602 1.0 5.8433e-01 1.2 1.09e+09 7.1 0.0e+00 0.0e+00 4.5e+04 0 1 0 0 7 2 1 0 0 82 5122
VecNorm 1668 1.0 2.3309e-01 14.7 8.08e+07 7.1 0.0e+00 0.0e+00 3.3e+03 0 0 0 0 1 1 0 0 0 6 952
VecScale 1668 1.0 7.4274e-03 1.8 4.04e+07 7.1 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 14942
VecCopy 66 1.0 1.2158e-03 3.4 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
VecSet 5092 1.0 2.4682e-02 2.9 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
VecAXPY 24207 1.0 2.5010e-01 8.7 1.17e+09 7.1 0.0e+00 0.0e+00 0.0e+00 0 1 0 0 0 0 2 0 0 0 12880
VecMAXPY 1668 1.0 2.5375e-01 7.9 1.17e+09 7.1 0.0e+00 0.0e+00 0.0e+00 0 1 0 0 0 0 2 0 0 0 12636
VecScatterBegin 16636 1.0 4.5890e-01 2.0 0.00e+00 0.0 2.9e+05 4.1e+03 0.0e+00 0 0 25 27 0 1 0 99 100 0 0
VecScatterEnd 16636 1.0 1.6224e+00 5.2 1.52e+05 1.5 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 3 0 0 0 0 1
VecNormalize 1668 1.0 2.3966e-01 9.5 1.21e+08 7.1 0.0e+00 0.0e+00 3.3e+03 0 0 0 0 1 1 0 0 0 6 1389
MatMult 1646 1.0 2.9477e+00 1.1 3.69e+09 4.6 1.7e+05 2.5e+03 3.3e+03 1 2 14 9 1 8 6 57 35 6 4422
MatMultAdd 6606 1.0 2.7769e+00 1.2 3.70e+09 4.5 1.7e+05 2.5e+03 0.0e+00 1 2 14 9 0 8 6 57 35 0 4740
MatMultTrAdd 22 1.0 1.3069e-01 8.7 2.44e+07 2.5 4.4e+02 1.8e+04 0.0e+00 0 0 0 0 0 0 0 0 1 0 930
MatSolve 3336 1.0 2.8659e+01 1.0 2.97e+10 1.4 1.3e+05 6.2e+03 1.2e+02 12 33 11 17 0 85 90 43 65 0 6592
MatLUFactorSym 2 1.0 1.5557e-01 1.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
MatLUFactorNum 23 1.0 8.7225e-01 1.0 3.18e+07 1.2 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 3 0 0 0 0 259
MatCreateSubMat 2 1.0 1.1005e-04 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 4.0e+00 0 0 0 0 0 0 0 0 0 0 0
PCSetUp 22 1.0 3.3416e-01 1.0 1.96e+07 1.5 4.8e+01 1.1e+04 2.3e+01 0 0 0 0 0 1 0 0 0 0 373
PCSetUpOnBlocks 22 1.0 6.9556e-01 1.0 1.39e+07 1.2 0.0e+00 0.0e+00 4.0e+00 0 0 0 0 0 2 0 0 0 0 145
PCApply 1668 1.0 2.9229e+01 1.0 2.97e+10 1.4 1.3e+05 6.2e+03 3.5e+03 12 33 11 17 1 87 90 43 65 6 6468
KSPSolve_FS_0 1668 1.0 2.5460e+01 1.0 2.79e+10 1.4 7.5e+04 8.2e+03 1.7e+03 11 30 6 14 0 76 83 26 51 3 6844
KSPSolve_FS_1 1668 1.0 3.6003e+00 1.0 2.23e+09 1.3 5.0e+04 3.2e+03 1.8e+03 2 3 4 4 0 11 7 17 13 3 4107
KSPSetUp 1 1.0 9.0738e-03 73.3 0.00e+00 0.0 0.0e+00 0.0e+00 5.0e+00 0 0 0 0 0 0 0 0 0 0 0
KSPSolve 22 1.0 3.2667e+01 1.0 3.17e+10 1.3 2.9e+05 4.1e+03 5.5e+04 14 36 25 27 9 98 100 100 99 100 6386
KSPGMRESOrthog 1602 1.0 7.4955e-01 1.2 2.18e+09 7.1 0.0e+00 0.0e+00 4.5e+04 0 1 0 0 7 2 3 0 0 82 7987
--- Event Stage 10: residual
SFPack 22 1.0 4.6019e-04 2.3 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFUnpack 22 1.0 1.2925e-05 4.2 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
VecNorm 22 1.0 1.1884e-01 13.1 3.61e+06 2.6 0.0e+00 0.0e+00 2.2e+01 0 0 0 0 0 30 2 0 0 100 154
VecAYPX 22 1.0 3.0130e-03 4.0 1.80e+06 2.6 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 1 1 0 0 0 3043
VecScatterBegin 22 1.0 1.2048e-03 3.5 0.00e+00 0.0 8.4e+02 1.3e+04 0.0e+00 0 0 0 0 0 0 0 100 100 0 0
VecScatterEnd 22 1.0 7.8711e-02 232.6 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 14 0 0 0 0 0
MatMult 22 1.0 1.6806e-01 2.1 1.64e+08 2.5 8.4e+02 1.3e+04 0.0e+00 0 0 0 0 0 69 97 100 100 0 4992
--- Event Stage 11: generateCMI
BuildTwoSided 8 1.0 1.1009e-01 70.7 0.00e+00 0.0 2.1e+02 4.0e+00 8.0e+00 0 0 0 0 0 23 0 48 0 21 0
BuildTwoSidedF 3 1.0 1.1001e-01 65.3 0.00e+00 0.0 1.6e+02 3.3e+04 3.0e+00 0 0 0 0 0 23 0 35 96 8 0
SFSetGraph 5 1.0 6.0274e-05 1.7 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFSetUp 5 1.0 4.6477e-04 2.0 0.00e+00 0.0 2.7e+02 7.9e+02 5.0e+00 0 0 0 0 0 0 0 61 4 13 0
SFReduceBegin 1 1.0 9.0370e-06 1.4 0.00e+00 0.0 2.0e+01 3.0e+02 0.0e+00 0 0 0 0 0 0 0 4 0 0 0
SFReduceEnd 1 1.0 2.8460e-06 1.4 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFPack 1 1.0 1.5600e-07 1.6 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFUnpack 1 1.0 6.5300e-07 3.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
VecNorm 1 1.0 2.3518e-04 13.2 0.00e+00 0.0 0.0e+00 0.0e+00 1.0e+00 0 0 0 0 0 0 0 0 0 3 0
MatAssemblyBegin 3 1.0 1.1399e-01 17.4 0.00e+00 0.0 1.6e+02 3.3e+04 3.0e+00 0 0 0 0 0 25 0 35 96 8 0
MatAssemblyEnd 3 1.0 9.6655e-03 1.1 0.00e+00 0.0 1.6e+02 1.3e+03 1.5e+01 0 0 0 0 0 3 0 35 4 39 0
MatTranspose 1 1.0 9.6320e-03 1.0 0.00e+00 0.0 1.4e+02 2.1e+04 9.0e+00 0 0 0 0 0 3 0 31 54 24 0
--- Event Stage 12: generatePMI
--- Event Stage 13: genprbPMI
BuildTwoSided 628 1.3 8.8015e-03 1.5 0.00e+00 0.0 2.0e+03 4.0e+00 5.3e+02 0 0 0 0 0 1 0 27 1 3 0
BuildTwoSidedF 314 1.3 8.0004e-03 1.5 0.00e+00 0.0 0.0e+00 0.0e+00 2.6e+02 0 0 0 0 0 1 0 0 0 2 0
SFSetGraph 314 1.3 2.3065e-04 1.3 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFSetUp 314 1.3 2.1381e-03 1.4 0.00e+00 0.0 4.1e+03 1.2e+02 2.6e+02 0 0 0 0 0 0 0 55 92 2 0
MatAssemblyBegin 2562 1.1 9.0187e-03 1.5 0.00e+00 0.0 0.0e+00 0.0e+00 2.6e+02 0 0 0 0 0 1 0 0 0 2 0
MatAssemblyEnd 2562 1.1 2.0708e-02 1.4 0.00e+00 0.0 4.1e+03 1.2e+02 1.3e+03 0 0 0 0 0 2 0 55 92 8 0
MatCreateSubMats 1934 1.0 6.8563e-01 1.2 0.00e+00 0.0 4.1e+03 1.2e+02 1.1e+04 0 0 0 0 2 68 0 55 92 68 0
--- Event Stage 14: genstdCMI
BuildTwoSided 22 1.0 9.0537e-02 4.8 0.00e+00 0.0 3.4e+02 4.0e+00 2.2e+01 0 0 0 0 0 14 0 25 0 24 0
BuildTwoSidedF 13 1.0 9.0227e-02 4.8 0.00e+00 0.0 2.4e+02 3.6e+04 1.3e+01 0 0 0 0 0 14 0 18 16 14 0
SFSetGraph 9 1.0 6.4168e-05 1.5 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFSetUp 9 1.0 7.4993e-04 2.0 0.00e+00 0.0 4.5e+02 9.5e+02 9.0e+00 0 0 0 0 0 0 0 33 1 10 0
SFReduceBegin 2 1.0 1.8681e-05 1.3 0.00e+00 0.0 4.2e+01 4.9e+03 0.0e+00 0 0 0 0 0 0 0 3 0 0 0
SFReduceEnd 2 1.0 1.6644e-05 1.6 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFPack 3 1.0 7.4100e-07 2.1 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFUnpack 3 1.0 1.0247e-05 1.7 5.43e+02 1.6 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 342
VecScatterBegin 1 1.0 4.9916e-05 2.2 0.00e+00 0.0 3.8e+01 7.4e+02 0.0e+00 0 0 0 0 0 0 0 3 0 0 0
VecScatterEnd 1 1.0 6.0967e-04 70.8 5.43e+02 1.6 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 6
MatMultTranspose 1 1.0 1.5998e-03 1.0 1.15e+06 2.6 3.8e+01 7.4e+02 0.0e+00 0 0 0 0 0 0 1 3 0 0 3578
MatConvert 2 1.0 4.3251e-03 9.8 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
MatAssemblyBegin 26 1.0 9.3141e-02 4.3 0.00e+00 0.0 2.4e+02 3.6e+04 1.3e+01 0 0 0 0 0 15 0 18 16 14 0
MatAssemblyEnd 26 1.0 2.6481e-02 1.1 2.58e+04 1.7 3.6e+02 6.1e+02 3.7e+01 0 0 0 0 0 6 0 27 0 40 6
MatZeroEntries 2 1.0 3.7047e-03 4.4 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 1 0 0 0 0 0
MatTranspose 6 1.0 3.0006e-02 1.6 0.00e+00 0.0 2.9e+02 1.4e+04 1.8e+01 0 0 0 0 0 5 0 22 7 19 0
MatMatMultSym 2 1.0 8.9209e-02 1.0 0.00e+00 0.0 2.9e+02 4.7e+04 2.0e+01 0 0 0 0 0 21 0 22 26 22 0
MatMatMultNum 2 1.0 6.3902e-02 1.0 6.80e+07 2.5 5.8e+01 1.5e+05 2.0e+00 0 0 0 0 0 15 44 4 16 2 5322
MatRARtSym 1 1.0 9.6765e-02 1.0 0.00e+00 0.0 4.0e+02 4.1e+04 2.9e+01 0 0 0 0 0 22 0 30 30 31 0
MatRARtNum 1 1.0 5.9881e-02 1.0 6.64e+07 2.5 5.8e+01 1.5e+05 2.0e+00 0 0 0 0 0 14 42 4 16 2 5501
MatTrnMatMultSym 1 1.0 7.6622e-02 1.0 0.00e+00 0.0 2.7e+02 1.2e+04 1.3e+01 0 0 0 0 0 18 0 20 6 14 0
MatTrnMatMultNum 1 1.0 3.5669e-02 1.0 2.05e+07 2.6 7.6e+01 3.3e+04 1.0e+00 0 0 0 0 0 8 13 6 5 1 2948
MatGetLocalMat 10 1.0 5.5626e-02 2.1 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 10 0 0 0 0 0
MatGetBrAoCol 8 1.0 6.4916e-02 2.1 0.00e+00 0.0 4.6e+02 9.6e+04 0.0e+00 0 0 0 1 0 10 0 35 82 0 0
--- Event Stage 15: solvePMI
BuildTwoSided 24332 1.0 2.5878e-01 1.5 0.00e+00 0.0 2.7e+04 4.0e+00 2.4e+04 0 0 2 0 4 1 0 20 0 16 0
SFSetGraph 24332 1.0 1.9474e-02 2.1 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFSetUp 24332 1.0 3.0006e-01 1.3 0.00e+00 0.0 5.4e+04 1.8e+01 2.4e+04 0 0 5 0 4 2 0 40 3 16 0
SFPack 32025 1.1 6.6922e-03 1.5 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFUnpack 32025 1.1 2.3891e-03 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
VecCopy 20850 2.9 7.1138e-03 2.5 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
VecAYPX 24175 2.3 1.3289e-02 2.1 5.04e+06 2.5 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 1874
VecScatterBegin 32025 1.1 1.2784e-01 1.4 0.00e+00 0.0 7.6e+04 3.6e+02 0.0e+00 0 0 6 1 0 1 0 56 93 0 0
VecScatterEnd 32025 1.1 5.1963e-02 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
MatMult 24175 2.3 2.0385e-01 1.9 1.37e+08 2.4 2.6e+04 4.1e+02 0.0e+00 0 0 2 0 0 1 13 19 37 0 3361
MatSolve 24175 2.3 1.2673e+01 9.7 8.81e+08 2.5 3.1e+04 4.4e+02 6.6e+02 2 1 3 0 0 26 85 23 48 0 346
MatLUFactorSym 967 2.3 2.6745e-01 1.8 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 1 0 0 0 0 0
MatLUFactorNum 967 2.3 3.9942e-01 2.1 1.28e+07 2.5 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 2 1 0 0 0 159
PCSetUp 967 2.3 6.7788e-01 1.9 1.28e+07 2.5 0.0e+00 0.0e+00 2.6e+02 0 0 0 0 0 3 1 0 0 0 94
PCApply 24175 2.3 1.2693e+01 9.6 8.81e+08 2.5 3.1e+04 4.4e+02 6.6e+02 2 1 3 0 0 26 85 23 48 0 345
KSPSetUp 967 2.3 3.9855e-04 2.5 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
KSPSolve 24175 2.3 1.2738e+01 9.4 8.81e+08 2.5 3.1e+04 4.4e+02 3.9e+03 2 1 3 0 1 26 85 23 48 3 344
--- Event Stage 16: resetCMI
--- Event Stage 17: PEupdateCMI
BuildTwoSided 1958 1.0 1.9379e-02 1.3 0.00e+00 0.0 4.9e+02 4.0e+00 2.0e+03 0 0 0 0 0 3 0 4 0 34 0
BuildTwoSidedF 991 1.0 5.8963e-03 1.1 0.00e+00 0.0 0.0e+00 0.0e+00 9.9e+02 0 0 0 0 0 1 0 0 0 17 0
SFSetGraph 967 1.0 1.2741e-03 2.2 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFSetUp 967 1.0 1.5936e-02 1.1 0.00e+00 0.0 9.9e+02 2.3e+02 9.7e+02 0 0 0 0 0 3 0 7 2 17 0
SFPack 24175 1.0 1.6678e-03 1.1 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFUnpack 24175 1.0 1.0706e-03 1.1 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
VecScatterBegin 24175 1.0 2.6274e-02 1.4 0.00e+00 0.0 1.2e+04 9.3e+02 0.0e+00 0 0 1 0 0 4 0 93 98 0 0
VecScatterEnd 24175 1.0 8.6421e-03 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 1 0 0 0 0 0
MatAssemblyBegin 991 1.0 7.8690e-03 1.0 0.00e+00 0.0 0.0e+00 0.0e+00 9.9e+02 0 0 0 0 0 1 0 0 0 17 0
MatAssemblyEnd 991 1.0 8.8338e-04 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
--- Event Stage 18: updateCMI
BuildTwoSided 55 1.0 5.7917e-01 23.5 0.00e+00 0.0 4.6e+02 4.0e+00 5.5e+01 0 0 0 0 0 35 0 23 0 100 0
BuildTwoSidedF 55 1.0 5.7941e-01 23.2 0.00e+00 0.0 9.2e+02 4.3e+04 5.5e+01 0 0 0 1 0 35 0 46 19 100 0
SFPack 11 1.0 5.6954e-04 1.3 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFUnpack 11 1.0 2.2720e-06 1.8 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
VecScatterBegin 11 1.0 7.3393e-04 1.3 0.00e+00 0.0 2.2e+02 1.8e+04 0.0e+00 0 0 0 0 0 0 0 11 2 0 0
VecScatterEnd 11 1.0 1.1336e-02 44.8 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
MatMult 11 1.0 2.2682e-02 1.4 1.21e+07 2.5 2.2e+02 1.8e+04 0.0e+00 0 0 0 0 0 2 1 11 2 0 2651
MatAssemblyBegin 55 1.0 5.9361e-01 15.4 0.00e+00 0.0 9.2e+02 4.3e+04 5.5e+01 0 0 0 1 0 36 0 46 19 100 0
MatAssemblyEnd 55 1.0 2.6936e-02 1.7 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 2 0 0 0 0 0
MatTranspose 22 1.0 9.0290e-02 1.1 0.00e+00 0.0 9.2e+02 4.3e+04 2.2e+01 0 0 0 1 0 8 0 46 19 40 0
MatMatMultNum 11 1.0 2.6104e-01 1.0 2.01e+08 2.5 2.2e+02 3.1e+05 1.1e+01 0 0 0 2 0 24 21 11 33 20 3813
MatRARtNum 11 1.0 7.7935e-01 1.0 7.31e+08 2.5 1.1e+03 1.2e+05 3.3e+01 0 1 0 3 0 72 77 54 61 60 4649
MatGetLocalMat 33 1.0 1.0373e-01 2.8 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 6 0 0 0 0 0
MatGetBrAoCol 33 1.0 3.6744e-01 7.4 0.00e+00 0.0 8.6e+02 1.9e+05 0.0e+00 0 0 0 4 0 20 0 43 79 0 0
--- Event Stage 19: solveCMI
BuildTwoSided 14 1.0 2.8963e-03 3.6 0.00e+00 0.0 3.3e+02 4.0e+00 1.4e+01 0 0 0 0 0 0 0 0 0 0 0
SFSetGraph 14 1.0 1.8400e-03 3.3 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFSetUp 14 1.0 3.2496e-03 2.0 0.00e+00 0.0 6.5e+02 1.1e+03 1.4e+01 0 0 0 0 0 0 0 0 0 0 0
SFPack 8178 1.0 3.0207e-02 1.7 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFUnpack 8178 1.0 1.7554e-02 6.7 7.59e+04 1.5 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 28
VecDot 12042 1.0 4.0926e-01 1.6 5.83e+08 7.1 0.0e+00 0.0e+00 2.4e+04 0 0 0 0 4 1 2 0 0 88 3915
VecMDot 788 1.0 2.9636e-01 1.2 5.43e+08 7.1 0.0e+00 0.0e+00 2.2e+04 0 0 0 0 4 1 1 0 0 82 5038
VecNorm 820 1.0 1.1515e-01 14.6 3.97e+07 7.1 0.0e+00 0.0e+00 1.6e+03 0 0 0 0 0 0 0 0 0 6 948
VecScale 820 1.0 3.8266e-03 1.9 1.99e+07 7.1 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 14258
VecCopy 32 1.0 7.1528e-04 3.1 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
VecSet 2503 1.0 1.1924e-02 2.9 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
VecAXPY 12063 1.0 1.2886e-01 6.7 5.84e+08 7.1 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 2 0 0 0 12457
VecMAXPY 820 1.0 1.3051e-01 6.2 5.82e+08 7.1 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 2 0 0 0 12245
VecScatterBegin 8178 1.0 2.3742e-01 2.1 0.00e+00 0.0 1.4e+05 4.5e+03 0.0e+00 0 0 12 14 0 0 0 99 100 0 0
VecScatterEnd 8178 1.0 9.0039e-01 5.0 7.59e+04 1.5 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 2 0 0 0 0 1
VecNormalize 820 1.0 1.1844e-01 9.3 5.96e+07 7.1 0.0e+00 0.0e+00 1.6e+03 0 0 0 0 0 0 0 0 0 6 1382
MatMult 809 1.0 1.5082e+00 1.1 1.81e+09 4.6 8.3e+04 2.5e+03 1.6e+03 1 1 7 5 0 4 6 59 33 6 4248
MatMultAdd 3247 1.0 1.4236e+00 1.2 1.82e+09 4.5 8.3e+04 2.5e+03 0.0e+00 1 1 7 5 0 4 6 59 33 0 4545
MatMultTrAdd 11 1.0 1.3360e-01 17.3 1.22e+07 2.5 2.2e+02 1.8e+04 0.0e+00 0 0 0 0 0 0 0 0 1 0 455
MatSolve 1640 1.0 3.2236e+01 1.0 1.39e+10 1.4 5.7e+04 7.3e+03 6.0e+01 14 16 5 9 0 90 90 41 67 0 2873
MatLUFactorSym 2 1.0 3.5744e-01 1.0 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 1 0 0 0 0 0
MatLUFactorNum 12 1.0 1.1628e+00 1.0 2.25e+07 1.2 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 3 0 0 0 0 139
MatCreateSubMat 2 1.0 1.1445e-04 1.2 0.00e+00 0.0 0.0e+00 0.0e+00 4.0e+00 0 0 0 0 0 0 0 0 0 0 0
PCSetUp 11 1.0 2.6599e-01 1.0 9.30e+06 1.6 4.8e+01 1.1e+04 2.3e+01 0 0 0 0 0 1 0 0 0 0 228
PCSetUpOnBlocks 11 1.0 1.2561e+00 1.0 1.48e+07 1.4 0.0e+00 0.0e+00 4.0e+00 1 0 0 0 0 4 0 0 0 0 80
PCApply 820 1.0 3.2616e+01 1.0 1.39e+10 1.4 5.7e+04 7.3e+03 1.7e+03 14 16 5 9 0 91 90 41 67 6 2841
KSPSolve_FS_0 820 1.0 2.2055e+01 1.0 1.28e+10 1.4 3.4e+04 9.9e+03 8.2e+02 9 15 3 7 0 62 83 24 53 3 3870
KSPSolve_FS_1 820 1.0 1.0478e+01 1.0 1.12e+09 1.5 2.3e+04 3.6e+03 8.8e+02 4 1 2 2 0 29 7 16 13 3 699
KSPSetUp 1 1.0 7.4526e-04 6.0 0.00e+00 0.0 0.0e+00 0.0e+00 5.0e+00 0 0 0 0 0 0 0 0 0 0 0
KSPSolve 11 1.0 3.4378e+01 1.0 1.50e+10 1.3 1.4e+05 4.5e+03 2.7e+04 15 18 12 14 4 96 100 100 99 100 2977
KSPGMRESOrthog 788 1.0 3.8267e-01 1.2 1.09e+09 7.1 0.0e+00 0.0e+00 2.2e+04 0 1 0 0 4 1 3 0 0 82 7804
--- Event Stage 20: residualI
SFPack 25 1.0 6.0155e-04 1.4 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
SFUnpack 25 1.0 1.2346e-05 3.7 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 0 0 0 0 0 0
VecNorm 25 1.0 1.0071e-01 2.4 4.10e+06 2.6 0.0e+00 0.0e+00 2.5e+01 0 0 0 0 0 25 2 0 0 100 207
VecAYPX 25 1.0 2.0805e-03 2.0 2.05e+06 2.6 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 1 1 0 0 0 5007
VecScatterBegin 25 1.0 1.4292e-03 1.9 0.00e+00 0.0 9.5e+02 1.3e+04 0.0e+00 0 0 0 0 0 0 0 100 100 0 0
VecScatterEnd 25 1.0 5.4926e-02 19.5 0.00e+00 0.0 0.0e+00 0.0e+00 0.0e+00 0 0 0 0 0 7 0 0 0 0 0
MatMult 25 1.0 2.6402e-01 1.8 1.87e+08 2.5 9.5e+02 1.3e+04 0.0e+00 0 0 0 0 0 74 97 100 100 0 3611
from pathlib import Path
from IPython.display import display, Markdown
text = Path("proc_0_output.txt").read_text()
marker = "[MPI_MAX] Summary of timings"
idx = text.find(marker)
if idx != -1:
display(Markdown(f"```text\n{text[idx:]}\n```"))
else:
print(f"'{marker}' not found")
[MPI_MAX] Summary of timings | reps usr avg usr tot wall avg wall tot
-------------------------------------------------------------------------------------------------------------------------------
Build dofmap data | 24 0.005000 0.120000 0.004583 0.110000
Build sparsity | 1 0.010000 0.010000 0.010000 0.010000
Compute Node weighted graph partition (ParMETIS) | 1 0.040000 0.040000 0.030000 0.030000
Compute connectivity 0-0 | 3 0.000000 0.000000 0.000000 0.000000
Compute connectivity 1-0 | 2 0.000000 0.000000 0.000000 0.000000
Compute connectivity 2-0 | 4 0.005000 0.020000 0.002500 0.010000
Compute connectivity 3-0 | 1 0.000000 0.000000 0.000000 0.000000
Compute dof reordering map | 24 0.000833 0.020000 0.000833 0.020000
Compute entities of dim = 1 | 4 0.055000 0.220000 0.057500 0.230000
Compute entities of dim = 2 | 5 0.032000 0.160000 0.032000 0.160000
Compute entity permutations | 2 0.015000 0.030000 0.010000 0.020000
Compute graph partition (ParMETIS) | 2 0.020000 0.040000 0.015000 0.030000
Compute local part of mesh dual graph (mixed) | 9 0.018889 0.170000 0.024444 0.220000
Compute local-to-local map | 6 0.001667 0.010000 0.000000 0.000000
Compute non-local part of mesh dual graph | 3 0.006667 0.020000 0.006667 0.020000
Compute-local-to-global links for global/local adjacency list | 6 0.001667 0.010000 0.000000 0.000000
Distribute fixed-degree adjacency list to destination ranks | 3 0.006667 0.020000 0.006667 0.020000
Distribute row-wise data (scalable) | 6 0.001667 0.010000 0.000000 0.000000
Fine solve | 1 7.750000 7.750000 10.600000 10.600000
GPS: create_level_structure | 56 0.001143 0.040000 0.000571 0.020000
Gibbs-Poole-Stockmeyer ordering | 11 0.005714 0.040000 0.005714 0.040000
Init dofmap from element dofmap | 24 0.002500 0.060000 0.002917 0.070000
PLAZA: Enforce rules | 2 0.010000 0.020000 0.000000 0.000000
PLAZA: refine | 2 0.050000 0.100000 0.045000 0.090000
ParMETIS: call ParMETIS_V3_PartKway | 3 0.023333 0.070000 0.020000 0.060000
SparsityPattern::finalize | 12 0.035833 0.430000 0.036667 0.440000
TS solve | 1 41.100000 41.100000 42.630000 42.630000
TS solve Eimp | 1 54.740000 54.740000 57.480000 57.480000
Topology: create | 6 0.053333 0.320000 0.056667 0.340000
Topology: determine shared index ownership | 6 0.026667 0.160000 0.031667 0.190000
Topology: determine vertex ownership groups (owned, undetermined, unowned) | 6 0.003333 0.020000 0.003333 0.020000
~MPC: Apply lifting (C++) | 2 0.005000 0.010000 0.005000 0.010000
~MPC: Assemble matrix (C++) | 3 0.636667 1.910000 0.643333 1.930000
~MPC: Assemble vector (C++) | 3 0.000000 0.000000 0.000000 0.000000
~MPC: Create new index map with additional ghosts | 1 0.010000 0.010000 0.010000 0.010000
~MPC: Create sparsity pattern | 2 0.030000 0.060000 0.035000 0.070000
~MPC: Create sparsity pattern (Classic) | 2 0.010000 0.020000 0.010000 0.020000