David Bold

On 12/20/19 2:09 PM, Tom Nicholas wrote: > at that point it becomes easier to just use 'boutcore'. > > But you would have to import BOUT's c++ code for...

I don't know how to use `apply_ufunc` How do I ensure that I get always full 3D data sets? Is there a way to get several fields? Generally it works:...

I think a first step into nicer integration, would be adding a wrapper of some sort, that lets you easily get a Field3D from the dataset. What would be a...

How would that imply for the advect example? would that return advect(n,phi,B)? I think that is a bit confusing, and I like to keep it close to the C++ interface,...

> Once for each set of data files, or once ever? Once per run. MPI needs to be initialised exactly once, and only finalised at the end (I think this...

There are also two different expressions for `s_i` used: ```c++ BoutReal s_i = (nesheath > 1e-5) ? nisheath / nesheath : 0.0; // Concentration BoutReal s_i = clip(nisheath / floor(nesheath,...

I have written a small MMS test and have been running it for the div_a_grad_perp_nonorthog, but it does not seem to converge:  For reference, the Delp2 operator:  I...

Can you either change this PR to go into master or open a new one that goes into master: https://github.com/bendudson/hermes-3/compare/bendudson:master...bendudson:nonorthogonal-divops-mms-test-minimal?expand=1 Like this it is hard to see the changes you...

Thanks John, that also explains why the matrix is re-created on every iteration :-)

I will see whether `petsc3damg` will work for us, otherwise I will try to fix `petsc`. If I do not get it fixed, deleting sounds like the right thing to...