X-Git-Url: https://hackdaworld.org/gitweb/?p=physik%2Fising.git;a=blobdiff_plain;f=ising.c;h=e1444d891d0b818178c60564007fdac64f60c50e;hp=b6c3a8721e8d1bc52d2f96a4d46901827601c016;hb=62002d5d0469045bb0b6f1ff7832ed5029c5267d;hpb=f83bfccf12f96b6b3121a3d8bf1fbc2df1a07b9a

diff --git a/ising.c b/ising.c
index b6c3a87..e1444d8 100644
--- a/ising.c
+++ b/ising.c
@@ -37,6 +37,7 @@ int usage(void)
  puts("-s <value> \t spin interaction strength");
  puts("-r \t run in interactive mode (still in work)");
  puts("-d <value> \t refresh display rate");
+ puts("-t <value> maximal temperature");
 
  return 1;
 }
@@ -75,6 +76,7 @@ int main(int argc, char **argv)
  itt=I;
  s=S;
  strcpy(output_file,"");
+ max_T=0;
  /* parse argv */
  for(i=1;i<argc;i++)
  {
@@ -103,6 +105,9 @@ int main(int argc, char **argv)
     case 'd':
      dr=atoi(argv[++i]);
      break;
+    case 't':
+     max_T=atof(argv[++i]);
+     break;
     default:
      usage();
      return -1;
@@ -127,8 +132,8 @@ int main(int argc, char **argv)
 
  /* begin at T=0 M=1 situation */
  memset(atom,0,max_x*max_y*sizeof(unsigned char));
- 
- max_T=3.0*s;
+
+ if(max_T==0) max_T=3.0*s; 
 
  for(T=.05;T<max_T;T+=.05)
  {