projects
/
physik
/
ising.git
/ commitdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
| commitdiff |
tree
raw
|
patch
|
inline
| side by side (from parent 1:
f83bfcc
)
added max t option
author
hackbard
<hackbard>
Mon, 19 May 2003 01:09:52 +0000
(
01:09
+0000)
committer
hackbard
<hackbard>
Mon, 19 May 2003 01:09:52 +0000
(
01:09
+0000)
ising.c
patch
|
blob
|
history
diff --git
a/ising.c
b/ising.c
index
b6c3a87
..
e1444d8
100644
(file)
--- 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("-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;
}
return 1;
}
@@
-75,6
+76,7
@@
int main(int argc, char **argv)
itt=I;
s=S;
strcpy(output_file,"");
itt=I;
s=S;
strcpy(output_file,"");
+ max_T=0;
/* parse argv */
for(i=1;i<argc;i++)
{
/* 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 'd':
dr=atoi(argv[++i]);
break;
+ case 't':
+ max_T=atof(argv[++i]);
+ break;
default:
usage();
return -1;
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));
/* 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)
{
for(T=.05;T<max_T;T+=.05)
{