[cig-commits] commit:

[Mercurial]

changeset:   57:65081410b863
user:        DavidLee
date:        Fri Mar 14 04:09:43 2008 +0000
files:       Swarm/src/GaussParticleLayout.c
description:
fixed up a bug in the 5 point gaussian quadrature. is working properly now


diff -r 968c59781919 -r 65081410b863 Swarm/src/GaussParticleLayout.c
--- a/Swarm/src/GaussParticleLayout.c	Fri Mar 14 03:01:56 2008 +0000
+++ b/Swarm/src/GaussParticleLayout.c	Fri Mar 14 04:09:43 2008 +0000
@@ -354,7 +354,7 @@ void GaussParticleLayout_GetAbscissaAndW
 			abscissa[1]  = - sqrt( 245.0 - 14.0 * sqrt( 70.0 ) )/21.0;
 			abscissa[2]  = 0.0;
 			abscissa[3]  = - abscissa[1];
-			abscissa[4]  = - abscissa[2];
+			abscissa[4]  = - abscissa[0];
 
 			weight[0] = ( 322.0 - 13.0 * sqrt( 70.0 ) )/900.0;
 			weight[1] = ( 322.0 + 13.0 * sqrt( 70.0 ) )/900.0;