Re: [petsc-users] *****SPAM*****Create vector with evenly spaced values

[Florian Lindner]

Am 27.08.2014 14:12, schrieb Florian Lindner:
Hello,
is there a way to create a vector with evenly spaced values in an
interval, like the numpy.arange function?
http://docs.scipy.org/doc/numpy/reference/generated/numpy.arange.html
or the linspace function?

Ok, I gave it as try. Do I need the VecAssembly the vector?


void arange(size_t start, size_t stop)
    {
      // given a vector of appropriate size
      PetscScalar ierr;
      PetscScalar *a;
      ierr = PetscMalloc(sizeof(PetscScalar) * (stop-start), &a); 
CHKERRV(ierr);
      size_t pos = 0;
      for (size_t i = start; i < stop; i++) {
        a[pos] = i;
        pos++;
      }
      cout << "Array filled." << endl;
      VecReplaceArray(vector, a);
      // VecAssemblyBegin(vector); VecAssemblyEnd(vector);  Do I need 
that?
    }

Any comments?

Thanks,
Florian