*************************************************************
			 S3P

ACE3P Codes Source Date:   Mon Mar 6 23:26:49 2023 -0800
ACE3P Codes Source Branch: master
ACE3P Codes Source Tag:    503652f41066a31de4a9b7088dd1fd7572ada43a
Support Lib Source Date:   Fri Dec 2 09:43:53 2022 -0800
Support Lib Source Branch: master
Support Lib Source Tag:    2f7bd8bf8ec6eb3646dc05e32622a4475531a105
Compilation Date:          Mon 06 Mar 2023 11:38:55 PM PST

*************************************************************

         Copyright 2023, Stanford University 

Authors make no representations or warranties, expressed or 
implied. By way of example, but not limitation, authors make 
no representatinos or warranties of merchantibility or 
fitness for any particular purpose or that the use of the 
software componentns of documentation will not infringe any
patents, copyrights, trademarks or other rights.

The Authors shall not be held liable for any liability 
nor for any direct, indirect or consequential damages with 
respect to any claim by users or any third party on account 
of or arising from the use of this software.

*************************************************************

Advanced Computations Department
SLAC National Accelerator Laboratory
https://slacportal.slac.stanford.edu/sites/ard_public/acd/Pages/Default.aspx
Contact: ace3p@slac.stanford.edu

Thank you for citing ACE3P when publishing related results.

*************************************************************

Starting master process on nid004422
Number of MPI processes: 128
Number of compute nodes: 1
Number of processes per node: 128
Data precision: 64 bits
Compiler: 11.2.0 20210728 (Cray Inc.)
Boundary conditions:
0 = INTERIOR
1 = IMPEDANCE
2 = MAGNETIC
3 = MAGNETIC
4 = WAVEGUIDE
5 = WAVEGUIDE

Material 1 has the following properties:
   Eps_r: (9, -0.0009)
    Mu_r: (1, 0)
   Sigma: 0

Boundary 1 has surface conductivity of 5.7e+07 S/m
Surface impedance model 
Metal properties: 
  Epsilon = 1
  MU = 1
  Sigma = 5.7e+07
  Frequency = 1e+09

Read Mesh: /pscratch/sd/l/liling/cw23/tem3p/SrfCoupler/2D-TTF3-FPC-Vacuum.ncdf
Time for reading the model: 0.0430083

Using curved quadratic tetrahedrons
Setting global vector finite element basis order to p=2
Partitioning Method: parmetis
      Setup: Max:   0.001, Sum:   0.174, Balance:   1.042
   Matching: Max:   0.003, Sum:   0.341, Balance:   1.001
Contraction: Max:   0.002, Sum:   0.208, Balance:   1.001
   InitPart: Max:   0.010, Sum:   1.301, Balance:   1.000
    Project: Max:   0.000, Sum:   0.007, Balance:   1.816
 Initialize: Max:   0.000, Sum:   0.057, Balance:   1.088
      K-way: Max:   0.006, Sum:   0.722, Balance:   1.038
      Remap: Max:   0.000, Sum:   0.013, Balance:   1.004
      Total: Max:   0.024, Sum:   3.022, Balance:   1.000

***********************************************************
*	Total Number of Elements read: 	139754
*	Total Number of Elements used: 	139754
*	Total Number of DOFs: 	968476
***********************************************************

Time for setting up finite element framework: 0.250318
/********************************/
/* input parameters, KVC syntax */
/********************************/


        Mesh : { 
            ReplicatedElementDistribution : { 
                average : 604.25 
                stddev : 301.48963497938 
                min : 85 
                max : 1571 
                total : 77344 
            }
            ElementDistribution : { 
                average : 1091.828125 
                stddev : 27.201814834951 
                min : 914 
                max : 1137 
                total : 139754 
            }
            File : /pscratch/sd/l/liling/cw23/tem3p/SrfCoupler/2D-TTF3-FPC-Vacuum.ncdf 
            MeshCoords : 32081 
        }
/********************************/





Checking Mesh Quality:

TETRAHEDRAL ELEMENTS:
  number   = 139754

INVERTED SECOND-ORDER ELEMENTS:
  number   = 0 <- GOOD!

ASPECT RATIO:
  min      = 1.02098
  max      = 21.8139 <- WARNING
  average  = 1.59691
  std dev  = 0.439985

SHAPE MEASURE:
  min      = 0.00141074 <- WARNING
  max      = 1.44152
  average  = 0.813522
  std dev  = 0.127058

ELEMENT VOLUME:
  min      = 2.83777e-12
  max      = 2.02203e-08
  average  = 6.54132e-10
  std dev  = 1.82823e-09

BOUNDING BOX:
  min      = (-0.5867, 0.00577425, -3.12445e-17)
  max      = (-0.025, 0.057, 0.021813)

EDGE LENGTH:
  min      = 0.0001
  max      = 0.0071349
  average  = 0.00122426
  std dev  = 0.0010002

Time for checking the mesh quality: 0.0187707
ESNumerical2D: Solving for port modes
 Numerical port solution time: 0.105597
ESNumerical2D: Solving for port modes
 Numerical port solution time: 0.070764
ESNumerical2D: Solving for port modes
 Numerical port solution time: 0.0542821
ESNumerical2D: Solving for port modes
 Numerical port solution time: 0.0391182
factorizing the matrix using MUMPS...
Use 32 processors to do parallel reordering (ParMetis)
Partition of Processors: 0 4 8 12 16 20 24 28 32 36 40 44 48 52 56 60 64 68 72 76 80 84 88 92 96 100 104 108 112 116 120 124 128 
Processor 0: 0, 30950
Processor 1: 30950, 61436
Processor 2: 61436, 92152
Processor 3: 92152, 122584
Processor 4: 122584, 153362
Processor 5: 153362, 182578
Processor 6: 182578, 213008
Processor 7: 213008, 243548
Processor 8: 243548, 273732
Processor 9: 273732, 303528
Processor 10: 303528, 333664
Processor 11: 333664, 363766
Processor 12: 363766, 394048
Processor 13: 394048, 424324
Processor 14: 424324, 454746
Processor 15: 454746, 484878
Processor 16: 484878, 514436
Processor 17: 514436, 543686
Processor 18: 543686, 573710
Processor 19: 573710, 603574
Processor 20: 603574, 634482
Processor 21: 634482, 664904
Processor 22: 664904, 694562
Processor 23: 694562, 725298
Processor 24: 725298, 755800
Processor 25: 755800, 786060
Processor 26: 786060, 816836
Processor 27: 816836, 847534
Processor 28: 847534, 877850
Processor 29: 877850, 907894
Processor 30: 907894, 937728
Processor 31: 937728, 968476
total: 968476 32
Generate ordering using parmetis...
Finished generating ordering using parmetis 


Analysis step: 1.962718 seconds

	Maximal per-core estimated memory 	 431 MB
	Aggregated estimated memory 	 45859 MB
	Maximal per-core estimated memory if OOC 	 351 MB
	Aggregated estimated memory if OOC 	 42281 MB


Factorization step: 1.304456 seconds

Memory usage:
used mem per MPI process: min:  222.1 MB, max:  340.9 MB, avg:  261.4 MB, stddev:  21.19 MB, total: 3.345e+04 MB
used mem per node in GB : min:  54.72 GB, max:  54.72 GB, avg:  54.72 GB, stddev:      0 GB, total:  54.72 GB
used mem per node in %  : min:  10.87  %, max:  10.87  %, avg:  10.87  %, stddev: 7.133e-15  %

S3P results for Frequency: 1.3e+09 Hz
**************************************

Index mapping: 
#          0 : Port 4, Mode 0, Type: TEM (cutoff: 0 Hz)
#          1 : Port 5, Mode 0, Type: TEM (cutoff: 0 Hz)
Frequency[Hz]           S(0,0)          S(0,1)          S(1,0)          S(1,1)
1.30000000e+09  9.83622408e-02  9.93681603e-01  9.93681603e-01  9.82094619e-02 
      Setup: Max:   0.001, Sum:   0.119, Balance:   1.076
   Matching: Max:   0.002, Sum:   0.272, Balance:   1.001
Contraction: Max:   0.001, Sum:   0.099, Balance:   1.002
   InitPart: Max:   0.007, Sum:   0.928, Balance:   1.000
    Project: Max:   0.000, Sum:   0.007, Balance:   1.872
 Initialize: Max:   0.000, Sum:   0.055, Balance:   1.087
      K-way: Max:   0.005, Sum:   0.641, Balance:   1.000
      Remap: Max:   0.000, Sum:   0.013, Balance:   1.005
      Total: Max:   0.017, Sum:   2.158, Balance:   1.000
Computed Total Energy (normalized by Epsilon0/2): (474.38111,-0.0016598072).  Q from energy: 142902.48
Writing mode file s3p_results/mode.port4.mode0.f1.300000000E+09.mod
Computed Total Energy (normalized by Epsilon0/2): (478.13778,-0.0019530658).  Q from energy: 122406.98


Analysis step: 0.002731 seconds

	Maximal per-core estimated memory 	 1 MB
	Aggregated estimated memory 	 1 MB
	Maximal per-core estimated memory if OOC 	 1 MB
	Aggregated estimated memory if OOC 	 1 MB


Factorization step: 0.027342 seconds



Analysis step: 0.000316 seconds

	Maximal per-core estimated memory 	 1 MB
	Aggregated estimated memory 	 1 MB
	Maximal per-core estimated memory if OOC 	 1 MB
	Aggregated estimated memory if OOC 	 1 MB


Factorization step: 0.000445 seconds



Analysis step: 0.000405 seconds

	Maximal per-core estimated memory 	 0 MB
	Aggregated estimated memory 	 0 MB
	Maximal per-core estimated memory if OOC 	 0 MB
	Aggregated estimated memory if OOC 	 0 MB


Factorization step: 0.001010 seconds



Analysis step: 0.000182 seconds

	Maximal per-core estimated memory 	 0 MB
	Aggregated estimated memory 	 0 MB
	Maximal per-core estimated memory if OOC 	 0 MB
	Aggregated estimated memory if OOC 	 0 MB


Factorization step: 0.000271 seconds

*************************************************************
			 Tem3P

ACE3P Codes Source Date:   Mon Mar 6 23:26:49 2023 -0800
ACE3P Codes Source Branch: master
ACE3P Codes Source Tag:    503652f41066a31de4a9b7088dd1fd7572ada43a
Support Lib Source Date:   Fri Dec 2 09:43:53 2022 -0800
Support Lib Source Branch: master
Support Lib Source Tag:    2f7bd8bf8ec6eb3646dc05e32622a4475531a105
Compilation Date:          Mon 06 Mar 2023 11:38:55 PM PST

*************************************************************

         Copyright 2023, Stanford University 

Authors make no representations or warranties, expressed or 
implied. By way of example, but not limitation, authors make 
no representatinos or warranties of merchantibility or 
fitness for any particular purpose or that the use of the 
software componentns of documentation will not infringe any
patents, copyrights, trademarks or other rights.

The Authors shall not be held liable for any liability 
nor for any direct, indirect or consequential damages with 
respect to any claim by users or any third party on account 
of or arising from the use of this software.

*************************************************************

Advanced Computations Department
SLAC National Accelerator Laboratory
https://slacportal.slac.stanford.edu/sites/ard_public/acd/Pages/Default.aspx
Contact: ace3p@slac.stanford.edu

Thank you for citing ACE3P when publishing related results.

*************************************************************

Starting master process on nid004422
Number of MPI processes: 128
Number of compute nodes: 1
Number of processes per node: 128
Data precision: 64 bits
Compiler: 11.2.0 20210728 (Cray Inc.)
Beginning run at Sat Mar 11 13:51:39 2023

---------------------------------

Entering ThermoElasticProblem::solve()...

---------------------------------

In ThermoElasticProblem::solve(): begining tem3p iter 0 with 'runThermostaticProblem' = 1 with 'runThermoTransientProblem' = 0 with 'runElasticProblem' = 0
1 In StaticProblem::readKvcParameters1() of ThermostaticProblem, rank 0: instantiating ScalarDiscreteDomain...
2 In StaticProblem::readKvcParameters1() of ThermostaticProblem, rank 0: setting up ScalarDiscreteDomain...
Read Mesh: 2D-TTF3-FPC-Body.ncdf
Check Mesh: 2D-TTF3-FPC-Body.ncdf
Time for read model: 0.04120367299992722
Time for setBasisOrder: 0.001406965000114724
    Time for making right-orientation and calculating Adjacency: 0.09187002299995584
Partitioning Method: PARMETIS
      Setup: Max:   0.002, Sum:   0.187, Balance:   1.041
   Matching: Max:   0.003, Sum:   0.341, Balance:   1.000
Contraction: Max:   0.002, Sum:   0.205, Balance:   1.001
   InitPart: Max:   0.009, Sum:   1.145, Balance:   1.000
    Project: Max:   0.000, Sum:   0.009, Balance:   1.879
 Initialize: Max:   0.001, Sum:   0.060, Balance:   1.095
      K-way: Max:   0.006, Sum:   0.767, Balance:   1.017
      Remap: Max:   0.000, Sum:   0.018, Balance:   1.004
      Total: Max:   0.022, Sum:   2.864, Balance:   1.000
    Time for partitioning mesh: 0.08737734100031957
    Time for replicating remote meshes: 0.02352289999998902
        Time for numberingHierarchical: 0.0483747560001575

***********************************************************
*	Total Number of Elements read: 	315091
*	Total Number of Elements used: 	315091
*	Total Number of DOFs: 	101857
***********************************************************

        Time for numbering remote nodes: 0.005346340999949462
        Time for numbering remote edges: 0.01691337000011117
        Time for numbering remote faces: 0.02491872500013415
        Time for numbering remote volumes: 0.0003098279998994258
    Time for numbering DOF entities: 0.0959137190002366
Time for set model: 0.2996180279997134
Time for read midpoints: 0.06859559200029253
Time for creating higher-order elements: 0.001468224999825907
Problem Name: ThermostaticProblemShell number         1
   Shell boundary id 6
   Shell basis order 1
   Shell layer       1
   Shell thickness   1e-05
   Shell material    221
Shell number         2
   Shell boundary id 7
   Shell basis order 1
   Shell layer       1
   Shell thickness   0.0001
   Shell material    221
 loading RF solution for TEM3P...
Omega3P directory name: s3p_results
Boundary conditions:
0 = INTERIOR
1 = IMPEDANCE
2 = MAGNETIC
3 = MAGNETIC
4 = WAVEGUIDE
5 = WAVEGUIDE

Material 1 has the following properties:
   Eps_r: (9, -0.0009)
    Mu_r: (1, 0)
   Sigma: 0

Boundary 1 has surface conductivity of 57000000.00000001 S/m
Surface impedance model 
Metal properties: 
  Epsilon = 1
  MU = 1
  Sigma = 57000000
  Frequency = 1000000000

Read Mesh: /pscratch/sd/l/liling/cw23/tem3p/SrfCoupler/2D-TTF3-FPC-Vacuum.ncdf
Time for reading the model: 0.3262526749999779

Using curved quadratic tetrahedrons
Setting global vector finite element basis order to p=2
Partitioning Method: parmetis

***********************************************************
*	Total Number of Elements read: 	139754
*	Total Number of Elements used: 	139754
*	Total Number of DOFs: 	968476
***********************************************************

Time for setting up finite element framework: 3.378400882999813
/********************************/
/* input parameters, KVC syntax */
/********************************/


        Mesh : { 
            ReplicatedElementDistribution : { 
                stddev : 0.0 
                total : 0 
                min : 0 
                max : 0 
                average : 0.0 
            }
            File : /pscratch/sd/l/liling/cw23/tem3p/SrfCoupler/2D-TTF3-FPC-Vacuum.ncdf 
            ElementDistribution : { 
                stddev : 0.0 
                total : 139754 
                min : 139754 
                max : 139754 
                average : 139754.0 
            }
            MeshCoords : 32081 
        }
/********************************/




In StaticProblem::readKvcParameters1() of ThermostaticProblem, rank 0: setting map of boundaries...
 EM frequency : (1300000000,0)
 Powerinput : 437.5 :  : 437.5
The scaling factor is 437.5
Volumetric Heating Boundary 5 dielectri constants de(T) = 
 EM frequency : (1300000000,0)
 Powerinput : 437.5 :  : 437.5
The scaling factor is 437.5
EM Heating Boundary 6 surface resis R(T) = 
 temp from -1000 to 4 K(T) = 76.67*sqrt(1.55e-8/20)
 temp from 4 to 300 K(T) = 76.67*sqrt(2.821e-8*exp(-190.0/T)+1.55e-8/20)
 temp from 300 to 1000 K(T) = 76.67*sqrt(1.57494e-8)
 temp from 300 to 1000 K(T) = 76.67*sqrt(1.57494e-8)
 EM frequency : (1300000000,0)
 Powerinput : 437.5 :  : 437.5
The scaling factor is 437.5
EM Heating Boundary 7 surface resis R(T) = 
 temp from -1000 to 4 K(T) = 76.67*sqrt(1.55e-8/20)
 temp from 4 to 300 K(T) = 76.67*sqrt(2.821e-8*exp(-190.0/T)+1.55e-8/20)
 temp from 300 to 1000 K(T) = 76.67*sqrt(1.57494e-8)
 temp from 300 to 1000 K(T) = 76.67*sqrt(1.57494e-8)
 EM frequency : (1300000000,0)
 Powerinput : 437.5 :  : 437.5
The scaling factor is 437.5
EM Heating Boundary 8 surface resis R(T) = 
 temp from -1000 to 4 K(T) = 76.67*sqrt(1.55e-8/20)
 temp from 4 to 300 K(T) = 76.67*sqrt(2.821e-8*exp(-190.0/T)+1.55e-8/20)
 temp from 300 to 1000 K(T) = 76.67*sqrt(1.57494e-8)
 temp from 300 to 1000 K(T) = 76.67*sqrt(1.57494e-8)
 EM frequency : (1300000000,0)
 Powerinput : 437.5 :  : 437.5
The scaling factor is 437.5
EM Heating Boundary 9 surface resis R(T) = 
 temp from -1000 to 2.1 K(T) = 3.5276e-7
 temp from 2.1 to 4.4 K(T) = 1.6359e-6*T-2.9191e-6
 temp from 4.4 to 1000 K(T) = 4.279e-6
 temp from 4.4 to 1000 K(T) = 4.279e-6
In StaticProblem::readKvcParameters1() of ThermostaticProblem, rank 0: map of boundaries set succesfully
In StaticProblem::readKvcParameters1() of ThermostaticProblem, rank 0: setting linear solver...
Leaving StaticProblem::readKvcParameters1() of ThermostaticProblem
 Volume id: 221
 temp from -10000 to 2 K(T) = 350
 temp from 2 to 300 K(T) = 10^((2.2154+(-0.88068)*T^0.5+0.29505*T^1.0+(-0.048310)*T^1.5+0.0032070*T^2.0)/(1+(-0.47461)*T^0.5+0.13871*T^1.0+(-0.02043)*T^1.5+0.0012810*T^2.0))
 temp from 300 to 10000 K(T) = 396.3
 temp from 300 to 10000 K(T) = 396.3
 Volume id: 1
 temp from -100 to 2 K(T) = 10.
 temp from 2 to 3 K(T) = 10*T-10
 temp from 3 to 4 K(T) = 25*T-55.
 temp from 4 to 5 K(T) = 55*T-175
 temp from 5 to 6 K(T) = 50*T-150
 temp from 6 to 7 K(T) = 70*T-270
 temp from 7 to 8 K(T) = 50*T-130
 temp from 8 to 9 K(T) = 30*T+30
 temp from 9 to 3000 K(T) = 300.
 temp from 9 to 3000 K(T) = 300.
 Volume id: 2
 temp from -1000 to 2 K(T) = 0.06
 temp from 2 to 3 K(T) = 0.04*T-0.02
 temp from 3 to 4 K(T) = 0.1*T-0.2
 temp from 4 to 5 K(T) = 0.08*T-0.12
 temp from 5 to 6 K(T) = 0.07*T-0.07
 temp from 6 to 7 K(T) = 0.07*T-0.07
 temp from 7 to 8 K(T) = 0.08*T-0.14
 temp from 8 to 9 K(T) = 0.1*T-0.3
 temp from 9 to 10 K(T) = 0.05*T+0.15
 temp from 10 to 20 K(T) = 0.075*T-0.1
 temp from 20 to 30 K(T) = 0.06*T+0.2
 temp from 30 to 40 K(T) = 0.05*T+0.5
 temp from 40 to 50 K(T) = 0.05*T+0.5
 temp from 50 to 60 K(T) = 0.1*T-2.0
 temp from 60 to 70 K(T) = 0.05*T+1
 temp from 70 to 80 K(T) = 0.05*T+1
 temp from 80 to 90 K(T) = 0.07*T-0.6
 temp from 90 to 100 K(T) = 0.03*T+3.0
 temp from 100 to 200 K(T) = 0.03*T+3.0
 temp from 200 to 300 K(T) = 0.008*T+7.4
 temp from 300 to 10000 K(T) = 9.8
 temp from 300 to 10000 K(T) = 9.8
 Volume id: 3
 temp from -10000 to 2.1 K(T) = 0.1124
 temp from 2.1 to 300 K(T) = 10^((-1.4087)+(1.3982)*log10(T)+(0.2543)*(log10(T))^2+(-0.6260)*(log10(T))^3+(0.2334)*(log10(T))^4+(0.4256)*(log10(T))^5+(-0.4658)*(log10(T))^6+(0.165)*(log10(T))^7+(-0.0199)*(log10(T))^8)
 temp from 300 to 10000 K(T) = 15.30865
 temp from 300 to 10000 K(T) = 15.30865
 Volume id: 4
 temp from -10000 to 2 K(T) = 350
 temp from 2 to 300 K(T) = 10^((2.2154+(-0.88068)*T^0.5+0.29505*T^1.0+(-0.048310)*T^1.5+0.0032070*T^2.0)/(1+(-0.47461)*T^0.5+0.13871*T^1.0+(-0.02043)*T^1.5+0.0012810*T^2.0))
 temp from 300 to 10000 K(T) = 396.3
 temp from 300 to 10000 K(T) = 396.3
 Volume id: 5
 temp from -1000 to 2 K(T) = 0.1
 temp from 2 to 300 K(T) = 10^(-2.26487+5.73102*log10(T)-6.7601*(log10(T))^2+5.96778*(log10(T))^3-1.75221*(log10(T))^4-0.171282*(log10(T))^5-0.0177992*(log10(T))^6+0.0856919*(log10(T))^7-0.0159493*(log10(T))^8)
 temp from 300 to 10000 K(T) = 36
 temp from 300 to 10000 K(T) = 36
Entering StaticProblem::instantiateMatrix() of ThermostaticProblem
Entering assembleMatrix().
 Time for assembling thermal matrix: 0.07940910700017412

Entering ThermostaticProblem::solveWithNewton() for global iteration 0
Entering assembleMatrix().
 Time for assembling thermal matrix: 0.1049824970000373
 ThermostaticProblem::computeRHS 
 Time for assembling RHS: 77.4601

 Initial norm         : 2941.062976517171
 DOFS norm            : 4138.193204769444
 Picard limit         : 100
 Max Nonlinear Iter   : 500
 Absolute Tol         : 1e-09
 Relative Tol         : 1e-09

 Nonlinear Newton iteration 0
Entering assembleMatrix().
 Time for assembling thermal matrix: 0.1093818939998528
 ThermostaticProblem::computeRHS 
 Time for assembling RHS: 76.8506

 Newton residual norm : 12.80728725514524
 Initial norm         : 2941.062976517171
 DOFS norm            : 102091.7014813783

 Nonlinear Newton iteration 1
Entering assembleMatrix().
 Time for assembling thermal matrix: 0.1025023980000697
 ThermostaticProblem::computeRHS 
 Time for assembling RHS: 77.334

 Newton residual norm : 0.5354560505086164
 Initial norm         : 2941.062976517171
 DOFS norm            : 87059.88012526993

 Nonlinear Newton iteration 2
Entering assembleMatrix().
 Time for assembling thermal matrix: 0.102281266999853
 ThermostaticProblem::computeRHS 
 Time for assembling RHS: 76.8021

 Newton residual norm : 0.01844334888494861
 Initial norm         : 2941.062976517171
 DOFS norm            : 86461.01621543053

 Nonlinear Newton iteration 3
Entering assembleMatrix().
 Time for assembling thermal matrix: 0.1021798719998515
 ThermostaticProblem::computeRHS 
 Time for assembling RHS: 76.9991

 Newton residual norm : 0.003293372469854504
 Initial norm         : 2941.062976517171
 DOFS norm            : 86261.50081796799

 Nonlinear Newton iteration 4
Entering assembleMatrix().
 Time for assembling thermal matrix: 0.1020625679998375
 ThermostaticProblem::computeRHS 
 Time for assembling RHS: 77.1499

 Newton residual norm : 0.0006122174728530146
 Initial norm         : 2941.062976517171
 DOFS norm            : 86224.59992749234

 Nonlinear Newton iteration 5
Entering assembleMatrix().
 Time for assembling thermal matrix: 0.1018799940002282
 ThermostaticProblem::computeRHS 
 Time for assembling RHS: 77.2445

 Newton residual norm : 0.0001051732992730229
 Initial norm         : 2941.062976517171
 DOFS norm            : 86219.06944617006

 Nonlinear Newton iteration 6
Entering assembleMatrix().
 Time for assembling thermal matrix: 0.1022953290002988
 ThermostaticProblem::computeRHS 
 Time for assembling RHS: 77.0815

 Newton residual norm : 7.905809985234996e-06
 Initial norm         : 2941.062976517171
 DOFS norm            : 86219.14645450425

 Nonlinear Newton iteration 7
Entering assembleMatrix().
 Time for assembling thermal matrix: 0.1018686960001105
 ThermostaticProblem::computeRHS 
 Time for assembling RHS: 77.2728

 Newton residual norm : 3.299021151355816e-06
 Initial norm         : 2941.062976517171
 DOFS norm            : 86219.18121534205

 Nonlinear Newton iteration 8
Entering assembleMatrix().
 Time for assembling thermal matrix: 0.103043555000113
 ThermostaticProblem::computeRHS 
 Time for assembling RHS: 77.2272

 Newton residual norm : 5.887218503899363e-07
 Initial norm         : 2941.062976517171
 DOFS norm            : 86219.15635929929

Leaving ThermostaticProblem::solveWithNewton() for global iteration 0
Distribute remote flux called 
Calculating the total power flow through every sideset...

 ===================================== 
 Max T = 324.3520007076649
 Min T = 2
 P (SideSet 1) = 0.08725785613942286
 P (SideSet 2) = 0.0005649793014200234
 P (SideSet 3) = 0.03957574428881457
 P (SideSet 4) = 0.7813392638376635
 P (SideSet 5) = 0.1005412253626891
 P (SideSet 6) = 0
 P (SideSet 7) = 0
 P (SideSet 8) = -0.03257071919747352
 P (SideSet 9) = -0.002889115633092689
 P (SideSet 10) = -0.1809940299856579
 P (SideSet 11) = -0.2004387158296005
 ===================================== 

 Temperature output file : tem3p_results/Temperature.mod


---------------------------------

Leaving ThermoElasticProblem::solve()

---------------------------------

Leaving StaticProblem::destructor() of ThermostaticProblem
Final result: total run time from timer manager 'Total run' is 829.599 seconds
Ending run at Sat Mar 11 14:05:29 2023
Total run time = 830 seconds
