# namespace geom::quadrics

    namespace quadrics {

    // global functions

    double f(Eigen::VectorXf Coeff, Eigen::VectorXf Q);
    Eigen::VectorXf nablaf(Eigen::VectorXf Coeff, Eigen::VectorXf Q);
    Eigen::VectorXf SurfacePoint(Eigen::VectorXf P, Eigen::VectorXf Coeff, double epsilon);
    Eigen::VectorXf FootPoint(Eigen::VectorXf P, Eigen::VectorXf Coeff, double epsilon);
    Eigen::VectorXf SolveSystem(Eigen::MatrixXf points, Eigen::MatrixXf normals, Eigen::VectorXf wPoints, Eigen::VectorXf wNormals);

    } // namespace quadrics