geodesic.h

// -*-c++-*-
#ifndef STORMM_GEODESIC_H
#define STORMM_GEODESIC_H
 
#include <vector>
#include "copyright.h"
#include "Constants/symbol_values.h"
#include "math_enumerators.h"
#include "matrix_ops.h"
 
namespace stormm {
namespace stmath {
 
using symbols::pi;
using symbols::twopi;
using symbols::tetrahedral_angle;

template <typename T, typename T3>
int gatherTriangles(int midpt_llim, int midpt_hlim, int ring_llim, int ring_hlim, int nring,
                    T crit_rsq, std::vector<T3> *result, int nadd_pt,
                    const std::vector<double3> &scaffold);

template <typename T, typename T3>
std::vector<T3> surfaceDistribution(int *n, SpherePlacement method = SpherePlacement::POLYHEDRON);
 
template <typename T, typename T3>
std::vector<T3> surfaceDistribution(int n, SpherePlacement method = SpherePlacement::POLYHEDRON);
 
} // namespace stmath
} // namespace stormm
 
#include "geodesic.tpp"
 
#endif