cotexM3开方函数SQRT

a8200845

cotexM3开方函数SQRT [复制链接]

谁有cotex M3开方函数，是32位或64位整数开方的函数，多谢提供者，本人非常需要
keil c的函数库中只有double的开方函数，太慢了，我的50多个模拟量需要开方，光开方在72Mhz下都需要10多ms，这是不能接受的

iaw2005

                                  

ahoo

在线等待大家的好消息呀，急急的

哪位位大侠有呀

Abrolex

                                  

qianrj

                                  

nannan50972

                                   开方函数SQRT，记得在其它版面有人讨论过，速度号称很快

hailong201

                                  

huleisheng

                                   同意楼上，手工算开方是万能的，不管用哪个CPU（MCU，MPU，DSP）

977845

多谢大家的意见，手工开放是用十进制的，计算机怎么要实现这个手工开方，还不如用double库了
我要的是实质性的建议

luliku

                                 7楼这个建议就是一个实质性的建议，估计楼主是想要现成的程序吧。

searchaben

                                 http://www.power-bbs.com/zixun/jishuwenzhang2/20070903/18935.html

varrior

                                  

jjww

我已经找到了，手算的方法，我想使用TI的函数原型，可是M3没有数符号的指令
需要一位一位的移动太费劲了，所以我就使用了手算，循环16次，虽然比T1的32的条无跳转指令慢了一下，但还是可以接受的

多谢大家的指点

yanghui_45

                                  

lxxo

                                  

eugenew

                                   牛顿迭代，确实很快，对32位数据进行开方，大概只需3个循环；对64位数据进行开方，估计6个循环就够了

villaining

ARM优化那本书上有的。
俺老师用过

feixue112

                                 下次再来问问题，别人可能就没有这么热心了。

air_cc

                                  

andrewlzj

unsigned int Real_Sqrt_16(unsigned long int M)
{
    unsigned int N, i; 
    unsigned long tmp, ttp;   // 结果、循环计数 
    if (M == 0)               // 被开方数，开方结果也为0      
       return 0;     
    N = 0;     
    tmp = (M >> 30);          // 获取最高位：B[m-1]   
    M <<= 2;     
    if (tmp > 1)              // 最高位为1    
    {   
        N ++;                 // 结果当前位为1，否则为默认的0     
        tmp -= N;     
    } 
    for (i=15; i>0; i--)      // 求剩余的15位   
     {         
        N <<= 1;              // 左移一位 
        tmp <<= 2;
          tmp += (M >> 30);     // 假设        
            ttp = N;         
        ttp = (ttp<<1)+1;         
        M <<= 2;         
        if (tmp >= ttp)       // 假设成立        
         {             
        tmp -= ttp;             
        N ++;         
        }     
    }     
    
    return N; 
}