Full-grid C++ doesn't match Full-grid matlab.
For Vlasov43
./fk6d -p 4 -l 3 -d 3 -n 369 -w 10 -f -c 0.1
matlab gives ...
but C++ gives ...
Notes :
- Looks fine for Vlasov8, so that means the problem is likely in the Poisson solve.
- Vlasov4 (
-p 0
) C++ goes unstable after about 900 time steps. I'm running the matlab out to that now to check.
I may get to tracking this down today, but if not we can take a look tomorrow.