67    int mLen = MOD_LEN(pME);  in alm_mont_inv()  local
73    BNU_CHUNK_T* ps = pu+mLen;  in alm_mont_inv()
74    BNU_CHUNK_T* pv = ps+mLen;  in alm_mont_inv()
75    BNU_CHUNK_T* pt = pv+mLen;  in alm_mont_inv()
83    COPY_BNU(pu, pm, mLen);  in alm_mont_inv()
84    ZEXPAND_BNU(ps, 0, mLen); ps[0] = 1;  in alm_mont_inv()
85    COPY_BNU(pv, pa, mLen);  in alm_mont_inv()
86    ZEXPAND_BNU(pt, 0, mLen);  in alm_mont_inv()
88    while(!cpEqu_BNU_CHUNK(pv, mLen, 0)) {             // while(v>0) {  in alm_mont_inv()
90          cpLSR_BNU(pu, pu, mLen, 1);                  //       u = u/2;  in alm_mont_inv()
91          cpAdd_BNU(ps, ps, ps, mLen);                 //       s = 2*s;  in alm_mont_inv()
94          cpLSR_BNU(pv, pv, mLen, 1);                  //       v = v/2;  in alm_mont_inv()
95          /*ext +=*/ cpAdd_BNU(pt, pt, pt, mLen);      //       t = 2*t;  in alm_mont_inv()
98          int cmpRes = cpCmp_BNU(pu, mLen, pv, mLen);  in alm_mont_inv()
100             cpSub_BNU(pu, pu, pv, mLen);              //       u = (u-v);  in alm_mont_inv()
101             cpLSR_BNU(pu, pu, mLen, 1);               //       u = u/2;  in alm_mont_inv()
102             /*ext +=*/ cpAdd_BNU(pt, pt, ps, mLen);   //       t = t+s;  in alm_mont_inv()
103             cpAdd_BNU(ps, ps, ps, mLen);              //       s = 2*s;  in alm_mont_inv()
106             cpSub_BNU(pv, pv, pu, mLen);              //       v = (v-u);  in alm_mont_inv()
107             cpLSR_BNU(pv, pv, mLen, 1);               //       v = v/2;  in alm_mont_inv()
108             cpAdd_BNU(ps, ps, pt, mLen);              //       s = s+t;  in alm_mont_inv()
109             ext += cpAdd_BNU(pt, pt, pt, mLen);       //       t = 2*t;  in alm_mont_inv()
116    if(1!=cpEqu_BNU_CHUNK(pu, mLen, 1)) {  in alm_mont_inv()
121       ext -= cpSub_BNU(pr, pt, pm, mLen);             // if(t>mod) r = t-mod;  in alm_mont_inv()
122       cpMaskMove_gs(pr, pt, mLen, cpIsNonZero(ext));  // else r = t;  in alm_mont_inv()
123       cpSub_BNU(pr, pm, pr, mLen);                    // return  r= (mod - r) and k  in alm_mont_inv()