Question regarding two-phase-clsvof.h

[Andrés Castillo-Castellanos]

Hello (and happy new year)

I have a question about the hybrid approach. We (the user) set d, then the code instantiates phi, and f using 

 

event init (i = 0)
{
  vertex scalar phi[];
  foreach_vertex()
    phi[] = (d[] + d[-1] + d[0,-1] + d[-1,-1])/4.;
  fractions (phi, f);

}

which is fine in 2D. Shouldn't there be a special treatment for 3D? 

Best regards,