Changeset 32 in TRACY3 for trunk/tracy/tracy/src/mathlib.cc

Timestamp:

Apr 9, 2014, 3:50:11 PM (10 years ago)

Author:

zhangj

Message:

active the transport of the twiss functions and orbits of the transfer line.

File:

• trunk/tracy/tracy/src/mathlib.cc (modified) (4 diffs)

trunk/tracy/tracy/src/mathlib.cc

@@ -218 +218 @@
 }
 
-
+/***********************************************************
+ * void CopyVec(const int n, const Vector &a, Vector &b)
+ * 
+ *  Copy 6D vector a to b: b = a. 
+ * 
+ ***********************************************************/
 void CopyVec(const int n, const Vector &a, Vector &b)
 {
-  int i;
+  int i=0;
 
   for (i = 0; i < n; i++)
@@ -457 +462 @@
 
 
+/*********************************************************
+ * void MulLMat(const int n, const Matrix &A, Matrix &B)
+ * 
+ * Matrix multiplication.  B = A * B
+ * 
+ * 
+ **********************************************************/
 void MulLMat(const int n, const Matrix &A, Matrix &B)
 {
-  int i, j, k;
-  double x;
+  int i=0, j=0, k=0;
+  double x=0.0;
   Matrix C;
 
@@ -474 +486 @@
 }
 
-
+/*********************************************************
+ * void MulRMat(const int n, const Matrix &A, Matrix &B)
+ * 
+ * Matrix multiplication.  A = A * B
+ * 
+ * 
+ **********************************************************/
 void MulRMat(const int n, Matrix &A, const Matrix &B)
 {
-  int i, j, k;
-  double x;
+  int i=0, j=0, k=0;
+  double x=0.0;
   Matrix C;
 
@@ -492 +510 @@
 }
 
-
+/*********************************************************
+ * double TrMat(const int n, const Matrix &A)
+ * 
+ * Calculate the trance of the matrix A.
+ * 
+ *  return:
+ *          x.
+ * 
+ **********************************************************/
 double TrMat(const int n, const Matrix &A)
 {