void rjMcMC1DSampler::writeresults(FILE* fp)
{
writeheader(fp);
writemodel(fp, mHighestLikelihood);
writemodel(fp, mLowestMisfit);
pmap.writedata(fp);
nmap.writedata(fp);
for (size_t ci = 0; ci < nchains; ci++){
mChainInfo[ci].writeconvergencedata(ci,fp);
mChainInfo[ci].writemodelchain_binary(ci,fp);
}
}
nemo_main()
{
anisorad = getdparam("anisorad");
mcutoff = getdparam("mcutoff");
ntab = MXTB;
plummer();
writemodel(getparam("out"));
}
int main(int argc, string argv[])
{
initparam(argv, defv);
readgsp();
init_random(getiparam("seed"));
layout_body(bodyfields, Precision, NDIM);
gspmodel();
writemodel();
return (0);
}
int main(int argc, string argv[])
{
initparam(argv, defv);
readgsp();
init_random(getiparam("seed"));
layout_body(bodyfields, Precision, NDIM);
gspsphere();
if (getbparam("zerocm"))
snapcenter(btab, nbody, MassField.offset);
writemodel();
return (0);
}
int main(int argc, string argv[])
{
float tmp1, tmp2;
initparam(argv, defv);
if (sscanf(getparam("mdisk"), "%f,%f", &tmp1, &tmp2) == 2) {
mdisk1 = tmp1;
mdisk2 = tmp2;
if (sscanf(getparam("alpha"), "%f,%f", &tmp1, &tmp2) != 2)
error("%s: must specify two alpha values\n", getargv0());
alpha1 = tmp1;
alpha2 = tmp2;
if (sscanf(getparam("epsilon"), "%f,%f", &tmp1, &tmp2) != 2)
error("%s: must specify two epsilon values\n", getargv0());
epsilon1 = tmp1;
epsilon2 = tmp2;
} else {
mdisk1 = getdparam("mdisk");
mdisk2 = 0.0;
alpha1 = getdparam("alpha");
alpha2 = 12.0; /* nonzero value stops div by zero */
epsilon1 = getdparam("epsilon");
epsilon2 = -1.0;
}
zdisk = getdparam("zdisk");
mu = getdparam("mu");
r_mu = getdparam("r_mu");
eta = getdparam("eta");
rcut = getdparam("rcut");
readgsp();
setprof();
layout_body(bodyfields, Precision, NDIM);
ndisk = getiparam("ndisk");
disk = (bodyptr) allocate(ndisk * SizeofBody);
init_random(getiparam("seed"));
if (getiparam("nlist") > 0)
listdisk(getiparam("nlist"));
makedisk();
writemodel();
return (0);
}
int main(int argc, string argv[])
{
double alpha, rcut;
int model;
initparam(argv, defv);
alpha = getdparam("alpha");
rcut = getdparam("rcut");
model = (alpha <= 0.0 ? -1 : getiparam("model"));
ndisk = getiparam("ndisk");
init_random(getiparam("seed"));
readgsp();
setprof(model, alpha, rcut);
layout_body(bodyfields, Precision, NDIM);
disk = (bodyptr) allocate(ndisk * SizeofBody);
makedisk(getbparam("randspin"));
if (! getbparam("randspin")) // if spins not random
bodyfields[2] = NULL; // don't write AuxVec field
writemodel();
fflush(NULL);
return 0;
}
