Initialization with 3D ellipsoid results in cylinder

When using an ellipsoid to initialize a 3D cell, a cylinder is created instead:

ellipsoid_bug

This is fixed by changing line 430 in init_cell_objects.cpp (see ellipsoid.patch) :

double p  = (sqr(d.x))/sqr(a) + (sqr(d.y))/sqr(b) + (c>0 && d.z>0 ? (sqr(d.z))/sqr(c): 0.0);

to

double p  = (sqr(d.x))/sqr(a) + (sqr(d.y))/sqr(b) + (c>0 && abs(d.z)>0 ? (sqr(d.z))/sqr(c): 0.0);

Now, the initial cell indeed looks like an ellipsoid:

ellipsoid_fixed