Tue Jan 19 20:24:15 CST 2016 job name = music job id = 6391860 using 1 nodes, 1 MPI tasks, 1 MPI tasks per node using 32 OpenMP threads per MPI task __ __ __ __ ______ __ ______ /\ "-./ \ /\ \/\ \ /\ ___\ /\ \ /\ ___\ \ \ \-./\ \ \ \ \_\ \ \ \___ \ \ \ \ \ \ \____ \ \_\ \ \_\ \ \_____\ \/\_____\ \ \_\ \ \_____\ \/_/ \/_/ \/_____/ \/_____/ \/_/ \/_____/ this is music! version 1.53 - Opening log file 'ic_agora_m12c.conf_log.txt'. - Using k-space sampled transfer functions... - Selecting transfer function plug-in 'eisenstein'... - starting at a=0.00990099 - Selecting region generator plug-in 'convex_hull'... - region point file appears to contain 3 columns - Region center from convex hull centroid determined at (0.486153,0.482990,0.592764) - Selecting random number generator plug-in 'MUSIC'... - WARNING: The selected transfer function does not support distinct amplitudes for baryon and DM fields! Perturbation amplitudes will be identical! - refinement region is 'convex_hull', w/ bounding box left = [0.447725,0.453497,0.522155] right = [0.527847,0.513229,0.672325] ------------------------------------------------------------- - Domain will be shifted by (1, 2, -12) - Grid structure: Level 7 : offset = ( 0, 0, 0) size = ( 128, 128, 128) Level 8 : offset = ( 46, 48, 43) size = ( 70, 64, 86) Level 9 : offset = ( 13, 12, 12) size = ( 88, 78, 124) Level 10 : offset = ( 12, 13, 12) size = ( 126, 106, 198) Level 11 : offset = ( 13, 13, 13) size = ( 202, 160, 346) Level 12 : offset = ( 12, 12, 13) size = ( 354, 270, 640) Level 13 : offset = ( 13, 13, 12) size = ( 658, 490, 1232) ------------------------------------------------------------- - Finest level : extent = 4819.34 x 3588.87 x 9023.44 h-3 kpc**3 mtotgrid = 1.17809e+13 h-1 M_o particle mass = 24697.1 h-1 M_o baryon mass/cell = 4961.23 h-1 M_o dx = 7.32422 h-1 kpc ------------------------------------------------------------- - Will use levelmin = 9 to compute density field... - Selecting output plug-in 'gadget2'... ============================================================= GENERATING WHITE NOISE ------------------------------------------------------------- - Found 0 density constraint(s) to be obeyed. - Generating random numbers (2) with seed 95064 - Generating random numbers w/ sample cube size of 256 - Generating random numbers w/ sample cube size of 256 - Generating a constrained random number set with seed 31415 using coarse mode replacement... - Generating random numbers w/ sample cube size of 256 - Generating a constrained random number set with seed 27183 using coarse mode replacement... - Generating random numbers w/ sample cube size of 256 - Generating a constrained random number set with seed -27018710 using coarse mode replacement... - Generating random numbers w/ sample cube size of 256 - Generating a constrained random number set with seed -85591679 using coarse mode replacement... - Generating random numbers w/ sample cube size of 256 - Generating a constrained random number set with seed -20707859 using coarse mode replacement... ============================================================= COMPUTING VELOCITIES ------------------------------------------------------------- - Using k-space transfer function kernel. - Performing noise convolution on level 9 - Performing density convolution... (512, 512, 512) - Performing noise convolution on level 10... - Performing density convolution... (416, 416, 416) - Performing noise convolution on level 11... - Performing density convolution... (712, 712, 712) - Performing noise convolution on level 12... - Performing density convolution... (1296, 1296, 1296) - Performing noise convolution on level 13... - Performing density convolution... (2480, 2480, 2480) - Density calculation took 967.615s with 32 threads. - Top grid mean density is off by -1.50434e-07, correcting... ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 4.78325e-07 - Poisson solver took 202.967s with 32 threads. - Computing 2LPT term.... - Solving 2LPT Poisson equation ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 7.43338e-07 - Poisson solver took 197.814s with 32 threads. - Performing hybrid Poisson step... (2464, 2464, 2464) - velocity component 0 : sigma = 0.116461 - velocity component 1 : sigma = 0.178583 - velocity component 2 : sigma = 0.205331 ============================================================= COMPUTING BARYON VELOCITIES ------------------------------------------------------------- - Using k-space transfer function kernel. - Performing noise convolution on level 9 - Performing density convolution... (512, 512, 512) - Performing noise convolution on level 10... - Performing density convolution... (416, 416, 416) - Performing noise convolution on level 11... - Performing density convolution... (712, 712, 712) - Performing noise convolution on level 12... - Performing density convolution... (1296, 1296, 1296) - Performing noise convolution on level 13... - Performing density convolution... (2480, 2480, 2480) - Density calculation took 960.069s with 32 threads. - Top grid mean density is off by -1.56464e-07, correcting... ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 6 steps to 6.19464e-08 - Poisson solver took 232.645s with 32 threads. - Writing baryon potential ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 6 steps to 1.0069e-07 - Poisson solver took 217.021s with 32 threads. - Performing hybrid Poisson step... (2464, 2464, 2464) - velocity component 0 : sigma = 0.116459 - velocity component 1 : sigma = 0.178658 - velocity component 2 : sigma = 0.205348 ============================================================= COMPUTING DARK MATTER DISPLACEMENTS ------------------------------------------------------------- - Using k-space transfer function kernel. - Performing noise convolution on level 9 - Performing density convolution... (512, 512, 512) - Performing noise convolution on level 10... - Performing density convolution... (416, 416, 416) - Performing noise convolution on level 11... - Performing density convolution... (712, 712, 712) - Performing noise convolution on level 12... - Performing density convolution... (1296, 1296, 1296) - Performing noise convolution on level 13... - Performing density convolution... (2480, 2480, 2480) - Density calculation took 983.075s with 32 threads. - Top grid mean density is off by -1.50434e-07, correcting... ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 4.85036e-07 - Poisson solver took 202.391s with 32 threads. ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 7.90258e-07 - Poisson solver took 182.966s with 32 threads. - Performing hybrid Poisson step... (2464, 2464, 2464) - max. x-displacement of HR particles is -14.492625 [mean dx] - max. y-displacement of HR particles is 17.575317 [mean dx] - max. z-displacement of HR particles is 14.655553 [mean dx] ============================================================= COMPUTING BARYON DISPLACEMENTS ------------------------------------------------------------- - Using k-space transfer function kernel. - Performing noise convolution on level 9 - Performing density convolution... (512, 512, 512) - Performing noise convolution on level 10... - Performing density convolution... (416, 416, 416) - Performing noise convolution on level 11... - Performing density convolution... (712, 712, 712) - Performing noise convolution on level 12... - Performing density convolution... (1296, 1296, 1296) - Performing noise convolution on level 13... - Performing density convolution... (2480, 2480, 2480) - Density calculation took 977.122s with 32 threads. - Top grid mean density is off by -1.56464e-07, correcting... ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 6 steps to 5.98929e-08 - Poisson solver took 234.019s with 32 threads. ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 6 steps to 1.05803e-07 - Poisson solver took 222.064s with 32 threads. - Performing hybrid Poisson step... (2464, 2464, 2464) - Gadget2 : writing 311621707 particles to file... - type 0 : 150682480 [m=4.96123e-07] - type 1 : 150682480 [m=2.46971e-06] - type 2 : 10256747 [m=0] - Gadget2 : set initial gas temperature to 140.40 K/mu ============================================================= - Wrote output file 'ic_agora_m12c_ref13_rad7-chull.ics' using plugin 'gadget2'... - Done! wall time: 12318 sec = 0.14 day = 3:25:18 cpu time: 109.5 hr Tue Jan 19 23:49:33 CST 2016