123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189 |
- #include "Curvepoint.hpp"
- Curvepoint::Curvepoint()
- {
- curvepoint_fp_setneutral(point);
- }
- Curvepoint::Curvepoint(curvepoint_fp_t input)
- {
- curvepoint_fp_set(point, input);
- }
- Curvepoint Curvepoint::operator+(const Curvepoint& b) const
- {
- Curvepoint retval;
- if (*this == b)
- curvepoint_fp_double(retval.point, point);
- else
- curvepoint_fp_add_vartime(retval.point, point, b.point);
- return retval;
- }
- Curvepoint Curvepoint::operator*(const Scalar& exp) const
- {
- Curvepoint retval;
- exp.mult(retval.point, point);
- return retval;
- }
- bool Curvepoint::operator==(const Curvepoint& b) const
- {
- bool retval;
- curvepoint_fp_t affine_this_point, affine_b_point;
- curvepoint_fp_set(affine_this_point, point);
- curvepoint_fp_set(affine_b_point, b.point);
- if (!(fpe_isone(affine_this_point->m_z) || fpe_iszero(affine_this_point->m_z)))
- curvepoint_fp_makeaffine(affine_this_point);
- if (!(fpe_isone(affine_b_point->m_z) || fpe_iszero(affine_b_point->m_z)))
- curvepoint_fp_makeaffine(affine_b_point);
- retval = fpe_iseq(affine_this_point->m_x, affine_b_point->m_x);
- retval = retval && fpe_iseq(affine_this_point->m_y, affine_b_point->m_y);
- retval = retval || (fpe_iszero(affine_this_point->m_z) && fpe_iszero(affine_b_point->m_z));
- return retval;
- }
- bool Curvepoint::operator<(const Curvepoint& b) const
- {
- bool lessThan[2];
- bool equal[2];
- curvepoint_fp_t affine_this_point, affine_b_point;
- lessThan[0] = lessThan[1] = false;
- curvepoint_fp_set(affine_this_point, point);
- curvepoint_fp_set(affine_b_point, b.point);
- if (fpe_iszero(affine_this_point->m_z))
- {
- // this case would be equal
- if (fpe_iszero(affine_b_point->m_z))
- return false;
- // point at infinity is less than all other points
- return true;
- }
- if (fpe_iszero(affine_b_point->m_z))
- return false;
- // already checked for the point at infinity, so we don't have to redo that here
- if (!fpe_isone(affine_this_point->m_z))
- curvepoint_fp_makeaffine(affine_this_point);
- if (!fpe_isone(affine_b_point->m_z))
- curvepoint_fp_makeaffine(affine_b_point);
- for (int i = 11; i >= 0; i--)
- {
- if (affine_this_point->m_x->v[i] > affine_b_point->m_x->v[i])
- {
- lessThan[0] = false;
- break;
- }
- if (affine_this_point->m_x->v[i] < affine_b_point->m_x->v[i])
- {
- lessThan[0] = true;
- break;
- }
- }
- for (int i = 11; i >= 0; i--)
- {
- if (affine_this_point->m_y->v[i] > affine_b_point->m_y->v[i])
- {
- lessThan[1] = false;
- break;
- }
- if (affine_this_point->m_y->v[i] < affine_b_point->m_y->v[i])
- {
- lessThan[1] = true;
- break;
- }
- }
- equal[0] = fpe_iseq(affine_this_point->m_x, affine_b_point->m_x);
- equal[1] = fpe_iseq(affine_this_point->m_y, affine_b_point->m_y);
- // sort is lesser x value first, and then lesser y value second if x's are equal
- return equal[0] ? (equal[1] ? false : lessThan[1]) : lessThan[0];
- }
- bool Curvepoint::operator>(const Curvepoint& b) const
- {
- return !(*this < b);
- }
- bool Curvepoint::operator<=(const Curvepoint& b) const
- {
- return (*this == b) || (*this < b);
- }
- bool Curvepoint::operator>=(const Curvepoint& b) const
- {
- return (*this == b) || !(*this < b);
- }
- bool Curvepoint::operator!=(const Curvepoint& b) const
- {
- return !(*this == b);
- }
- void Curvepoint::make_affine()
- {
- if (!(fpe_isone(point->m_z) || fpe_iszero(point->m_z)))
- curvepoint_fp_makeaffine(point);
- }
- std::ostream& operator<<(std::ostream& os, const Curvepoint& output)
- {
- Curvepoint affine_out = output;
- affine_out.make_affine();
-
- os << Fpe(affine_out.point->m_x) << Fpe(affine_out.point->m_y) << Fpe(affine_out.point->m_z);
- return os;
- }
- std::istream& operator>>(std::istream& is, Curvepoint& input)
- {
- Fpe x, y, z;
- is >> x >> y >> z;
- fpe_set(input.point->m_x, x.data);
- fpe_set(input.point->m_y, y.data);
- fpe_set(input.point->m_z, z.data);
- return is;
- }
- size_t CurvepointHash::operator()(const Curvepoint& x) const
- {
- if (fpe_iszero(x.point->m_z))
- {
- return 0;
- }
- size_t retval;
- std::hash<double> hasher;
- Curvepoint affine_x = x;
- affine_x.make_affine();
- for (int j = 0; j < 12; j++)
- {
- retval ^= hasher(affine_x.point->m_x->v[j]);
- retval ^= hasher(affine_x.point->m_y->v[j]);
- }
- return retval;
- }
|