int main(int, char **) {
FILE *fout = fopen("papers/fuzzy-fmt/figs/wallsfillingfracInfo.txt", "w");
fclose(fout);
const double etas[] = { 0.1, 0.4 };
const double temps[] = { 0.0, 0.01, 0.02, 0.03 };
for (double eta = 0.4; eta > 0; eta-=0.3) {
for (unsigned int i = 0; i<sizeof(temps)/sizeof(temps[0]); i++) {
const double temp = temps[i];
Functional f = HardFluid(1,0);
if (temp > 0) f = SoftFluid(1, 1, 0);
const double mu = find_chemical_potential(OfEffectivePotential(f), (temp)?temp:1, eta/(4*M_PI/3));
printf("mu is %g\n", mu);
if (temp > 0) f = SoftFluid(1, 1, mu);
else f = HardFluid(1, mu);
const char *name = "hard";
if (temp > 0) name = "soft";
run_walls(eta, name, f, temp);
}
}
// Just create this file so make knows we have run.
if (!fopen("papers/fuzzy-fmt/figs/walls.dat", "w")) {
printf("Error creating walls.dat!\n");
return 1;
}
return 0;
}
int main(int, char **) {
FILE *out = fopen("papers/fuzzy-fmt/figs/bh-homogeneous.dat", "w");
const double Tmax = 10.0, dT = 0.01, Tmin = dT;
fprintf(out, "# n_reduced");
for (double T = Tmin; T<= Tmax + dT/2; T += dT) {
fprintf(out, "\tp(kT=%g)/nkT", T);
}
fprintf(out, "\n0");
for (double T = Tmin; T<= Tmax + dT/2; T += dT) {
fprintf(out, "\t%g", T);
}
fprintf(out, "\n");
const double dn = 0.01, nmax = 2.5;
for (double n_reduced = dn; n_reduced <= nmax; n_reduced += dn) {
fprintf(out, "%g", n_reduced);
for (double T = Tmin; T<= Tmax + dT/2; T += dT) {
const double temp = T;
double rad = R_BH(temp);
Functional f = HardFluid(rad,0);
double usekT = temp;
//if (temp == 0) usekT = 1.0;
const double n = n_reduced*pow(2,-5.0/2.0);
// return *reduced* pressure!
fprintf(out, "\t%g", pressure(OfEffectivePotential(f), usekT, n)/pow(2,-5.0/2.0));
}
fprintf(out, "\n");
}
fclose(out);
return 0;
}