Program: hc2
Copyright (C) by Sherryl Manalo
This program is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 2 of the License, or (at your option) any later version.
This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
Contact the author at: sm@antitachyon.com
The program hc2 calculates the critical field as a function of temperature for a certain external field angle.
----------------------- compile ----------------------
To compile the program:
make hc2
----------------------- run --------------------------
To run the program:
hc2 hc2.cfg temperaturefile outputname
hc2.cfg:config file temperaturefile: one-column file with temperature points, preferably from high temperatures to lower i.e. dtemp.conf outputname: main name of files in which the data are stored outputname.lst, outputname.orig (outputname.gap if calcgap=1)
---------------------- hc2.cfg ---------------------- rm: parameter rm
mod: 0 - s-wave
1 - d-wave
phi: in-plane angle of field (for d-wave)
hkapt: starting value of field
ntemp: number of temperatures to iterate
Q-num: quantum number (s-wave),
starting quantum number (d-wave - 0,1,2,3)
Dim.-D: dimension of D-matrix (d-wave)
Max.-iter: maximum number of iterations
tolerance: tolerance of iteration
convergence: 0 if everything's normal, 1 if problems with conv.
slow: 0 if everything's normal, 1 if problems with conv.
sdamp: damping factor of iteration
stable: stabilizing factor < 1, use if slow = 1
up: 1 - add a value to next fieldvalue to iterate
0 - don't ...
upsize: this value is added
full: 0 - normal
TOLE: tolerance of integration
maxit-integ: maximum number of iterations in integration
calcgapv: 0 - calculate H_crit vs t
1 - calculate Gap-eq. vs. H (Gap-eq. = 0 => H_crit)
rangeh-max: maximum numbers of mesh points for calcgapv
stepsize: stepsize for calcgapv
doit: 0 - normal calculation (one iteration of hcomDelt)
1 - in case of problems (more than 1 iteration, see below)
domax: number of iterations of hcomDelt
range: takes effect if calcgapv = 1
gapread: 1 - read gap-values from file
0 - don't ..
compl: 1 - print gap-values in .lst-file (this could get long)
0 - don't ..