Hello Chris and all,
The holding field routines are now available as promised in the meeting.
You may be able to check out at exp/analysis/utils. You may be able to
retrieve three files:
HoldField_Tosca.h ... A class definition file which interpolates
the holding field from a TOSCA field map.
HoldField_BiotSavart.h ... A class definition file which calculates the
holding field completely independent from HoldField_Tosca.h
show_hf.C ... a sample usage of the classes and
draws 2d field graph from one of the methods.
HoldField_Tosca.h actually needs two tosca table files (8.1MB each),
which are located at /scratch/bud18/shino/magnet/tosca/*.table. The
program use a "Liner Tetrahedral method," which guarantees to be
continuous but not smooth, unlike "SUBROUTINE SNAK" in blast_proc.f,
which uses a Taylor expansion up to 2nd order, which is smooth but not
actually continuous (a bit of a problem).
HoldField_BiotSavart.h is a very compact and fast theoretical package
and due to a semi phenomenological model. The difference from the TOSCA
is within ~20 Gauss. There is no boundary limit for the computation.
Thank you,
Aki
This archive was generated by hypermail 2.1.2 : Mon Feb 24 2014 - 14:07:32 EST