static int
adwrite_fir (long stageid, double gnom, Fir * fir)
{
int i;
char *desc = "",
symmetry;
double group_correction = 0.0;
int n2print;
int errors = 0;
get_symmetry (fir, &n2print, &symmetry);
fprintf (stdout, "FIR2 %2ld %10.2e %4d %8.3f %c %4d %-25.25s\n",
stageid, gnom, fir->dec_factor,
group_correction, symmetry, n2print, desc);
if (fir->nden > 0) {
complain (0, "fir filter has denominator coefficients, autoDRM doesn't allow");
errors++;
}
for (i = 0; i < n2print; i++) {
if (i > 0 && (i % 5) == 0)
fprintf (stdout, "\n");
fprintf (stdout, " %15.8e",
fir->num_coefs[i]);
}
fprintf (stdout, "\n");
return errors;
}
int spg_get_symmetry(int rotation[][3][3],
double translation[][3],
const int max_size,
SPGCONST double lattice[3][3],
SPGCONST double position[][3],
const int types[],
const int num_atom,
const double symprec)
{
sym_set_angle_tolerance(-1.0);
return get_symmetry(rotation,
translation,
max_size,
lattice,
position,
types,
num_atom,
symprec);
}