Posts in category parameter

Meeting Recap 12/17/2012

  • Spent some time working on Christina's setup and improving refinement criteria. 'I' shape due to density gradients produced due to bending magnetic field lines and strong shear. In 3D this I is a cylinder. So I modified module to use three refinement objects.
    • First maintain refinement to level 3 within a 2 pc wide region around the initial interface (nested inside of coarser patches)
    • Next allow for density gradient based refinement to level 3 within a 10 pc region centered on the initial interface (nested inside of coarser patches)
    • And finally allow for jeans length based refinement (64 cells per jeans length) all the way to level 4 so that collapsing structures will be better resolved.

  • 3D MHD with self-gravity on bluestreak runs at 1475 cell updates/second/core
  • On 32 bluestreak cores, at 323 + 3 it was spending 76% of the time advancing, 23% solving elliptic equations, and < 1% everything else… and was slotted to finish in 8 days (30 Myr).
  • In principle if everything scales up - the initial resolution of 1283+4 could require between 256 and 4096 times as many resources depending on the filling fraction of the highest level…
  • At 256 this would give 16 days on ¼ of bluestreak (4000) which seems doable? This is 1.5 Mhrs of CPU time on bluestreak - or probably .4 Mhrs on kraken
  • Things to work on over break…
    • MHD convergence - fast waves now converge - slow waves still do not.
    • PPM convergence
    • FLD solver
    • Documentation ???
    • Developer's manual ???
  • Bash script for running parameter studies…
    interporders=( 1 2 3 )
    Gmxs=( 100 200 400 )
    Waves=( 1 2 3 4 5 6 7 )
    mkdir out
    mkdir testdata
    for (( i=0;i<${#interporders[@]};i++)); do                                                                                                            
        interporder=${interporders[i]}
        ../../scripts/subst.s InterpOrder  $interporder solver.data
        for (( j=0;j<${#Waves[@]};j++)); do                                                                                                               
            wave=${Waves[j]}
            ../../scripts/subst.s WaveFamily $wave problem.data
            for (( m=0;m<${#Gmxs[@]};m++)); do                                                                                                            
                Gmx=${Gmxs[m]}
                ../../scripts/subst.s Gmx $Gmx,1,1 global.data
                mpirun -n 1 ../../astrobear > output.out
                grep dx output.out | awk '{print $3, $4}' >> scaling_${interporder}_${j}.curve
                mv out/data.curve data_${interporder}_${j}_${Gmx}.curve
                cat output.out
            done
        done
    done