From b0162a19fd9d77c7facee1f64d45edf1363287cc Mon Sep 17 00:00:00 2001
From: Ton Roosendaal <ton@blender.org>
Date: Sun, 19 Dec 2004 14:19:20 +0000
Subject: [PATCH] Bug fix #2026

Slightly altered rules for calculating vertexnormals. By only averaging
face normals from faces actually set 'smooth', the result looks much more
nice (for example on a cylinder with caps solid). Vertex normals not being
used by smooth faces are set to the face normal direction.

Shows both in editor as rendering.
---
 .../intern/convertBlenderScene.c              | 141 ++++++++++--------
 source/blender/src/editmesh_lib.c             | 114 ++++++++------
 2 files changed, 145 insertions(+), 110 deletions(-)

diff --git a/source/blender/renderconverter/intern/convertBlenderScene.c b/source/blender/renderconverter/intern/convertBlenderScene.c
index 927f1b8bca1..f43dca9ecaa 100644
--- a/source/blender/renderconverter/intern/convertBlenderScene.c
+++ b/source/blender/renderconverter/intern/convertBlenderScene.c
@@ -129,7 +129,6 @@ static Material *give_render_material(Object *ob, int nr);
 static void split_u_renderfaces(int startvlak, int startvert, int usize, int plek, int cyclu);
 static void split_v_renderfaces(int startvlak, int startvert, int usize, int vsize, int plek, int cyclu, int cyclv);
 static int contrpuntnormr(float *n, float *puno);
-static void normalenrender(int startvert, int startvlak);
 static void as_addvert(VertRen *v1, VlakRen *vlr);
 static void as_freevert(VertRen *ver);
 static void autosmooth(int startvert, int startvlak, int degr);
@@ -453,7 +452,7 @@ static int contrpuntnormr(float *n, float *puno)
 
 /* ------------------------------------------------------------------------- */
 
-static void normalenrender(int startvert, int startvlak)
+static void calc_vertexnormals(int startvert, int startvlak)
 {
 	VlakRen *vlr;
 	VertRen *ver, *adrve1, *adrve2, *adrve3, *adrve4;
@@ -469,34 +468,35 @@ static void normalenrender(int startvert, int startvlak)
 	/* calculate cos of angles and point-masses */
 	for(a= startvlak; a<R.totvlak; a++) {
 		vlr= RE_findOrAddVlak(a);
+		if(vlr->flag & ME_SMOOTH) {
+			adrve1= vlr->v1;
+			adrve2= vlr->v2;
+			adrve3= vlr->v3;
+			adrve4= vlr->v4;
 
-		adrve1= vlr->v1;
-		adrve2= vlr->v2;
-		adrve3= vlr->v3;
-		adrve4= vlr->v4;
+			VecSubf(n1, adrve2->co, adrve1->co);
+			Normalise(n1);
+			VecSubf(n2, adrve3->co, adrve2->co);
+			Normalise(n2);
+			if(adrve4==0) {
+				VecSubf(n3, adrve1->co, adrve3->co);
+				Normalise(n3);
 
-		VecSubf(n1, adrve2->co, adrve1->co);
-		Normalise(n1);
-		VecSubf(n2, adrve3->co, adrve2->co);
-		Normalise(n2);
-		if(adrve4==0) {
-			VecSubf(n3, adrve1->co, adrve3->co);
-			Normalise(n3);
+				*(tfl++)= saacos(-n1[0]*n3[0]-n1[1]*n3[1]-n1[2]*n3[2]);
+				*(tfl++)= saacos(-n1[0]*n2[0]-n1[1]*n2[1]-n1[2]*n2[2]);
+				*(tfl++)= saacos(-n2[0]*n3[0]-n2[1]*n3[1]-n2[2]*n3[2]);
+			}
+			else {
+				VecSubf(n3, adrve4->co, adrve3->co);
+				Normalise(n3);
+				VecSubf(n4, adrve1->co, adrve4->co);
+				Normalise(n4);
 
-			*(tfl++)= saacos(-n1[0]*n3[0]-n1[1]*n3[1]-n1[2]*n3[2]);
-			*(tfl++)= saacos(-n1[0]*n2[0]-n1[1]*n2[1]-n1[2]*n2[2]);
-			*(tfl++)= saacos(-n2[0]*n3[0]-n2[1]*n3[1]-n2[2]*n3[2]);
-		}
-		else {
-			VecSubf(n3, adrve4->co, adrve3->co);
-			Normalise(n3);
-			VecSubf(n4, adrve1->co, adrve4->co);
-			Normalise(n4);
-
-			*(tfl++)= saacos(-n4[0]*n1[0]-n4[1]*n1[1]-n4[2]*n1[2]);
-			*(tfl++)= saacos(-n1[0]*n2[0]-n1[1]*n2[1]-n1[2]*n2[2]);
-			*(tfl++)= saacos(-n2[0]*n3[0]-n2[1]*n3[1]-n2[2]*n3[2]);
-			*(tfl++)= saacos(-n3[0]*n4[0]-n3[1]*n4[1]-n3[2]*n4[2]);
+				*(tfl++)= saacos(-n4[0]*n1[0]-n4[1]*n1[1]-n4[2]*n1[2]);
+				*(tfl++)= saacos(-n1[0]*n2[0]-n1[1]*n2[1]-n1[2]*n2[2]);
+				*(tfl++)= saacos(-n2[0]*n3[0]-n2[1]*n3[1]-n2[2]*n3[2]);
+				*(tfl++)= saacos(-n3[0]*n4[0]-n3[1]*n4[1]-n3[2]*n4[2]);
+			}
 		}
 	}
 
@@ -511,47 +511,66 @@ static void normalenrender(int startvert, int startvlak)
 	tfl= adrco;
 	for(a=startvlak; a<R.totvlak; a++) {
 		vlr= RE_findOrAddVlak(a);
+		
+		if(vlr->flag & ME_SMOOTH) {
 
-		adrve1= vlr->v1;
-		adrve2= vlr->v2;
-		adrve3= vlr->v3;
-		adrve4= vlr->v4;
+			adrve1= vlr->v1;
+			adrve2= vlr->v2;
+			adrve3= vlr->v3;
+			adrve4= vlr->v4;
 
-		temp= adrve1->n;
-		fac= *(tfl++);
-		if( vlr->flag & R_NOPUNOFLIP);
-		else if( contrpuntnormr(vlr->n, temp) ) fac= -fac ;
-		*(temp++) +=fac*vlr->n[0];
-		*(temp++) +=fac*vlr->n[1];
-		*(temp)   +=fac*vlr->n[2];
-
-		temp= adrve2->n;
-		fac= *(tfl++);
-		if( vlr->flag & R_NOPUNOFLIP);
-		else  if( contrpuntnormr(vlr->n, temp) ) fac= -fac ;
-		*(temp++) +=fac*vlr->n[0];
-		*(temp++) +=fac*vlr->n[1];
-		*(temp)   +=fac*vlr->n[2];
-
-		temp= adrve3->n;
-		fac= *(tfl++);
-		if( vlr->flag & R_NOPUNOFLIP);
-		else if( contrpuntnormr(vlr->n, temp) ) fac= -fac ;
-		*(temp++) +=fac*vlr->n[0];
-		*(temp++) +=fac*vlr->n[1];
-		*(temp)   +=fac*vlr->n[2];
-
-		if(adrve4) {
-			temp= adrve4->n;
+			temp= adrve1->n;
 			fac= *(tfl++);
 			if( vlr->flag & R_NOPUNOFLIP);
 			else if( contrpuntnormr(vlr->n, temp) ) fac= -fac ;
 			*(temp++) +=fac*vlr->n[0];
 			*(temp++) +=fac*vlr->n[1];
 			*(temp)   +=fac*vlr->n[2];
+
+			temp= adrve2->n;
+			fac= *(tfl++);
+			if( vlr->flag & R_NOPUNOFLIP);
+			else  if( contrpuntnormr(vlr->n, temp) ) fac= -fac ;
+			*(temp++) +=fac*vlr->n[0];
+			*(temp++) +=fac*vlr->n[1];
+			*(temp)   +=fac*vlr->n[2];
+
+			temp= adrve3->n;
+			fac= *(tfl++);
+			if( vlr->flag & R_NOPUNOFLIP);
+			else if( contrpuntnormr(vlr->n, temp) ) fac= -fac ;
+			*(temp++) +=fac*vlr->n[0];
+			*(temp++) +=fac*vlr->n[1];
+			*(temp)   +=fac*vlr->n[2];
+
+			if(adrve4) {
+				temp= adrve4->n;
+				fac= *(tfl++);
+				if( vlr->flag & R_NOPUNOFLIP);
+				else if( contrpuntnormr(vlr->n, temp) ) fac= -fac ;
+				*(temp++) +=fac*vlr->n[0];
+				*(temp++) +=fac*vlr->n[1];
+				*(temp)   +=fac*vlr->n[2];
+			}
 		}
 	}
-
+	/* do solid faces */
+	for(a=startvlak; a<R.totvlak; a++) {
+		vlr= RE_findOrAddVlak(a);
+		if((vlr->flag & ME_SMOOTH)==0) {
+			float *f1= vlr->v1->n;
+			if(f1[0]==0.0 && f1[1]==0.0 && f1[2]==0.0) VECCOPY(f1, vlr->n);
+			f1= vlr->v2->n;
+			if(f1[0]==0.0 && f1[1]==0.0 && f1[2]==0.0) VECCOPY(f1, vlr->n);
+			f1= vlr->v3->n;
+			if(f1[0]==0.0 && f1[1]==0.0 && f1[2]==0.0) VECCOPY(f1, vlr->n);
+			if(vlr->v4) {
+				f1= vlr->v4->n;
+				if(f1[0]==0.0 && f1[1]==0.0 && f1[2]==0.0) VECCOPY(f1, vlr->n);
+			}			
+		}
+	}
+	
 	/* normalise vertex normals */
 	for(a=startvert; a<R.totvert; a++) {
 		ver= RE_findOrAddVert(a);
@@ -1548,7 +1567,7 @@ static void init_render_mesh(Object *ob)
 	}
 	
 	if(do_puno) {
-		normalenrender(totverto, totvlako);
+		calc_vertexnormals(totverto, totvlako);
 		do_puno= 0;
 	}
 	
@@ -1560,7 +1579,7 @@ static void init_render_mesh(Object *ob)
 		do_puno= 1;
 	}
 	
-	if(do_puno) normalenrender(totverto, totvlako);
+	if(do_puno) calc_vertexnormals(totverto, totvlako);
 	
 	mesh_modifier(ob, 'e');  // end
 }
@@ -3151,7 +3170,7 @@ static void do_displacement(Object *ob, int startface, int numface, int startver
 	}
 	
 	/* Recalc vertex normals */
-	normalenrender(startvert, startface);
+	calc_vertexnormals(startvert, startface);
 }
 
 static void displace_render_face(VlakRen *vlr, float *scale)
diff --git a/source/blender/src/editmesh_lib.c b/source/blender/src/editmesh_lib.c
index 452c1585a13..f1fda764d9b 100644
--- a/source/blender/src/editmesh_lib.c
+++ b/source/blender/src/editmesh_lib.c
@@ -1284,61 +1284,77 @@ void vertexnormals(int testflip)
 	/* calculate cosine angles and add to vertex normal */
 	efa= em->faces.first;
 	while(efa) {
-		VecSubf(n1, efa->v2->co, efa->v1->co);
-		VecSubf(n2, efa->v3->co, efa->v2->co);
-		Normalise(n1);
-		Normalise(n2);
+		if(efa->flag & ME_SMOOTH) {
+			VecSubf(n1, efa->v2->co, efa->v1->co);
+			VecSubf(n2, efa->v3->co, efa->v2->co);
+			Normalise(n1);
+			Normalise(n2);
 
-		if(efa->v4==0) {
-			VecSubf(n3, efa->v1->co, efa->v3->co);
-			Normalise(n3);
+			if(efa->v4==0) {
+				VecSubf(n3, efa->v1->co, efa->v3->co);
+				Normalise(n3);
+				
+				co[0]= saacos(-n3[0]*n1[0]-n3[1]*n1[1]-n3[2]*n1[2]);
+				co[1]= saacos(-n1[0]*n2[0]-n1[1]*n2[1]-n1[2]*n2[2]);
+				co[2]= saacos(-n2[0]*n3[0]-n2[1]*n3[1]-n2[2]*n3[2]);
+				
+			}
+			else {
+				VecSubf(n3, efa->v4->co, efa->v3->co);
+				VecSubf(n4, efa->v1->co, efa->v4->co);
+				Normalise(n3);
+				Normalise(n4);
+				
+				co[0]= saacos(-n4[0]*n1[0]-n4[1]*n1[1]-n4[2]*n1[2]);
+				co[1]= saacos(-n1[0]*n2[0]-n1[1]*n2[1]-n1[2]*n2[2]);
+				co[2]= saacos(-n2[0]*n3[0]-n2[1]*n3[1]-n2[2]*n3[2]);
+				co[3]= saacos(-n3[0]*n4[0]-n3[1]*n4[1]-n3[2]*n4[2]);
+			}
 			
-			co[0]= saacos(-n3[0]*n1[0]-n3[1]*n1[1]-n3[2]*n1[2]);
-			co[1]= saacos(-n1[0]*n2[0]-n1[1]*n2[1]-n1[2]*n2[2]);
-			co[2]= saacos(-n2[0]*n3[0]-n2[1]*n3[1]-n2[2]*n3[2]);
+			temp= efa->v1->no;
+			if(testflip && contrpuntnorm(efa->n, temp) ) co[0]= -co[0];
+			temp[0]+= co[0]*efa->n[0];
+			temp[1]+= co[0]*efa->n[1];
+			temp[2]+= co[0]*efa->n[2];
 			
-		}
-		else {
-			VecSubf(n3, efa->v4->co, efa->v3->co);
-			VecSubf(n4, efa->v1->co, efa->v4->co);
-			Normalise(n3);
-			Normalise(n4);
+			temp= efa->v2->no;
+			if(testflip && contrpuntnorm(efa->n, temp) ) co[1]= -co[1];
+			temp[0]+= co[1]*efa->n[0];
+			temp[1]+= co[1]*efa->n[1];
+			temp[2]+= co[1]*efa->n[2];
 			
-			co[0]= saacos(-n4[0]*n1[0]-n4[1]*n1[1]-n4[2]*n1[2]);
-			co[1]= saacos(-n1[0]*n2[0]-n1[1]*n2[1]-n1[2]*n2[2]);
-			co[2]= saacos(-n2[0]*n3[0]-n2[1]*n3[1]-n2[2]*n3[2]);
-			co[3]= saacos(-n3[0]*n4[0]-n3[1]*n4[1]-n3[2]*n4[2]);
-		}
-		
-		temp= efa->v1->no;
-		if(testflip && contrpuntnorm(efa->n, temp) ) co[0]= -co[0];
-		temp[0]+= co[0]*efa->n[0];
-		temp[1]+= co[0]*efa->n[1];
-		temp[2]+= co[0]*efa->n[2];
-		
-		temp= efa->v2->no;
-		if(testflip && contrpuntnorm(efa->n, temp) ) co[1]= -co[1];
-		temp[0]+= co[1]*efa->n[0];
-		temp[1]+= co[1]*efa->n[1];
-		temp[2]+= co[1]*efa->n[2];
-		
-		temp= efa->v3->no;
-		if(testflip && contrpuntnorm(efa->n, temp) ) co[2]= -co[2];
-		temp[0]+= co[2]*efa->n[0];
-		temp[1]+= co[2]*efa->n[1];
-		temp[2]+= co[2]*efa->n[2];
-		
-		if(efa->v4) {
-			temp= efa->v4->no;
-			if(testflip && contrpuntnorm(efa->n, temp) ) co[3]= -co[3];
-			temp[0]+= co[3]*efa->n[0];
-			temp[1]+= co[3]*efa->n[1];
-			temp[2]+= co[3]*efa->n[2];
-		}
-		
+			temp= efa->v3->no;
+			if(testflip && contrpuntnorm(efa->n, temp) ) co[2]= -co[2];
+			temp[0]+= co[2]*efa->n[0];
+			temp[1]+= co[2]*efa->n[1];
+			temp[2]+= co[2]*efa->n[2];
+			
+			if(efa->v4) {
+				temp= efa->v4->no;
+				if(testflip && contrpuntnorm(efa->n, temp) ) co[3]= -co[3];
+				temp[0]+= co[3]*efa->n[0];
+				temp[1]+= co[3]*efa->n[1];
+				temp[2]+= co[3]*efa->n[2];
+			}
+		}		
 		efa= efa->next;
 	}
-
+	/* check solid faces */
+	for(efa= em->faces.first; efa; efa= efa->next) {
+		if((efa->flag & ME_SMOOTH)==0) {
+			f1= efa->v1->no;
+			if(f1[0]==0.0 && f1[1]==0.0 && f1[2]==0.0) VECCOPY(f1, efa->n);
+			f1= efa->v2->no;
+			if(f1[0]==0.0 && f1[1]==0.0 && f1[2]==0.0) VECCOPY(f1, efa->n);
+			f1= efa->v3->no;
+			if(f1[0]==0.0 && f1[1]==0.0 && f1[2]==0.0) VECCOPY(f1, efa->n);
+			if(efa->v4) {
+				f1= efa->v4->no;
+				if(f1[0]==0.0 && f1[1]==0.0 && f1[2]==0.0) VECCOPY(f1, efa->n);
+			}
+		}
+	}
+	
 	/* normalise vertex normals */
 	eve= em->verts.first;
 	while(eve) {