Changeset 4782e78 for pcp/src/ions.c

Timestamp:

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

Author:

Frederik Heber <heber@…>

Children:

6590194

Parents:

9a2c8c

git-author:

Frederik Heber <heber@…> (04/18/08 16:17:36)

git-committer:

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

Message:

CorrectVelocity(): FixedIon may still move, there's just no force acting on them

File:

• pcp/src/ions.c (modified) (3 diffs)

pcp/src/ions.c ¶

@@ -787 +787 @@
 /** Shifts center of gravity of ion velocities (rather momentums).
  * First sums up ion speed U times their IonMass (summed up for each
- * dimension), then reduces the velocity by temp per Ions::TotalMass (so here
+ * dimension) in temp, then reduces the velocity by temp per Ions::TotalMass (so here
  * total number of Ions is included)
  * \param *P Problem at hand
- * \sa CorrectForces
+ * \sa CorrectForces()
  */
 void CorrectVelocity(struct Problem *P)
@@ -803 +803 @@
     for (ia=0; ia < I->I[is].Max_IonsOfType; ia++) {    // ... all ions per type ...
       U = &I->I[is].U[NDIM*ia];
-      if (I->I[is].IMT[ia] == MoveIon) {        // ... if it's movable
+      //if (I->I[is].IMT[ia] == MoveIon) {      // ... even FixedIon moves, only not by other's forces
                                 for (d=0; d<NDIM; d++) 
                                   UTemp[d] += I->I[is].IonMass*U[d];
-      }
+      //}
     }
   }
@@ -812 +812 @@
     for (ia=0; ia < I->I[is].Max_IonsOfType; ia++) {
       U = &I->I[is].U[NDIM*ia];
-      if (I->I[is].IMT[ia] == MoveIon) {
+      //if (I->I[is].IMT[ia] == MoveIon) {
                                 for (d=0; d<NDIM; d++) 
                                   U[d] -= UTemp[d]/I->TotalMass;                // shift by mean velocity over mass and number of ions
-      }
+      //}
     }
   }