processMesh.m

[hbbdyzf]

Hello, gael
there is bug in processMesh.m, ''triangleNew(i).air = (1/2)*norm(cross(A-C,B-C)); '' triangleNew(i).air is always more than 0.

[gBringout]

Hi,

Their is definitively something weird with those lines:

CoilDesign/streamBEM/otherFunctions/processMesh.m

Lines 26 to 27 in 223d8c2

	%calculate the air, according to: http://geomalgorithms.com/a01-_area.html
	triangleNew(i).air = (1/2)*norm(cross(A-C,B-C));

The norm() function juste send positive value. This cannot be a signed results.
Moreover, the idea of the triangle orientation via its normal is defined in the link for a 2D triangles. Here we are working with 3D triangles.

This should be investigated