sqr
radius = sqrt((float)(sqr(x0 - x) + sqr(y0 - y)));
delta = sqr(Xs + 1)
+ sqr(Ys - 1)
-sqr(xs)
-sqr(ys);
Ys = Yc = sqrt((float)(sqr(xs) + sqr(ys)));
delta = sqr(Xs + 1) + sqr(Ys - 1) - sqr(xs) - sqr(ys);
mp_err mp_sqr(const mp_int *a, mp_int *sqr)
ARGCHK(a != NULL && sqr != NULL, MP_BADARG);
if (a == sqr) {
if (ix > MP_ALLOC(sqr)) {
MP_USED(sqr) = 1;
MP_CHECKOK( s_mp_grow(sqr, ix) );
MP_USED(sqr) = ix;
MP_DIGIT(sqr, 0) = 0;
s_mp_sqr_comba_4(a, sqr);
s_mp_sqr_comba_8(a, sqr);
s_mp_sqr_comba_16(a, sqr);
s_mp_sqr_comba_32(a, sqr);
s_mpv_mul_d(pa, count, d, MP_DIGITS(sqr) + 1);
s_mpv_mul_d_add(pa, count, d, MP_DIGITS(sqr) + ix);
MP_DIGIT(sqr, MP_USED(sqr)-1) = 0; /* above loop stopped short of this. */
s_mp_mul_2(sqr);
MP_DIGIT(sqr, 1) = 0;
s_mpv_sqr_add_prop(MP_DIGITS(a), MP_USED(a), MP_DIGITS(sqr));
SIGN(sqr) = ZPOS;
s_mp_clamp(sqr);