Changeset 13d7f6

Timestamp:

Apr 21, 2008, 2:19:23 PM (17 years ago)

Author:

Frederik Heber <heber@…>

Children:

460e95

Parents:

9b2f5d

git-author:

Frederik Heber <heber@…> (04/18/08 14:16:44)

git-committer:

Frederik Heber <heber@…> (04/21/08 14:19:23)

Message:

ResetForces() added

File:

• pcp/src/ions.c (modified) (1 diff)

pcp/src/ions.c

@@ -579 +579 @@
 }
 
+/** Resets the forces array.
+ * \param *P Problem at hand
+ */
+void ResetForces(struct Problem *P)
+{
+  struct Ions *I = &P->Ion;
+  double *FIon, *FIonL, *FIonNL, *FEwald;
+  int is, ia, d;
+  for (is=0; is < I->Max_Types; is++) { // go through all types ...
+    for (ia=0; ia < I->I[is].Max_IonsOfType; ia++) { // .. and each ion of the type
+      FIon = &I->I[is].FIon[NDIM*ia];
+      FIonL = &I->I[is].FIonL[NDIM*ia];
+      FIonNL = &I->I[is].FIonNL[NDIM*ia];
+      FEwald = &I->I[is].FEwald[NDIM*ia];
+      for (d=0; d<NDIM; d++) {
+        FIon[d] = 0.;                   // zero all as a sign that's been moved
+        FIonL[d] = 0.;
+        FIonNL[d] = 0.;
+        FEwald[d] = 0.;
+      }
+    }
+  }
+};
+
 /** Changes ion's velocity according to acting force.
  * IonType::U is changed by the weighted force of actual step and last one