Skip to content
rhaas80 edited this page Jun 7, 2022 · 50 revisions

Getting started

Con2Prim

As a first step we will use Reprimand to reduce the complexity of code we implement. Instructions for this are found here

NERSC Hackathon

CarpetX ideas

  • ideas: call scheduled functions with full MFab: done
  • ideas: use OpenMP tasking or similar to have a more task based scheduling
  • work item ideas

CarpetX implementation tasks

  • vertex centered grid functions by default: done
  • test multi-node GPU support: done
  • checkpoint / restart
  • test a qc0 run or similar to compare results and speed to Carpet w/o subcycling in time
  • subcycling in time, likely using runge-kutta prolongation
    • allow flexibiliy of subcycling vs no subcycling on a per level case, currently partially implemented

Workshop tasks

  • Get the image setup to run code [Steve]
  • Get the shocktube/KHI notebook working [Lorenzo]
  • Slides [Federico]
  • Fix AMR [Erik]

Workshop Timeline

  • 16/05 Create workshop branch (code freeze) with AMR
  • 16/05 Ready the KHI thorn for AsterX
  • 21/05 Final pre-compiled Cactus on K8
  • 21/05 Ready visualization of plots
  • 24/05 Copying pre-compiled Cactus for each user
  • 24/05 Presentation ready (common slides)
  • 26/05 Ready tutorial Jupyter notebook to be shared with participants
  • 07/06 Mock run with ~12

Workshop tutorial schedule

The authorative source of information is the program page.

title date start time end time needs GPU
New user tutorial June 13 2:05PM PDT 3:50PM PDT no
Thorn writing June 13 3:20PM PDT 6:05PM PDT no
Kuibit June 16 4:20PM PDT 6:05PM PDT no
GPus June 17 4:20PM PDT 6:05PM PDT yes

CarpetX desired feature list

Science runs would be say a MHD driven supernova with tracer particles and a neutrino evolution scheme. This would require these abilities:

  • cell centered grid with AMR for fluid
  • vertex centered grid with AMR for metric
  • staggered grid for magnetic field
  • subcycling in time
    • Runge-Kutta prolongation into mesh refinement boundaries to avoid "buffer" zones (which AMReX may not support)
  • some sort refluxing for fluid due to subcycling in time
  • some sort of monopole avoiding prolongation scheme for B
  • some sort of scheme that conserves magnetic flux at mesh refinement boundaries
  • interpolators to get values of fields at particle locations
  • reduction operators for min, max, sum etc.
  • checkpoint and recovery
  • 3d output for visualization: done
  • GPU support
    • for RHS calculation: done
    • for AMR operators
    • for everything else
Clone this wiki locally