int lindholm(double g, double w, double x, double y, double z,
    double *Hx, double *Hy, double *Hz);