Ok, sorry, that
gendir command doesn’t exist any more. It’s called
dirgen now. I’ll edit my previous post to correct that.
Otherwise, the thing to focus on in your output is the b=1000 results, specifically the nearest-neighbour angles: you have an average of 14°, with a min & max of 4° & 29°. That’s pretty bad, ideally the min & max would be fairly tight around the mean. This indicates badly distributed directions – as also shown by the condition numbers for the SH fit, these should ideally be close to 1.
Often, poorly distributed directions arise because the people designing the direction set use a unipolar electrostatic repulsion model, rather than the correct bipolar model (which assumes antipodal symmetry). Looking at the results from the unipolar model though, the results are just as bad, so it’s unlikely to have been optimised for a unipolar model either. Which makes it very likely to have been designed for overplus (where the DW gradients use at least two gradient axes concurrently to boost the effective gradient strength by ✓2).
Not a lot you can do about this after the fact, unfortunately… If eddy won’t take it, then you won’t be able to correct for motion and susceptibility distortions that way. Maybe using TORTOISE (as you’re already doing, judging from your other posts) is your best bet…