cimport cposit
-# special values needed for some initializations
+# special values and C helpers
-cdef posit8_one = cposit.ui32_to_p8(1)
-cdef posit16_one = cposit.ui32_to_p16(1)
+cdef _p8_one = cposit.ui32_to_p8(1)
+cdef _p16_one = cposit.ui32_to_p16(1)
+
+cdef inline cposit.posit8_t _p8_neg(cposit.posit8_t f):
+ f.v = -f.v
+ return f
+
+cdef inline cposit.posit16_t _p16_neg(cposit.posit16_t f):
+ f.v = -f.v
+ return f
+
+cdef inline cposit.posit8_t _p8_abs(cposit.posit8_t f):
+ f.v = <uint8_t> abs(<int8_t> f.v)
+ return f
+
+cdef inline cposit.posit16_t _p16_abs(cposit.posit16_t f):
+ f.v = <uint16_t> abs(<int16_t> f.v)
+ return f
cdef class Posit8:
# arithmetic
+ cpdef Posit8 neg(self):
+ cdef cposit.posit8_t f = _p8_neg(self._c_posit)
+ return Posit8.from_c_posit(f)
+
+ def __neg__(self):
+ return self.neg()
+
+ cpdef Posit8 abs(self):
+ cdef cposit.posit8_t f = _p8_abs(self._c_posit)
+ return Posit8.from_c_posit(f)
+
+ def __abs__(self):
+ return self.abs()
+
cpdef Posit8 round(self):
cdef cposit.posit8_t f = cposit.p8_roundToInt(self._c_posit)
return Posit8.from_c_posit(f)
def __round__(self):
- return self.roundToInt()
+ return self.round()
cpdef Posit8 add(self, Posit8 other):
cdef cposit.posit8_t f = cposit.p8_add(self._c_posit, other._c_posit)
cdef cposit.posit8_t f = cposit.p8_mulAdd(self._c_posit, a2._c_posit, a3._c_posit)
return Posit8.from_c_posit(f)
- cpdef Posit8 qma(self, Posit8 a1, Posit8 a2):
+ cpdef Posit8 fam(self, Posit8 a1, Posit8 a2):
cdef cposit.posit8_t f = cposit.p8_mulAdd(a1._c_posit, a2._c_posit, self._c_posit)
return Posit8.from_c_posit(f)
# in-place arithmetic
+ cpdef void ineg(self):
+ self._c_posit = _p8_neg(self._c_posit)
+
+ cpdef void iabs(self):
+ self._c_posit = _p8_abs(self._c_posit)
+
cpdef void iround(self):
self._c_posit = cposit.p8_roundToInt(self._c_posit)
cpdef void ifma(self, Posit8 a2, Posit8 a3):
self._c_posit = cposit.p8_mulAdd(self._c_posit, a2._c_posit, a3._c_posit)
- cpdef void iqma(self, Posit8 a1, Posit8 a2):
+ cpdef void ifam(self, Posit8 a1, Posit8 a2):
self._c_posit = cposit.p8_mulAdd(a1._c_posit, a2._c_posit, self._c_posit)
cpdef void idiv(self, Posit8 other):
cpdef to_quire(self):
cdef cposit.quire8_t f
cposit.q8_clr(f)
- f = cposit.q8_fdp_add(f, self._c_posit, posit8_one)
+ f = cposit.q8_fdp_add(f, self._c_posit, _p8_one)
return Quire8.from_c_quire(f)
else:
f = float(value)
cposit.q8_clr(self._c_quire)
- self._c_quire = cposit.q8_fdp_add(self._c_quire, cposit.convertDoubleToP8(f), posit8_one)
+ self._c_quire = cposit.q8_fdp_add(self._c_quire, cposit.convertDoubleToP8(f), _p8_one)
def __float__(self):
return cposit.convertP8ToDouble(cposit.q8_to_p8(self._c_quire))
# arithmetic
- cpdef Quire8 qma(self, Posit8 a1, Posit8 a2):
+ cpdef Quire8 qam(self, Posit8 a1, Posit8 a2):
cdef cposit.quire8_t f = cposit.q8_fdp_add(self._c_quire, a1._c_posit, a2._c_posit)
return Quire8.from_c_quire(f)
- cpdef Quire8 qms(self, Posit8 a1, Posit8 a2):
+ cpdef Quire8 qsm(self, Posit8 a1, Posit8 a2):
cdef cposit.quire8_t f = cposit.q8_fdp_sub(self._c_quire, a1._c_posit, a2._c_posit)
return Quire8.from_c_quire(f)
- cpdef void iqma(self, Posit8 a1, Posit8 a2):
+ cpdef void iqam(self, Posit8 a1, Posit8 a2):
self._c_quire = cposit.q8_fdp_add(self._c_quire, a1._c_posit, a2._c_posit)
- cpdef void iqms(self, Posit8 a1, Posit8 a2):
+ cpdef void iqsm(self, Posit8 a1, Posit8 a2):
self._c_quire = cposit.q8_fdp_sub(self._c_quire, a1._c_posit, a2._c_posit)
cpdef void iclr(self):
# external, non-method arithmetic
+cpdef Posit8 p8_neg(Posit8 a1):
+ cdef cposit.posit8_t f = _p8_neg(a1._c_posit)
+ return Posit8.from_c_posit(f)
+
+cpdef Posit8 p8_abs(Posit8 a1):
+ cdef cposit.posit8_t f = _p8_abs(a1._c_posit)
+ return Posit8.from_c_posit(f)
+
cpdef Posit8 p8_round(Posit8 a1):
cdef cposit.posit8_t f = cposit.p8_roundToInt(a1._c_posit)
return Posit8.from_c_posit(f)
cdef cposit.posit8_t f = cposit.p8_mulAdd(a1._c_posit, a2._c_posit, a3._c_posit)
return Posit8.from_c_posit(f)
-cpdef Posit8 p8_qma(Posit8 a3, Posit8 a1, Posit8 a2):
+cpdef Posit8 p8_fam(Posit8 a3, Posit8 a1, Posit8 a2):
cdef cposit.posit8_t f = cposit.p8_mulAdd(a1._c_posit, a2._c_posit, a3._c_posit)
return Posit8.from_c_posit(f)
cpdef Quire8 p8_to_q8(Posit8 a1):
cdef cposit.quire8_t f
cposit.q8_clr(f)
- f = cposit.q8_fdp_add(f, a1._c_posit, posit8_one)
+ f = cposit.q8_fdp_add(f, a1._c_posit, _p8_one)
return Quire8.from_c_quire(f)
-cpdef Quire8 q8_qma(Quire8 a3, Posit8 a1, Posit8 a2):
+cpdef Quire8 q8_qam(Quire8 a3, Posit8 a1, Posit8 a2):
cdef cposit.quire8_t f = cposit.q8_fdp_add(a3._c_quire, a1._c_posit, a2._c_posit)
return Quire8.from_c_quire(f)
-cpdef Quire8 q8_qms(Quire8 a3, Posit8 a1, Posit8 a2):
+cpdef Quire8 q8_qsm(Quire8 a3, Posit8 a1, Posit8 a2):
cdef cposit.quire8_t f = cposit.q8_fdp_sub(a3._c_quire, a1._c_posit, a2._c_posit)
return Quire8.from_c_quire(f)
# arithmetic
+ cpdef Posit16 neg(self):
+ cdef cposit.posit16_t f = _p16_neg(self._c_posit)
+ return Posit16.from_c_posit(f)
+
+ def __neg__(self):
+ return self.neg()
+
+ cpdef Posit16 abs(self):
+ cdef cposit.posit16_t f = _p16_abs(self._c_posit)
+ return Posit16.from_c_posit(f)
+
+ def __abs__(self):
+ return self.abs()
+
cpdef Posit16 round(self):
cdef cposit.posit16_t f = cposit.p16_roundToInt(self._c_posit)
return Posit16.from_c_posit(f)
def __round__(self):
- return self.roundToInt()
+ return self.round()
cpdef Posit16 add(self, Posit16 other):
cdef cposit.posit16_t f = cposit.p16_add(self._c_posit, other._c_posit)
cdef cposit.posit16_t f = cposit.p16_mulAdd(self._c_posit, a2._c_posit, a3._c_posit)
return Posit16.from_c_posit(f)
- cpdef Posit16 qma(self, Posit16 a1, Posit16 a2):
+ cpdef Posit16 fam(self, Posit16 a1, Posit16 a2):
cdef cposit.posit16_t f = cposit.p16_mulAdd(a1._c_posit, a2._c_posit, self._c_posit)
return Posit16.from_c_posit(f)
# in-place arithmetic
+ cpdef void ineg(self):
+ self._c_posit = _p16_neg(self._c_posit)
+
+ cpdef void iabs(self):
+ self._c_posit = _p16_abs(self._c_posit)
+
cpdef void iround(self):
self._c_posit = cposit.p16_roundToInt(self._c_posit)
cpdef void ifma(self, Posit16 a2, Posit16 a3):
self._c_posit = cposit.p16_mulAdd(self._c_posit, a2._c_posit, a3._c_posit)
- cpdef void iqma(self, Posit16 a1, Posit16 a2):
+ cpdef void ifam(self, Posit16 a1, Posit16 a2):
self._c_posit = cposit.p16_mulAdd(a1._c_posit, a2._c_posit, self._c_posit)
cpdef void idiv(self, Posit16 other):
cpdef to_quire(self):
cdef cposit.quire16_t f
cposit.q16_clr(f)
- f = cposit.q16_fdp_add(f, self._c_posit, posit16_one)
+ f = cposit.q16_fdp_add(f, self._c_posit, _p16_one)
return Quire16.from_c_quire(f)
else:
f = float(value)
cposit.q16_clr(self._c_quire)
- self._c_quire = cposit.q16_fdp_add(self._c_quire, cposit.convertDoubleToP16(f), posit16_one)
+ self._c_quire = cposit.q16_fdp_add(self._c_quire, cposit.convertDoubleToP16(f), _p16_one)
def __float__(self):
return cposit.convertP16ToDouble(cposit.q16_to_p16(self._c_quire))
# arithmetic
- cpdef Quire16 qma(self, Posit16 a1, Posit16 a2):
+ cpdef Quire16 qam(self, Posit16 a1, Posit16 a2):
cdef cposit.quire16_t f = cposit.q16_fdp_add(self._c_quire, a1._c_posit, a2._c_posit)
return Quire16.from_c_quire(f)
- cpdef Quire16 qms(self, Posit16 a1, Posit16 a2):
+ cpdef Quire16 qsm(self, Posit16 a1, Posit16 a2):
cdef cposit.quire16_t f = cposit.q16_fdp_sub(self._c_quire, a1._c_posit, a2._c_posit)
return Quire16.from_c_quire(f)
- cpdef void iqma(self, Posit16 a1, Posit16 a2):
+ cpdef void iqam(self, Posit16 a1, Posit16 a2):
self._c_quire = cposit.q16_fdp_add(self._c_quire, a1._c_posit, a2._c_posit)
- cpdef void iqms(self, Posit16 a1, Posit16 a2):
+ cpdef void iqsm(self, Posit16 a1, Posit16 a2):
self._c_quire = cposit.q16_fdp_sub(self._c_quire, a1._c_posit, a2._c_posit)
cpdef void iclr(self):
# external, non-method arithmetic
+cpdef Posit16 p16_neg(Posit16 a1):
+ cdef cposit.posit16_t f = _p16_neg(a1._c_posit)
+ return Posit16.from_c_posit(f)
+
+cpdef Posit16 p16_abs(Posit16 a1):
+ cdef cposit.posit16_t f = _p16_abs(a1._c_posit)
+ return Posit16.from_c_posit(f)
+
cpdef Posit16 p16_round(Posit16 a1):
cdef cposit.posit16_t f = cposit.p16_roundToInt(a1._c_posit)
return Posit16.from_c_posit(f)
cdef cposit.posit16_t f = cposit.p16_mulAdd(a1._c_posit, a2._c_posit, a3._c_posit)
return Posit16.from_c_posit(f)
-cpdef Posit16 p16_qma(Posit16 a3, Posit16 a1, Posit16 a2):
+cpdef Posit16 p16_fam(Posit16 a3, Posit16 a1, Posit16 a2):
cdef cposit.posit16_t f = cposit.p16_mulAdd(a1._c_posit, a2._c_posit, a3._c_posit)
return Posit16.from_c_posit(f)
cpdef Quire16 p16_to_q16(Posit16 a1):
cdef cposit.quire16_t f
cposit.q16_clr(f)
- f = cposit.q16_fdp_add(f, a1._c_posit, posit16_one)
+ f = cposit.q16_fdp_add(f, a1._c_posit, _p16_one)
return Quire16.from_c_quire(f)
-cpdef Quire16 q16_qma(Quire16 a3, Posit16 a1, Posit16 a2):
+cpdef Quire16 q16_qam(Quire16 a3, Posit16 a1, Posit16 a2):
cdef cposit.quire16_t f = cposit.q16_fdp_add(a3._c_quire, a1._c_posit, a2._c_posit)
return Quire16.from_c_quire(f)
-cpdef Quire16 q16_qms(Quire16 a3, Posit16 a1, Posit16 a2):
+cpdef Quire16 q16_qsm(Quire16 a3, Posit16 a1, Posit16 a2):
cdef cposit.quire16_t f = cposit.q16_fdp_sub(a3._c_quire, a1._c_posit, a2._c_posit)
return Quire16.from_c_quire(f)