vecsetvec.c
Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024 static char const Id[] = "$Id: vecsetvec.c 2681 2009-07-31 04:30:53Z scholz $\n\n";
00025 static char const Td[] = "$Today$\n\n";
00026
00027 #include "util.h"
00028
00029 int VecSetVec(Vec vec, PetscReal *vec3, int dim)
00030 {
00031 MagparFunctionInfoBegin;
00032
00033 int vsize;
00034 ierr = VecGetLocalSize(vec,&vsize);CHKERRQ(ierr);
00035 if (vsize%dim!=0) {
00036 SETERRQ3(PETSC_ERR_ARG_CORRUPT,
00037 "vsize=%i, dim=%i, %i!=0!\n",
00038 vsize,dim,vsize%dim
00039 );
00040 }
00041
00042 PetscReal *ta_vec;
00043 ierr = VecGetArray(vec,&ta_vec);CHKERRQ(ierr);
00044 for (int j=0; j<vsize/dim; j++) {
00045 for (int i=0; i<dim; i++) {
00046 ta_vec[dim*j+i]=vec3[i];
00047 }
00048 }
00049 ierr = VecRestoreArray(vec,&ta_vec);CHKERRQ(ierr);
00050
00051 MagparFunctionInfoReturn(0);
00052 }
00053