Bug in wedgePolyPatchC

su_junwei · February 12, 2009, 07:09

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

henry · February 12, 2009, 07:13

> 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

February 12, 2009, 07:09	Dear Foamers? There seems	#1
su_junwei Senior Member su junwei Join Date: Mar 2009 Location: Xi'an China Posts: 151 Rep Power: 20	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

February 12, 2009, 07:13	> If the face in this patch is	#2
henry Senior Member Join Date: Mar 2009 Posts: 854 Rep Power: 22	> 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