[cig-commits] commit: Bug fix.

[Mercurial]

changeset:   97:d717a6f624dc
user:        JulianGiordani
date:        Wed May 28 02:09:22 2008 +0000
files:       MaterialPoints/src/SwarmAdvectorD.c
description:
Bug fix.
Was always assuming the velocity vector was of size 3, even in 2D. But that's not the case anymore.
Fixed by a conditional statement.


diff -r 8ceb22bb0e67 -r d717a6f624dc MaterialPoints/src/SwarmAdvectorD.c
--- a/MaterialPoints/src/SwarmAdvectorD.c	Wed May 14 05:44:12 2008 +0000
+++ b/MaterialPoints/src/SwarmAdvectorD.c	Wed May 28 02:09:22 2008 +0000
@@ -207,6 +207,7 @@ Bool _SwarmAdvectorD_TimeDeriv( void* sw
 	coord = Variable_GetPtrDouble( self->variable, array_I );
 
 	result = FieldVariable_InterpolateValueAt( velocityField, coord, timeDeriv );
+	if( velocityField->dim == 3 )
 	        timeDeriv[ K_AXIS ] = 0.0;