hexternal.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: hexternal.c 2744 2009-08-26 21:44:29Z scholz $\n\n";
00025 static char const Td[] = "$Today: " __FILE__ " " __DATE__ " " __TIME__ " $\n\n";
00026
00027 #include "field.h"
00028 #include "util/util.h"
00029
00030 #ifdef ADDONS
00031 #include "addons/addons.h"
00032 #endif
00033
00034 #ifdef SPINTORQ
00035 #include "spintorq/spintorq.h"
00036 #endif
00037
00038 static int doinit=1;
00039 static Vec VHextsum=PETSC_NULL;
00040 static PetscReal hextsum=0;
00041
00042 int Hexternal_Init(GridData *gdata)
00043 {
00044 MagparFunctionLogBegin;
00045
00046 PetscPrintf(PETSC_COMM_WORLD,"Hexternal on\n");
00047
00048 ierr = VecDuplicate(gdata->M,&VHextsum);CHKERRQ(ierr);
00049 ierr = VecDuplicate(gdata->M,&gdata->VHext);CHKERRQ(ierr);
00050
00051 MagparFunctionLogReturn(0);
00052 }
00053
00054
00055 int Hexternal(GridData *gdata,Vec VHtotsum)
00056 {
00057 MagparFunctionInfoBegin;
00058
00059
00060 if (doinit) {
00061 ierr = Hexternal_Init(gdata);CHKERRQ(ierr);
00062 doinit=0;
00063 }
00064
00065 hextsum=0;
00066 ierr = VecZeroEntries(VHextsum);CHKERRQ(ierr);
00067
00068 PetscReal hext;
00069 hext=0;ierr = Hext_ho(gdata,VHextsum,&hext);CHKERRQ(ierr); hextsum += hext;
00070 hext=0;ierr = Hext_cu(gdata,VHextsum,&hext);CHKERRQ(ierr); hextsum += hext;
00071 hext=0;ierr = Hext_py(gdata,VHextsum,&hext);CHKERRQ(ierr); hextsum += hext;
00072 #ifdef ADDONS
00073 hext=0;ierr = Hext_ac(gdata,VHextsum,&hext);CHKERRQ(ierr); hextsum += hext;
00074 hext=0;ierr = Hext_co(gdata,VHextsum,&hext);CHKERRQ(ierr); hextsum += hext;
00075 hext=0;ierr = Hext_do(gdata,VHextsum,&hext);CHKERRQ(ierr); hextsum += hext;
00076 hext=0;ierr = Hext_ht1(gdata,VHextsum,&hext);CHKERRQ(ierr);hextsum += hext;
00077 hext=0;ierr = Hext_ht2(gdata,VHextsum,&hext);CHKERRQ(ierr);hextsum += hext;
00078 hext=0;ierr = Hext_mp(gdata,VHextsum,&hext);CHKERRQ(ierr); hextsum += hext;
00079 hext=0;ierr = Hext_pd(gdata,VHextsum,&hext);CHKERRQ(ierr); hextsum += hext;
00080 hext=0;ierr = Hext_rg(gdata,VHextsum,&hext);CHKERRQ(ierr); hextsum += hext;
00081 hext=0;ierr = Hext_tf(gdata,VHextsum,&hext);CHKERRQ(ierr); hextsum += hext;
00082 #endif
00083 #ifdef SPINTORQ
00084 hext=0;ierr = Hcurrent(gdata,VHextsum,&hext);CHKERRQ(ierr);hextsum += hext;
00085 hext=0;ierr = Hcurrent_stack(gdata,VHextsum,&hext);CHKERRQ(ierr);hextsum += hext;
00086 #endif
00087
00088
00089 ierr = VecCopy(VHextsum,gdata->VHext);CHKERRQ(ierr);
00090
00091 ierr = VecAXPY(VHtotsum,1,VHextsum);CHKERRQ(ierr);
00092
00093 MagparFunctionProfReturn(0);
00094 }
00095
00096
00097 int Hexternal_Energy(GridData *gdata,Vec VMom,PetscReal *energy)
00098 {
00099 PetscReal e;
00100
00101 MagparFunctionInfoBegin;
00102
00103 ierr = VecDot(VMom,VHextsum,&e);CHKERRQ(ierr);
00104 e *= -1.0;
00105
00106 *energy=e;
00107
00108 MagparFunctionProfReturn(0);
00109 }
00110
00111
00112 PetscReal Hexternal_hext()
00113 {
00114 MagparFunctionInfoBegin;
00115 MagparFunctionInfoReturn(hextsum);
00116 }
00117