m_apm_sqrt
extern void m_apm_sqrt(M_APM, int, M_APM);
m_apm_sqrt(r, places, tmp2);
m_apm_sqrt(b0, dplaces, MM_0_5); /* sqrt(0.5) */
m_apm_sqrt(b1, dplaces, tmp1);
m_apm_sqrt(tmpB0, dplaces, tmp2);
m_apm_sqrt(tmp9, (places + 6), tmp8);
m_apm_sqrt(tmp1, (places + 6), tmp2);
m_apm_sqrt(tmp1, (places + 6), tmp2);