[Bf-blender-cvs] SVN commit: /data/svn/bf-blender [16443] trunk/blender/source/blender/src/ meshlaplacian.c:

[Brecht Van Lommel]

Revision: 16443
          http://projects.blender.org/plugins/scmsvn/viewcvs.php?view=rev&root=bf-blender&revision=16443
Author:   blendix
Date:     2008-09-09 22:00:57 +0200 (Tue, 09 Sep 2008)

Log Message:
-----------

Fix for bug #17443: make bone heat weighting solve a bit less
sensitive to poorly shaped geometry, should succeed solving
in more cases now.

Modified Paths:
--------------
    trunk/blender/source/blender/src/meshlaplacian.c

Modified: trunk/blender/source/blender/src/meshlaplacian.c
===================================================================
--- trunk/blender/source/blender/src/meshlaplacian.c	2008-09-09 19:00:16 UTC (rev 16442)
+++ trunk/blender/source/blender/src/meshlaplacian.c	2008-09-09 20:00:57 UTC (rev 16443)
@@ -178,6 +178,7 @@
 		varea[i1] += (obtuse == 1)? area: area*0.5;
 		varea[i2] += (obtuse == 2)? area: area*0.5;
 		varea[i3] += (obtuse == 3)? area: area*0.5;
+		//printf("area %f\n", area);
 	}
 	else {
 		len1= VecLenf(v2, v3);
@@ -191,7 +192,10 @@
 		varea[i1] += (t2 + t3)*0.25f;
 		varea[i2] += (t1 + t3)*0.25f;
 		varea[i3] += (t1 + t2)*0.25f;
+		//printf("varea %f %f %f\n", t1, t2, t3);
 	}
+
+	//printf("triangle area %f %f %f\n", t1, t2, t3);
 }
 
 static void laplacian_triangle_weights(LaplacianSystem *sys, int f, int i1, int i2, int i3)
@@ -204,7 +208,7 @@
 	v3= sys->verts[i3];
 
 	/* instead of *0.5 we divided by the number of faces of the edge, it still
-	   needs to be varified that this is indeed the correct thing to do! */
+	   needs to be verified that this is indeed the correct thing to do! */
 	t1= cotan_weight(v1, v2, v3)/laplacian_edge_count(sys->edgehash, i2, i3);
 	t2= cotan_weight(v2, v3, v1)/laplacian_edge_count(sys->edgehash, i3, i1);
 	t3= cotan_weight(v3, v1, v2)/laplacian_edge_count(sys->edgehash, i1, i2);
@@ -229,7 +233,7 @@
 	}
 }
 
-LaplacianSystem *laplacian_system_construct_begin(int totvert, int totface)
+LaplacianSystem *laplacian_system_construct_begin(int totvert, int totface, int lsq)
 {
 	LaplacianSystem *sys;
 
@@ -248,6 +252,8 @@
 	/* create opennl context */
 	nlNewContext();
 	nlSolverParameteri(NL_NB_VARIABLES, totvert);
+	if(lsq)
+		nlSolverParameteri(NL_LEAST_SQUARES, NL_TRUE);
 
 	sys->context= nlGetCurrent();
 
@@ -292,7 +298,7 @@
 	for(a=0; a<totvert; a++) {
 		if(sys->areaweights) {
 			if(sys->varea[a] != 0.0f)
-				sys->varea[a]= 0.5f/sys->varea[a];
+				sys->varea[a]= 0.5f/sys->varea[a]; //MAX2(sys->varea[a], 0.001f);
 		}
 		else
 			sys->varea[a]= 1.0f;
@@ -631,7 +637,7 @@
 	}
 
 	/* create laplacian */
-	sys = laplacian_system_construct_begin(me->totvert, totface);
+	sys = laplacian_system_construct_begin(me->totvert, totface, 1);
 
 	sys->heat.mesh= me;
 	sys->heat.verts= verts;
@@ -933,7 +939,7 @@
 	}
 
 	/* create laplacian */
-	sys = laplacian_system_construct_begin(totvert, totface);
+	sys = laplacian_system_construct_begin(totvert, totface, 0);
 
 	sys->rigid.mesh= em;
 	sys->rigid.R = MEM_callocN(sizeof(float)*3*3*totvert, "RigidDeformR");