distortvec.c

Go to the documentation of this file.
00001 /*
00002     This file is part of magpar.
00003 
00004     Copyright (C) 2002-2010 Werner Scholz
00005 
00006     www:   http://www.magpar.net/
00007     email: magpar(at)magpar.net
00008 
00009     magpar is free software; you can redistribute it and/or modify
00010     it under the terms of the GNU General Public License as published by
00011     the Free Software Foundation; either version 2 of the License, or
00012     (at your option) any later version.
00013 
00014     magpar is distributed in the hope that it will be useful,
00015     but WITHOUT ANY WARRANTY; without even the implied warranty of
00016     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
00017     GNU General Public License for more details.
00018 
00019     You should have received a copy of the GNU General Public License
00020     along with magpar; if not, write to the Free Software
00021     Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
00022 */
00023 
00024 static char const Id[] = "$Id: distortvec.c 2979 2010-02-22 20:14:26Z scholz $\n\n";
00025 static char const Td[] = "$Today: " __FILE__ "  " __DATE__ "  "  __TIME__  " $\n\n";
00026 
00027 #include "util.h"
00028 
00029 int DistortVec(PetscRandom rctx, Vec v, PetscReal distpar)
00030 {
00031   MagparFunctionInfoBegin;
00032 
00033   Vec VecRnd;
00034   ierr = VecDuplicate(v,&VecRnd);CHKERRQ(ierr);
00035 
00036   /* FIXME: distributions generated by this method are not uniform on the sphere!
00037      see/fix also init/magset.c: init_mag 8
00038    */
00039   ierr = VecSetRandom(VecRnd,rctx);CHKERRQ(ierr);
00040 
00041   ierr = VecScale(VecRnd,2.0);CHKERRQ(ierr);
00042   ierr = VecShift(VecRnd,-1.0);CHKERRQ(ierr);
00043 
00044   ierr = VecAXPY(v,distpar,VecRnd);CHKERRQ(ierr);
00045 
00046   ierr = VecDestroy(VecRnd);CHKERRQ(ierr);
00047 
00048   MagparFunctionInfoReturn(0);
00049 }
00050 

magpar - Parallel Finite Element Micromagnetics Package
Copyright (C) 2002-2010 Werner Scholz