[petsc-users] DMPlexCreateLocalVec Within a function

[Matthew Kury]

Hello, 

	I have done some reading and I have not been able to find out what the proper way to create and use vectors within a function.

	If I have a function that passes in dm and V then: 
+++
		DMPlexGetLocalVec(dm, &VLocal) ;
		DMGlobalToLocalBegin(dm,V,INSERT_VALUES,VLocal); 
		DMGlobalToLocalEnd(dm,V,INSERT_VALUES,VLocal) ;

		VecGetArray(VLoc,VLocalArray);

		“Stuff with VLocalArray”

		VecRestoreArray(VLoc,VLocalArray);
		
		DMLocalToGlobalBegin(dm,VLocal,ADD_VALUES,V);
		DMLocalToGlobalEnd(dm,VLocal,ADD_VALUES,V);

		DMRestoreLocalVector(dm,&VLocal);
+++

	The second time I call the function I get errors with the VecGetArray() function.  I would appreciate any advice as to what is going wrong. 

Thank you, 

Matthew