restraint_util.h

// -*-c++-*-
#ifndef STORMM_RESTRAINT_UTIL_H
#define STORMM_RESTRAINT_UTIL_H
 
#include "copyright.h"
#include "DataTypes/stormm_vector_types.h"
 
namespace stormm {
namespace restraints {

template <typename T>
Vec3<T> restraintDelta(const Vec2<T> init_k, const Vec2<T> final_k, const Vec4<T> init_r,
                       const Vec4<T> final_r, const Vec2<T> mixwt, T dr);

template <typename T>
Vec2<T> computeRestraintMixture(int step_number, int init_step, int final_step);
 
} // namespace restraints
} // namespace stormm
 
#include "restraint_util.tpp"
 
#endif