Intro This is my personal website and research blog. Here you will find what I’m working on, some information about me, and my CV/resume. If you’re interested in my research, contact info, CV/res...

MHD Piecewise Linear Reconstruction This week I started on the MHD piecewise linear reconstruction in the characteristic variables (PLMC). The steps I’ve completed are as follows: Do some rese...

MHD Advecting Field Loop This is done and in PR #273. Last week I discovered that the issues I was having with it were just due to some errors in how I was computing the mean magnetic pressure du...

Angled Linear Waves After a lot more work this week I discovered that my implementation of linear waves from the vector potential did not in fact pass with grid aligned or angled waves. Since this...

Circularly Polarized Alfvên Wave Test (CPAW) I finished the CPAW test. There were a few other bugs in the initial conditions but I got those sorted and cleaned up the code, wrote the tests, and su...

Circularly Polarized Alfvên Wave Test (CPAW) I worked this week to add the CPAW initial conditions and test based on the CPAW test in Gardiner & Stone 2008 and the Athena++ source code. The im...

MHD Debugging I’m pretty sure that MHD is finally debugged! I found a couple of bugs and typos in initial conditions in the shock tubes. There were typos in the papers I was referencing, the fixed...

MHD Debugging I found the issue I was struggling with last week with some help from Alwin Mao. It turns out that since MHD is face centered it needs the ghost cell closest to the grid to always ha...

MHD Debugging It looks like the weird edge effect issues I was seeing might just be oscillations due to the higher order Van Leer method as it’s present in hydro builds as well. It turns out that...

HLLD Refactor I finished refactoring the HLLD solver to use structs instead of scalar variables this week (PR #237). It uses a few less registers now, is much more readable, and several hundred li...