Dear Foamers? There seems
 

Dear Foamers?

There seems to be a bug in the implementation of initTransforms() in wedgePolyPatch.C(line 44)

If the face[0] in this patch is in xy plane, patchNormal_ will equal {0 0 1} (after normal treatment) .
centerNormal_ is {0 0 0.5} before normal treatment and {0 0 1} after normal treatment.

While evaluating the axis_ using
axis_ = centreNormal_ ^ patchNormal_;
we get a zero vector for axis_.

In the following normal treatment for axis_, we will get a zero division exception.

It seems that the method for evaluation of centreNormal_ can't always give us a feasible result.

Is it right?

Junwei su

> If the face in this patch is
 

> If the face[0] in this patch is in xy plane

This should not be the case. Take a look at the documentation on how to setup wedge cases.

H