Tue May 12 23:41:17 CDT 2015 job name = music job id = 5252342 using 1 nodes, 1 MPI tasks, 1 MPI tasks per node using 32 OpenMP threads per MPI task __ __ __ __ ______ __ ______ /\ "-./ \ /\ \/\ \ /\ ___\ /\ \ /\ ___\ \ \ \-./\ \ \ \ \_\ \ \ \___ \ \ \ \ \ \ \____ \ \_\ \ \_\ \ \_____\ \/\_____\ \ \_\ \ \_____\ \/_/ \/_/ \/_____/ \/_____/ \/_/ \/_____/ this is music! version 1.51 - Opening log file 'ic_agora_m12i.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 - unapplying shift of previous zoom region to region particles : [-16,10,-14] = (-0.125000,0.078125,-0.109375) - Region center from convex hull centroid determined at (0.630174,0.405907,0.611845) - 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.593396,0.379671,0.557484] right = [0.664305,0.437849,0.652260] ------------------------------------------------------------- - Domain will be shifted by (-16, 11, -13) - Grid structure: Level 7 : offset = ( 0, 0, 0) size = ( 128, 128, 128) Level 8 : offset = ( 51, 50, 49) size = ( 54, 52, 62) Level 9 : offset = ( 9, 10, 10) size = ( 72, 66, 84) Level 10 : offset = ( 9, 10, 9) size = ( 106, 92, 130) Level 11 : offset = ( 10, 10, 10) size = ( 174, 146, 222) Level 12 : offset = ( 10, 9, 9) size = ( 310, 256, 408) Level 13 : offset = ( 10, 9, 9) size = ( 582, 476, 778) ------------------------------------------------------------- - Finest level : extent = 4262.7 x 3486.33 x 5698.24 h-3 kpc**3 mtotgrid = 6.39229e+12 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... (280, 280, 280) - Performing noise convolution on level 11... - Performing density convolution... (464, 464, 464) - Performing noise convolution on level 12... - Performing density convolution... (832, 832, 832) - Performing noise convolution on level 13... - Performing density convolution... (1576, 1576, 1576) - Density calculation took 244.022s with 32 threads. - Top grid mean density is off by 1.40706e-09, correcting... ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 3.90537e-07 - Poisson solver took 126.29s with 32 threads. - Computing 2LPT term.... - Solving 2LPT Poisson equation ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 6 steps to 1.76669e-07 - Poisson solver took 148.39s with 32 threads. - Performing hybrid Poisson step... (1556, 1556, 1556) - velocity component 0 : sigma = 0.127951 - velocity component 1 : sigma = 0.159753 - velocity component 2 : sigma = 0.20011 ============================================================= 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... (280, 280, 280) - Performing noise convolution on level 11... - Performing density convolution... (464, 464, 464) - Performing noise convolution on level 12... - Performing density convolution... (832, 832, 832) - Performing noise convolution on level 13... - Performing density convolution... (1576, 1576, 1576) - Density calculation took 224.498s with 32 threads. - Top grid mean density is off by 1.91738e-09, correcting... ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 4.91458e-07 - Poisson solver took 122.547s with 32 threads. - Writing baryon potential ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 8.20623e-07 - Poisson solver took 121.788s with 32 threads. - Performing hybrid Poisson step... (1556, 1556, 1556) - velocity component 0 : sigma = 0.127884 - velocity component 1 : sigma = 0.159715 - velocity component 2 : sigma = 0.200088 ============================================================= 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... (280, 280, 280) - Performing noise convolution on level 11... - Performing density convolution... (464, 464, 464) - Performing noise convolution on level 12... - Performing density convolution... (832, 832, 832) - Performing noise convolution on level 13... - Performing density convolution... (1576, 1576, 1576) - Density calculation took 220.671s with 32 threads. - Top grid mean density is off by 1.40706e-09, correcting... ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 3.94046e-07 - Poisson solver took 123.594s with 32 threads. ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 6 steps to 1.72147e-07 - Poisson solver took 147.156s with 32 threads. - Performing hybrid Poisson step... (1556, 1556, 1556) - max. x-displacement of HR particles is -8.395323 [mean dx] - max. y-displacement of HR particles is 10.657855 [mean dx] - max. z-displacement of HR particles is 13.258490 [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... (280, 280, 280) - Performing noise convolution on level 11... - Performing density convolution... (464, 464, 464) - Performing noise convolution on level 12... - Performing density convolution... (832, 832, 832) - Performing noise convolution on level 13... - Performing density convolution... (1576, 1576, 1576) - Density calculation took 219.795s with 32 threads. - Top grid mean density is off by 1.91738e-09, correcting... ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 4.83969e-07 - Poisson solver took 126.233s with 32 threads. ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 8.36883e-07 - Poisson solver took 125.619s with 32 threads. - Performing hybrid Poisson step... (1556, 1556, 1556) - Gadget2 : writing 146541875 particles to file... - type 0 : 70514272 [m=4.96123e-07] - type 1 : 70514272 [m=2.46971e-06] - type 2 : 5513331 [m=0] - Gadget2 : set initial gas temperature to 140.40 K/mu ============================================================= - Wrote output file 'ic_agora_m12i_ref13_rad4-chull.ics' using plugin 'gadget2'... - Done! wall time: 4123 sec = 0.05 day = 1:08:42 cpu time: 36.6 hr Wed May 13 00:50:00 CDT 2015