************************************************************* Omega3P 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 nid004330 Number of MPI processes: 16 Number of compute nodes: 1 Number of processes per node: 16 Data precision: 64 bits Compiler: 11.2.0 20210728 (Cray Inc.) Boundary conditions: 0 = INTERIOR 3 = EXTERIOR 4 = EXTERIOR 6 = EXTERIOR 7 = WAVEGUIDE 8 = WAVEGUIDE Boundary 6 has surface conductivity of 58000000 S/m Boundary 3 has surface conductivity of 25000000 S/m Boundary 4 has surface conductivity of 25000000 S/m Surface impedance model Metal properties: Epsilon = 1 MU = 1 Sigma = 58000000 Frequency = 1000000000 Surface impedance model Metal properties: Epsilon = 1 MU = 1 Sigma = 25000000 Frequency = 1000000000 Surface impedance model Metal properties: Epsilon = 1 MU = 1 Sigma = 25000000 Frequency = 1000000000 Read Mesh: /pscratch/sd/l/lge/cw23/examples/track3p/Muon201MHz/Muon.ncdf Time for reading the model: 0.09019938400001593 Using curved quadratic tetrahedrons Setting global vector finite element basis order to p=2 Partitioning Method: parmetis *********************************************************** * Total Number of Elements read: 207642 * Total Number of Elements used: 207642 * Total Number of DOFs: 1247482 *********************************************************** Time for setting up finite element framework: 0.3744195500000842 /********************************/ /* input parameters, KVC syntax */ /********************************/ Mesh : { MeshCoords : 41889 ElementDistribution : { stddev : 834.31559776062 total : 207642 average : 12977.625 min : 10875 max : 13549 } ReplicatedElementDistribution : { stddev : 1054.7222521435 total : 42251 average : 2640.6875 min : 704 max : 4405 } File : /pscratch/sd/l/lge/cw23/examples/track3p/Muon201MHz/Muon.ncdf } /********************************/ Checking Mesh Quality: TETRAHEDRAL ELEMENTS: number = 207642 INVERTED SECOND-ORDER ELEMENTS: number = 0 <- GOOD! ASPECT RATIO: min = 1.034209861834297 max = 6.697910803072359 average = 1.62123524992641 std dev = 0.2320653641298694 SHAPE MEASURE: min = 0.04414136048337203 max = 1.174909388399401 average = 0.7968964944353626 std dev = 0.113471712772405 ELEMENT VOLUME: min = 3.989060984263698e-10 max = 2.509613025947056e-05 average = 1.940466344412714e-06 std dev = 3.168249239435707e-06 BOUNDING BOX: min = (-0.6277102851844776, -1.135710285184472, -0.2076001051193858) max = (0.6277102851844983, 1.135710285184472, 0.2365920719214952) EDGE LENGTH: min = 0.001125211976716972 max = 0.07781505221802698 average = 0.01959344394022806 std dev = 0.01412502299302757 Time for checking the mesh quality: 0.07288602699986768 Time for save/load ComputationalMesh: 0.00128531599989401 Total Volume of the structure is : 0.4029223126865449 Calling complex solver ESNumerical2D: Solving for port modes Numerical port solution time: 2.388147253999932 ESNumerical2D: Solving for port modes Numerical port solution time: 2.519377071000008 Reducing the number of unique cutoff frequencies from 2 to 1 (this can be done since they are closer than 0.1% apart from each other This allows the use of the SOAR solver for quadratic eigenvalue systems Starting SOAR to search for 1 eigenpairs. /********************************/ /* input parameters, KVC syntax */ /********************************/ FactorizationType : 1 Tolerance : 1e-08 MaxRestarts : 10 MaxBadIterations : 200 EigenvalueLocation : AboveShift ncv : 198 MaxIterations : 200 /********************************/ ========================================================== SOAR Parameters: Shift is (3.143767532927523,0) Search Space is 198 ========================================================== Preparing for solving linear system ... factorizing the matrix using MUMPS... Use 16 processors to do parallel reordering (ParMetis) Partition of Processors: 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 Processor 0: 0, 72212 Processor 1: 72212, 152430 Processor 2: 152430, 233250 Processor 3: 233250, 309530 Processor 4: 309530, 389064 Processor 5: 389064, 469072 Processor 6: 469072, 549698 Processor 7: 549698, 628302 Processor 8: 628302, 695700 Processor 9: 695700, 776780 Processor 10: 776780, 857832 Processor 11: 857832, 937928 Processor 12: 937928, 1018970 Processor 13: 1018970, 1099178 Processor 14: 1099178, 1165206 Processor 15: 1165206, 1247482 total: 1247482 16 Generate ordering using parmetis... Finished generating ordering using parmetis Memory usage: used mem per MPI process: min: 1988 MB, max: 3835 MB, avg: 2376 MB, stddev: 588.9 MB, total: 3.801e+04 MB used mem per node in GB : min: 72.78 GB, max: 72.78 GB, avg: 72.78 GB, stddev: 0 GB, total: 72.78 GB used mem per node in % : min: 14.46 %, max: 14.46 %, avg: 14.46 %, stddev: 1.835e-15 % ---------------------------------------------------------- Iteration: 0 Solving Projected System... The 0th ritz pair: Frequency: (17212442281.93695,50692.12007301897) Residual: 12049.28577787428 norm(v): 41.19694716399505 tolerance: required: 1e-08 reached: 292.480064843373 ---------------------------------------------------------- Iteration: 1 Solving Projected System... The 0th ritz pair: Frequency: (267888542.4130364,45011240.52650037) Residual: 31.67158296559618 norm(v): 20.4743113741728 tolerance: required: 1e-08 reached: 1.546893684812673 ---------------------------------------------------------- Iteration: 2 Solving Projected System... The 0th ritz pair: Frequency: (558944452.0496403,21916061.09047264) Residual: 75.52536742469175 norm(v): 14.38275177943631 tolerance: required: 1e-08 reached: 5.251106921880825 ---------------------------------------------------------- Iteration: 3 Solving Projected System... The 0th ritz pair: Frequency: (214388836.790329,8615210.650737597) Residual: 4.616687696842686 norm(v): 14.03148442863454 tolerance: required: 1e-08 reached: 0.329023470062886 ---------------------------------------------------------- Iteration: 4 Solving Projected System... The 0th ritz pair: Frequency: (202308352.0723563,42368.48282094699) Residual: 0.8774007150230536 norm(v): 10.61340235960618 tolerance: required: 1e-08 reached: 0.08266912770238274 ---------------------------------------------------------- Iteration: 5 Solving Projected System... The 0th ritz pair: Frequency: (201909491.9836876,734.1162158851755) Residual: 0.1057151904208474 norm(v): 10.59684220435239 tolerance: required: 1e-08 reached: 0.009976103105265408 ---------------------------------------------------------- Iteration: 6 Solving Projected System... The 0th ritz pair: Frequency: (201900860.5408965,2098.870125692776) Residual: 0.02236222444618413 norm(v): 10.59774887159466 tolerance: required: 1e-08 reached: 0.002110091937177528 ---------------------------------------------------------- Iteration: 7 Solving Projected System... The 0th ritz pair: Frequency: (201901016.7995947,1957.746355543227) Residual: 0.001499290392870216 norm(v): 10.59758919883792 tolerance: required: 1e-08 reached: 0.0001414746660527869 ---------------------------------------------------------- Iteration: 8 Solving Projected System... The 0th ritz pair: Frequency: (201901011.6284615,1975.459675415349) Residual: 0.0002335733474021258 norm(v): 10.59760192545578 tolerance: required: 1e-08 reached: 2.204020768520047e-05 ---------------------------------------------------------- Iteration: 9 Solving Projected System... Memory usage: used mem per MPI process: min: 1990 MB, max: 3838 MB, avg: 2379 MB, stddev: 588.7 MB, total: 3.806e+04 MB used mem per node in GB : min: 72.81 GB, max: 72.81 GB, avg: 72.81 GB, stddev: 0 GB, total: 72.81 GB used mem per node in % : min: 14.47 %, max: 14.47 %, avg: 14.47 %, stddev: 5.504e-15 % The 0th ritz pair: Frequency: (201901009.5454043,1970.122189028231) Residual: 3.705026707977634e-05 norm(v): 10.59759946897525 tolerance: required: 1e-08 reached: 3.49609996001849e-06 ---------------------------------------------------------- Iteration: 10 Solving Projected System... The 0th ritz pair: Frequency: (201901010.550669,1970.961625659273) Residual: 2.49593058470018e-06 norm(v): 10.59759953950982 tolerance: required: 1e-08 reached: 2.355184846714473e-07 ---------------------------------------------------------- Iteration: 11 Solving Projected System... The 0th ritz pair: Frequency: (201901010.5462654,1971.002863172307) Residual: 2.227560712682635e-07 norm(v): 10.59759956467093 tolerance: required: 1e-08 reached: 2.101948369618175e-08 ---------------------------------------------------------- Iteration: 12 Solving Projected System... The 0th ritz pair: Frequency: (201901010.5458142,1970.998041451227) Residual: 1.804020304509099e-08 norm(v): 10.59759956240736 tolerance: required: 1e-08 reached: 1.702291442402167e-09 Linear Solver Preparation Time: 27.55956726281252 Solver Time: 0 ========================================================== SOAR Adaptive found 1 eigenvalues COMMIT MODE: 0 FREQ = (201901010.5458142,1970.998041451227) k = (4.231532278794588,4.130906433452068e-05) Q = 51217.96326016552