[cig-commits] commit: removing the PESKy wrapper crap and replacing with PETSc native calls

[Mercurial]

changeset:   123:d7db09e1378c
parent:      119:480cfbde2b53
user:        dave at tango.vpac.org
date:        Thu Oct 16 17:28:18 2008 +1100
files:       MaterialPoints/src/ParticleFeVariable.c
description:
removing the PESKy wrapper crap and replacing with PETSc native calls


diff -r 480cfbde2b53 -r d7db09e1378c MaterialPoints/src/ParticleFeVariable.c
--- a/MaterialPoints/src/ParticleFeVariable.c	Fri Oct 03 08:05:25 2008 +1000
+++ b/MaterialPoints/src/ParticleFeVariable.c	Thu Oct 16 17:28:18 2008 +1100
@@ -296,13 +296,16 @@ void ParticleFeVariable_Update( void* ma
 	printf( "***\n*** Updating %s\n***\n", self->name );
 
 	/* Initialise Vectors */
-	Vector_Zero( self->assemblyVector->vector );
-	Vector_Zero( self->massMatrix->vector );
+	//Vector_Zero( self->assemblyVector->vector );
+	//Vector_Zero( self->massMatrix->vector );
+	VecSet( self->assemblyVector->vector, 0.0 );
+	VecSet( self->massMatrix->vector, 0.0 );
 
 	ForceVector_Assemble( self->assemblyVector );
 	ForceVector_Assemble( self->massMatrix );
 
-	Vector_PointwiseDivide( self->assemblyVector->vector, self->assemblyVector->vector, self->massMatrix->vector );
+	//Vector_PointwiseDivide( self->assemblyVector->vector, self->assemblyVector->vector, self->massMatrix->vector );
+	VecPointwiseDivide( self->assemblyVector->vector, self->assemblyVector->vector, self->massMatrix->vector );
 
 	SolutionVector_UpdateSolutionOntoNodes( self->assemblyVector );
 }