Mon Jun 27 14:52:26 CDT 2016 job name = music job id = 7248487 using 1 nodes, 1 MPI tasks per node, 1 MPI tasks total using 32 OpenMP threads per MPI task using 32 CPUs total __ __ __ __ ______ __ ______ /\ "-./ \ /\ \/\ \ /\ ___\ /\ \ /\ ___\ \ \ \-./\ \ \ \ \_\ \ \ \___ \ \ \ \ \ \ \____ \ \_\ \ \_\ \ \_____\ \/\_____\ \ \_\ \ \_____\ \/_/ \/_/ \/_____/ \/_____/ \/_/ \/_____/ this is music! version 1.53 - Opening log file 'ic_agora_m11q.conf_log.txt'. - Using k-space sampled transfer functions... - Selecting transfer function plug-in 'eisenstein'... - starting at a=0.01 - Selecting region generator plug-in 'convex_hull'... - region point file appears to contain 3 columns - unapplying shift of previous zoom region to region particles : [-18,-10,6] = (-0.140625,-0.078125,0.046875) - Region center from convex hull centroid determined at (0.639073,0.577560,0.447254) - 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.619380,0.549677,0.433606] right = [0.660605,0.607909,0.463038] ------------------------------------------------------------- - Domain will be shifted by (-17, -10, 6) - Grid structure: Level 7 : offset = ( 0, 0, 0) size = ( 128, 128, 128) Level 8 : offset = ( 50, 48, 49) size = ( 60, 64, 58) Level 9 : offset = ( 12, 13, 13) size = ( 70, 78, 64) Level 10 : offset = ( 13, 12, 12) size = ( 90, 106, 78) Level 11 : offset = ( 13, 12, 13) size = ( 128, 162, 104) Level 12 : offset = ( 13, 12, 13) size = ( 206, 276, 158) Level 13 : offset = ( 12, 13, 12) size = ( 364, 502, 266) Level 14 : offset = ( 13, 12, 12) size = ( 678, 956, 482) ------------------------------------------------------------- - Finest level : extent = 2482.91 x 3500.98 x 1765.14 h-3 kpc**3 mtotgrid = 1.15822e+12 h-1 M_o particle mass = 3087.14 h-1 M_o baryon mass/cell = 620.154 h-1 M_o dx = 3.66211 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... - Generating random numbers w/ sample cube size of 256 - Generating a constrained random number set with seed -79280828 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... (232, 232, 232) - Performing noise convolution on level 11... - Performing density convolution... (344, 344, 344) - Performing noise convolution on level 12... - Performing density convolution... (568, 568, 568) - Performing noise convolution on level 13... - Performing density convolution... (1024, 1024, 1024) - Performing noise convolution on level 14... - Performing density convolution... (1928, 1928, 1928) - Density calculation took 426.858s with 32 threads. - Top grid mean density is off by -1.5194e-07, correcting... ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 8.54928e-07 - Poisson solver took 159.348s with 32 threads. - Computing 2LPT term.... - Solving 2LPT Poisson equation ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 9.1482e-07 - Poisson solver took 162.468s with 32 threads. - Performing hybrid Poisson step... (1912, 1912, 1912) - velocity component 0 : sigma = 0.0616284 - velocity component 1 : sigma = 0.0894673 - velocity component 2 : sigma = 0.0861595 ============================================================= 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... (232, 232, 232) - Performing noise convolution on level 11... - Performing density convolution... (344, 344, 344) - Performing noise convolution on level 12... - Performing density convolution... (568, 568, 568) - Performing noise convolution on level 13... - Performing density convolution... (1024, 1024, 1024) - Performing noise convolution on level 14... - Performing density convolution... (1928, 1928, 1928) - Density calculation took 590.606s with 32 threads. - Top grid mean density is off by -1.55502e-07, correcting... ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 8.80521e-07 - Poisson solver took 159.795s with 32 threads. - Writing baryon potential ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 6 steps to 8.60007e-08 - Poisson solver took 193.715s with 32 threads. - Performing hybrid Poisson step... (1912, 1912, 1912) - velocity component 0 : sigma = 0.0615937 - velocity component 1 : sigma = 0.0894547 - velocity component 2 : sigma = 0.0862657 ============================================================= 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... (232, 232, 232) - Performing noise convolution on level 11... - Performing density convolution... (344, 344, 344) - Performing noise convolution on level 12... - Performing density convolution... (568, 568, 568) - Performing noise convolution on level 13... - Performing density convolution... (1024, 1024, 1024) - Performing noise convolution on level 14... - Performing density convolution... (1928, 1928, 1928) - Density calculation took 444.108s with 32 threads. - Top grid mean density is off by -1.5194e-07, correcting... ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 7.48367e-07 - Poisson solver took 157.403s with 32 threads. ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 9.22689e-07 - Poisson solver took 145.118s with 32 threads. - Performing hybrid Poisson step... (1912, 1912, 1912) - max. x-displacement of HR particles is -7.723090 [mean dx] - max. y-displacement of HR particles is 14.621579 [mean dx] - max. z-displacement of HR particles is -13.082934 [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... (232, 232, 232) - Performing noise convolution on level 11... - Performing density convolution... (344, 344, 344) - Performing noise convolution on level 12... - Performing density convolution... (568, 568, 568) - Performing noise convolution on level 13... - Performing density convolution... (1024, 1024, 1024) - Performing noise convolution on level 14... - Performing density convolution... (1928, 1928, 1928) - Density calculation took 440.034s with 32 threads. - Top grid mean density is off by -1.55502e-07, correcting... ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 9.889e-07 - Poisson solver took 158.127s with 32 threads. ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 6 steps to 8.91291e-08 - Poisson solver took 173.162s with 32 threads. - Performing hybrid Poisson step... (1912, 1912, 1912) - Gadget2 : writing 268678472 particles to file... - type 0 : 129840464 [m=6.20154e-08] - type 1 : 129840464 [m=3.08714e-07] - type 2 : 8997544 [m=0] - Gadget2 : set initial gas temperature to 137.63 K/mu ============================================================= - Wrote output file 'ic_agora_m11q_ref14_rad5-chull.ics' using plugin 'gadget2'... - Done! wall time = 5748 sec = 0.07 day = 1:35:47 cpu time = 51.1 hr Mon Jun 27 16:28:14 CDT 2016