Tue Jan 19 20:21:01 CST 2016 job name = music job id = 6391845 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_m12m.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 : [0,14,-12] = (0.000000,0.109375,-0.093750) - Region center from convex hull centroid determined at (0.497590,0.383508,0.586124) - 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.458308,0.358734,0.520541] right = [0.538674,0.406962,0.657671] ------------------------------------------------------------- - Domain will be shifted by (0, 14, -11) - Grid structure: Level 7 : offset = ( 0, 0, 0) size = ( 128, 128, 128) Level 8 : offset = ( 47, 48, 44) size = ( 68, 60, 82) Level 9 : offset = ( 12, 12, 12) size = ( 88, 72, 116) Level 10 : offset = ( 12, 13, 12) size = ( 126, 92, 184) Level 11 : offset = ( 12, 12, 12) size = ( 202, 136, 318) Level 12 : offset = ( 12, 12, 12) size = ( 354, 224, 586) Level 13 : offset = ( 13, 13, 12) size = ( 658, 396, 1124) ------------------------------------------------------------- - Finest level : extent = 4819.34 x 2900.39 x 8232.42 h-3 kpc**3 mtotgrid = 8.68629e+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... (384, 384, 384) - Performing noise convolution on level 11... - Performing density convolution... (656, 656, 656) - Performing noise convolution on level 12... - Performing density convolution... (1192, 1192, 1192) - Performing noise convolution on level 13... - Performing density convolution... (2264, 2264, 2264) - Density calculation took 778.444s with 32 threads. - Top grid mean density is off by -1.5886e-07, correcting... ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 5.28212e-07 - Poisson solver took 151.393s with 32 threads. - Computing 2LPT term.... - Solving 2LPT Poisson equation ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 6 steps to 1.20764e-07 - Poisson solver took 176.23s with 32 threads. - Performing hybrid Poisson step... (2248, 2248, 2248) - velocity component 0 : sigma = 0.169224 - velocity component 1 : sigma = 0.154493 - velocity component 2 : sigma = 0.370977 ============================================================= 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... (384, 384, 384) - Performing noise convolution on level 11... - Performing density convolution... (656, 656, 656) - Performing noise convolution on level 12... - Performing density convolution... (1192, 1192, 1192) - Performing noise convolution on level 13... - Performing density convolution... (2264, 2264, 2264) - Density calculation took 793.577s with 32 threads. - Top grid mean density is off by -1.53192e-07, correcting... ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 4.18645e-07 - Poisson solver took 137.17s with 32 threads. - Writing baryon potential ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 6 steps to 8.85975e-08 - Poisson solver took 168.336s with 32 threads. - Performing hybrid Poisson step... (2248, 2248, 2248) - velocity component 0 : sigma = 0.16925 - velocity component 1 : sigma = 0.15457 - velocity component 2 : sigma = 0.370917 ============================================================= 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... (384, 384, 384) - Performing noise convolution on level 11... - Performing density convolution... (656, 656, 656) - Performing noise convolution on level 12... - Performing density convolution... (1192, 1192, 1192) - Performing noise convolution on level 13... - Performing density convolution... (2264, 2264, 2264) - Density calculation took 779.173s with 32 threads. - Top grid mean density is off by -1.5886e-07, correcting... ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 7.81634e-07 - Poisson solver took 150.128s with 32 threads. ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 6 steps to 1.92417e-07 - Poisson solver took 177.132s with 32 threads. - Performing hybrid Poisson step... (2248, 2248, 2248) - max. x-displacement of HR particles is -14.191147 [mean dx] - max. y-displacement of HR particles is 7.888302 [mean dx] - max. z-displacement of HR particles is 17.294173 [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... (384, 384, 384) - Performing noise convolution on level 11... - Performing density convolution... (656, 656, 656) - Performing noise convolution on level 12... - Performing density convolution... (1192, 1192, 1192) - Performing noise convolution on level 13... - Performing density convolution... (2264, 2264, 2264) - Density calculation took 791.155s with 32 threads. - Top grid mean density is off by -1.53192e-07, correcting... ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 5 steps to 4.18479e-07 - Poisson solver took 149.125s with 32 threads. ------------------------------------------------------------- - Invoking multi-grid Poisson solver... - Converged in 6 steps to 8.67132e-08 - Poisson solver took 164.398s with 32 threads. - Performing hybrid Poisson step... (2248, 2248, 2248) - Gadget2 : writing 291816596 particles to file... - type 0 : 140876456 [m=4.96123e-07] - type 1 : 140876456 [m=2.46971e-06] - type 2 : 10063684 [m=0] - Gadget2 : set initial gas temperature to 140.40 K/mu ============================================================= - Wrote output file 'ic_agora_m12m_ref13_rad5-chull.ics' using plugin 'gadget2'... - Done! wall time: 9272 sec = 0.11 day = 2:34:31 cpu time: 82.4 hr Tue Jan 19 22:55:32 CST 2016