Particle Trajectory Code - Revision history

Oborn at 22:29, 17 June 2009

2009-06-17T22:29:07Z

2009-06-17T22:28:50Z

2009-06-11T20:05:24Z

2009-06-10T21:44:04Z

2009-06-10T21:42:47Z

2009-06-10T21:42:35Z

2009-06-10T21:40:36Z

2009-06-10T21:38:19Z

2009-06-10T21:37:50Z

2009-06-10T21:35:18Z

← Older revision		Revision as of 22:29, 17 June 2009
Line 1:		Line 1:
−	This is the old version of the code, which is now obsolete. The new version can be found [[New Particle Trajectory Code\|here]].	+	'''This is the old version of the code, which is now obsolete.''' The new version can be found [[New Particle Trajectory Code\|here]].

	For some reason I am unable to load the source code files onto the wiki, so I have just copied and pasted them below. The boldface headers give the appropriate name for each file. You will need all of the files in order to compile the program.		For some reason I am unable to load the source code files onto the wiki, so I have just copied and pasted them below. The boldface headers give the appropriate name for each file. You will need all of the files in order to compile the program.

← Older revision		Revision as of 22:28, 17 June 2009
Line 1:		Line 1:
		+	This is the old version of the code, which is now obsolete. The new version can be found [[New Particle Trajectory Code\|here]].
		+
	For some reason I am unable to load the source code files onto the wiki, so I have just copied and pasted them below. The boldface headers give the appropriate name for each file. You will need all of the files in order to compile the program.		For some reason I am unable to load the source code files onto the wiki, so I have just copied and pasted them below. The boldface headers give the appropriate name for each file. You will need all of the files in order to compile the program.

@@ Line 482: / Line 482: @@
          /* Break if electron hits the yoke */
-	if(ey1<=0.0 && ey2>=0.0 && ex1>ex2)
+	if(ey2<=0.0 && ey1>0.0 && ex2>0.0 && ex2<LMAGNET)
+	{
 	  flag=0;
@@ Line 494: / Line 494: @@
        /* Check if the electron ever made it out of the magnet */
-       if(ey1>-0.05 && ey1<HMAGNET+0.05 && ex1>-0.05 && ex1<LMAGNET+0.05)
+       if(flag!=0)
-        flag=5;
        /* Free the solver structures so they can be reused */
@@ Line 566: / Line 569: @@
        /* Check if the positron ever made it out of the magnet */
-       if(ey1>-0.05 && ey1<HMAGNET+0.05 && ex1>-0.05 && ex1<LMAGNET+0.05)
+       if(flag!=0)
-        flag=6;
        /* Free these variables so they can be reused later on */

@@ Line 1: / Line 1: @@
 For some reason I am unable to load the source code files onto the wiki, so I have just copied and pasted them below. The boldface headers give the appropriate name for each file. You will need all of the files in order to compile the program.
 ----
@@ Line 700: / Line 702: @@
 #endif
 </nowiki></pre>

@@ Line 4: / Line 4: @@
 '''*** build.sh ***'''
-<pre><nowiki><highlightSyntax>
+<pre><nowiki>
 #!/bin/sh
@@ Line 12: / Line 12: @@
 gcc -g -O2 -Wall -I. -o check_field fields.c check_field.c -lm
-</highlightSyntax></nowiki></pre>
+</nowiki></pre>
 ----

@@ Line 4: / Line 4: @@
 '''*** build.sh ***'''
-<nowiki>
+<pre><nowiki>
 #!/bin/sh
@@ Line 12: / Line 12: @@
 gcc -g -O2 -Wall -I. -o check_field fields.c check_field.c -lm
-</nowiki>
+</nowiki></pre>
 ----
 '''*** fields.c ***'''
 #include <trajectory.h>
 #include <math.h>
@@ Line 66: / Line 67: @@
    return field;
+}
 ----
 '''*** func.c ***'''
 #include <trajectory.h>
 #include <gsl/gsl_errno.h>
@@ Line 103: / Line 106: @@
    return GSL_SUCCESS;
+}
 ----
 '''*** init.c ***'''
 #include <stdio.h>
 #include <stdlib.h>
@@ Line 188: / Line 193: @@
    return 0;
+}
 ----
 '''*** readinput.c ***'''
 #include <stdio.h>
 #include <stdlib.h>
@@ Line 242: / Line 249: @@
    pars->eelectron=pars->eelectron*1.602E-13; /* Convert to Joules */
+}
 ----
 '''*** trajectory.c ***'''
 /* trajectory.c
   * A program to simulate the paths of an electron-positron pair
@@ Line 633: / Line 642: @@
    return 0;
+}
 ----
 '''*** trajectory.h ***'''
 #ifndef _TRAJECTORY_H_
 #define _TRAJECTORY_H_
@@ Line 688: / Line 699: @@
 #endif