0
  • 聊天消息
  • 系統(tǒng)消息
  • 評(píng)論與回復(fù)
登錄后你可以
  • 下載海量資料
  • 學(xué)習(xí)在線課程
  • 觀看技術(shù)視頻
  • 寫文章/發(fā)帖/加入社區(qū)
會(huì)員中心
創(chuàng)作中心

完善資料讓更多小伙伴認(rèn)識(shí)你,還能領(lǐng)取20積分哦,立即完善>

3天內(nèi)不再提示

浮點(diǎn)庫應(yīng)用,你也有困惑嗎?

恩智浦MCU加油站 ? 來源:未知 ? 2023-12-21 08:55 ? 次閱讀

今天小編想要給大家?guī)淼氖亲罱谡{(diào)試一個(gè)項(xiàng)目時(shí)候發(fā)現(xiàn)的怪事,同樣的函數(shù)庫在使用不同的IDE時(shí),得到的運(yùn)行結(jié)果竟然是不一致的。相信眼尖的讀者朋友已經(jīng)從標(biāo)題中猜出一二了,正是因?yàn)閹熘幸蕾嚵烁↑c(diǎn)計(jì)算庫所導(dǎo)致的問題。那么就請(qǐng)和小編一起,探究下詳細(xì)的來龍去脈吧!

項(xiàng)目背景

首先再詳細(xì)描述以下項(xiàng)目背景:我們使用了一個(gè)由GCC工具鏈構(gòu)建的函數(shù)庫,編譯器選項(xiàng)為-mfloat-abi=hard,即在編譯時(shí),使能了硬件浮點(diǎn)單元指令加速。但在我們將該庫與Keil項(xiàng)目鏈接后,發(fā)現(xiàn)結(jié)果不正確。當(dāng)然,我們最先懷疑的當(dāng)然是庫編譯的有問題啦。但是奇怪的事情發(fā)生了,我當(dāng)我們使用GCC工具鏈鏈接編譯相同的工程,并運(yùn)行后,得到了正確的結(jié)果。那么可以基本確定,庫應(yīng)該是沒有問題的。那么問題出在哪呢?

問題分析

由于小編所拿到的庫并不包含調(diào)試信息,只能通過一些技術(shù)手段進(jìn)行破解。通過反匯編庫代碼,我們發(fā)現(xiàn)庫依賴了一些浮點(diǎn)計(jì)算庫的C函數(shù),如sqrtf、expf等,為了簡便,讓我們稱之為xxf函數(shù),由于GCC并沒有提供其具體實(shí)現(xiàn),因此需要由庫的使用者鏈接這些函數(shù),而經(jīng)過進(jìn)一步的debug,我們發(fā)現(xiàn)這些函數(shù)導(dǎo)致了錯(cuò)誤的結(jié)果,換句話說,這些函數(shù)的返回值是錯(cuò)的。

問題調(diào)試

那就有意思了,為了簡單起見,讓我們先編寫一段簡單的測試代碼來復(fù)現(xiàn)這個(gè)問題,代碼很簡單,我們直接定義一個(gè)開平方根的函數(shù):

#include "math.h"

float calc_sqrt(float a){

return sqrtf(a);

}

使用gcc工具鏈對(duì)其進(jìn)行編譯:

arm-none-eabi-gcc.exe -mcpu=cortex-m7 -mthumb -mfpu=fpv5-sp-d16 -mfloat-abi=hard test.c -fshort-wchar -c -o test.o

隨后,隨便找到一個(gè)Keil的測試工程,我們這里選擇一個(gè)Hello_World示例工程,將編譯出來的.o文件添加到工程中:

wKgZomWDjhGANp5cAABcyFsEcJc683.png

并在主程序中添加調(diào)用代碼:

float calc_sqrt(float a)

volatile float a = calc_sqrt(4.0f);

編譯鏈接下載程序,并讓程序停在函數(shù)調(diào)用處:

wKgZomWDjhKAUXnTAADNmNIH8Rw407.png

單步進(jìn)入calc_sqrt函數(shù)內(nèi)部,到這里,我們可以發(fā)現(xiàn)對(duì)這些函數(shù)的調(diào)用順序是正確的。通過將參數(shù)傳遞給S0(對(duì)于float)如下所示,S0中寸的就是待計(jì)算的數(shù)據(jù)4.0f:

wKgZomWDjhKAcUdQAAFHqs3OJMk254.png

wKgZomWDjhKAbZ7dAAA743Nmdr0045.png

看起來好像沒有問題,再進(jìn)一步現(xiàn)在讓我們檢查由Keil鏈接的sqrtf函數(shù)的匯編實(shí)現(xiàn):

wKgZomWDjhKAF50tAADszXvbpi8487.png

相信大家發(fā)現(xiàn)了奇怪的事情了,鏈接的sqrtf將S0中的值傳遞給S0,而此時(shí)R0的值其實(shí)為0,但正如之前所說,浮點(diǎn)值已經(jīng)由庫的代碼傳遞給了S0。因此,由于S0中實(shí)際要計(jì)算的值被臨時(shí)替換了,就導(dǎo)致了一個(gè)錯(cuò)誤的結(jié)果。

wKgZomWDjhKAFognAABSVmp7PxU258.png

這里要強(qiáng)調(diào)一下,如果在Keil中直接調(diào)用sqrtf時(shí)候,或是使用keil編譯器所編譯出來的.o文件,Keil運(yùn)行時(shí)庫會(huì)使用“__hardfp_sqrtf”作為sqrtf的混淆名稱:

wKgZomWDjhKAZvr5AACqgxPX2Ic445.png

而因?yàn)槲覀兯褂脦靵碜訥CC工具鏈,因此Keil并不會(huì)對(duì)其進(jìn)行替換,而是會(huì)將C庫中叫做sqrtf的函數(shù)直接鏈接進(jìn)去,而這個(gè)函數(shù)的默認(rèn)實(shí)現(xiàn),是使用R0作為參數(shù)傳遞的寄存器。這也就導(dǎo)致,實(shí)際要被計(jì)算的數(shù)丟失,最終導(dǎo)致結(jié)果計(jì)算錯(cuò)誤。

那么怎么解決這個(gè)問題,讓keil不去鏈接這個(gè)奇怪的sqrtf呢,這就要用到Keil的一個(gè)小技巧了:

float $Sub$$sqrtf(float a){

return __builtin_sqrtf(a);

}

這樣一來,調(diào)用sqrtf函數(shù)的地方,就會(huì)編程對(duì)$Sub$$sqrtf的調(diào)用:

wKgZomWDjhKAMFbbAAE1zkGftN4411.png

wKgZomWDjhOACcEzAAIb8mxu1DA400.png

而結(jié)果也變?yōu)檎_的結(jié)果2.0了:

wKgZomWDjhOAJ6olAABTCNVByDY135.png

當(dāng)然,大家可能會(huì)想啊,那我開了gcc優(yōu)化之后,是不是就可以自動(dòng)展開sqrtf了呢?讓我們來看看:

arm-none-eabi-gcc.exe -mcpu=cortex-m7 -mthumb -mfpu=fpv5-sp-d16 -mfloat-abi=hard test.c -fshort-wchar -c -o3 -o test.o

代碼并沒有變化:

wKgZomWDjhOAE65GAAEt3Qti-14440.png

結(jié)論

小編想用這個(gè)例子和大家說明下,在涉及到跨工具鏈開發(fā)時(shí),一定要注意浮點(diǎn)庫的使用或依賴問題,由于不同編譯器對(duì)于浮點(diǎn)運(yùn)算的實(shí)現(xiàn)可能有些許不同,會(huì)導(dǎo)致意想不到的奇怪問題出現(xiàn)。最好的方案,還是根據(jù)不同的工具鏈都構(gòu)建一個(gè)專屬的庫來使用。

END

更多恩智浦AI-IoT市場和產(chǎn)品信息,邀您同時(shí)關(guān)注“NXP客?!蔽⑿殴娞?hào)

wKgZomWDjhOACI5YAABCdkRE230322.jpg ? ? ?

NXP客棧


恩智浦致力于打造安全的連接和基礎(chǔ)設(shè)施解決方案,為智慧生活保駕護(hù)航。

長按二維碼,關(guān)注我們

恩智浦MCU加油站


這是由恩智浦官方運(yùn)營的公眾號(hào),著重為您推薦恩智浦MCU的產(chǎn)品信息、開發(fā)技巧、教程文檔、培訓(xùn)課程等內(nèi)容。

wKgZomWDjhOAHK_FAAATNlPH08Y075.jpg ?

長按二維碼,關(guān)注我們


原文標(biāo)題:浮點(diǎn)庫應(yīng)用,你也有困惑嗎?

文章出處:【微信公眾號(hào):恩智浦MCU加油站】歡迎添加關(guān)注!文章轉(zhuǎn)載請(qǐng)注明出處。


聲明:本文內(nèi)容及配圖由入駐作者撰寫或者入駐合作網(wǎng)站授權(quán)轉(zhuǎn)載。文章觀點(diǎn)僅代表作者本人,不代表電子發(fā)燒友網(wǎng)立場。文章及其配圖僅供工程師學(xué)習(xí)之用,如有內(nèi)容侵權(quán)或者其他違規(guī)問題,請(qǐng)聯(lián)系本站處理。 舉報(bào)投訴
  • mcu
    mcu
    +關(guān)注

    關(guān)注

    146

    文章

    17171

    瀏覽量

    351486
  • 恩智浦
    +關(guān)注

    關(guān)注

    14

    文章

    5863

    瀏覽量

    107595

原文標(biāo)題:浮點(diǎn)庫應(yīng)用,你也有困惑嗎?

文章出處:【微信號(hào):NXP_SMART_HARDWARE,微信公眾號(hào):恩智浦MCU加油站】歡迎添加關(guān)注!文章轉(zhuǎn)載請(qǐng)注明出處。

收藏 人收藏

    評(píng)論

    相關(guān)推薦

    【RA-Eco-RA4E2-64PIN-V1.0開發(fā)板試用】RA4E2的DSP浮點(diǎn)性能的軟件浮點(diǎn)測試和硬件浮點(diǎn)測試對(duì)比

    , atan,等等基本操作。 當(dāng)然為了測試出硬件浮點(diǎn)運(yùn)算性能,這里很有必要在測試一下軟件浮點(diǎn)運(yùn)算能力,方便做個(gè)測試對(duì)比,怎么做軟件浮點(diǎn)運(yùn)算呢? 這個(gè)可以使用開源Qfplib,這是國外
    發(fā)表于 12-30 17:55

    FPGA中的浮點(diǎn)四則運(yùn)算是什么

    由于定點(diǎn)的四則運(yùn)算比較簡單,如加減法只要注意符號(hào)擴(kuò)展,小數(shù)點(diǎn)對(duì)齊等問題即可。在本文中,運(yùn)用在前一節(jié)中描述的自定義浮點(diǎn)格式FPGA中數(shù)的表示方法(下),完成浮點(diǎn)四則運(yùn)算的實(shí)現(xiàn)過程 1.自定義浮點(diǎn)格式加
    的頭像 發(fā)表于 11-16 12:51 ?285次閱讀
    FPGA中的<b class='flag-5'>浮點(diǎn)</b>四則運(yùn)算是什么

    FPGA中浮點(diǎn)四則運(yùn)算的實(shí)現(xiàn)過程

    由于定點(diǎn)的四則運(yùn)算比較簡單,如加減法只要注意符號(hào)擴(kuò)展,小數(shù)點(diǎn)對(duì)齊等問題即可。在本文中,運(yùn)用在前一節(jié)中描述的自定義浮點(diǎn)格式FPGA中數(shù)的表示方法(下),完成浮點(diǎn)四則運(yùn)算的實(shí)現(xiàn)過程 1.自定義浮點(diǎn)格式加
    的頭像 發(fā)表于 11-16 11:19 ?468次閱讀
    FPGA中<b class='flag-5'>浮點(diǎn)</b>四則運(yùn)算的實(shí)現(xiàn)過程

    【RA-Eco-RA2E1-48PIN-V1.0開發(fā)板試用】在M23內(nèi)核上使用qfplib浮點(diǎn)運(yùn)算進(jìn)行浮點(diǎn)運(yùn)算

    運(yùn)算。難道由于硬件不支持FPU, 常用的M0/M0+/M23/M3內(nèi)核就無緣浮點(diǎn)運(yùn)算了呢?答案是顯然不是的。 我們可以移植使用開源的qfplib來實(shí)現(xiàn)浮點(diǎn)運(yùn)算,這樣在一些不具備FPU的低性能MCU上
    發(fā)表于 11-05 22:07

    【RV32F與RV32D】之單精度和雙精度浮點(diǎn)數(shù)

    和所有其它31個(gè)f寄存器一樣,是一個(gè)可變寄存器。 IEEE 754-2008標(biāo)準(zhǔn)提供了幾種浮點(diǎn)運(yùn)算舍入的方法,這有助于確定誤差范圍和編寫數(shù)值。最準(zhǔn)確且最常見的舍入模式是舍入到最近的偶數(shù)(RNE
    發(fā)表于 08-29 12:24

    TMS320C6742定點(diǎn)和浮點(diǎn)DSP數(shù)據(jù)表

    電子發(fā)燒友網(wǎng)站提供《TMS320C6742定點(diǎn)和浮點(diǎn)DSP數(shù)據(jù)表.pdf》資料免費(fèi)下載
    發(fā)表于 08-07 09:57 ?0次下載
    TMS320C6742定點(diǎn)和<b class='flag-5'>浮點(diǎn)</b>DSP數(shù)據(jù)表

    TMS320C6746定點(diǎn)和浮點(diǎn)DSP數(shù)據(jù)表

    電子發(fā)燒友網(wǎng)站提供《TMS320C6746定點(diǎn)和浮點(diǎn)DSP數(shù)據(jù)表.pdf》資料免費(fèi)下載
    發(fā)表于 08-07 09:55 ?0次下載
    TMS320C6746定點(diǎn)和<b class='flag-5'>浮點(diǎn)</b>DSP數(shù)據(jù)表

    ESP8266_RTOS_SDK如何打印浮點(diǎn)數(shù)?

    嘗試使用 printf 打印浮點(diǎn)數(shù),例如, 浮點(diǎn)數(shù) d = 1.0; printf("d = %f", d); 但它打印出來: d = 此外,sprintf
    發(fā)表于 07-09 06:56

    STM32CubeIDE printf浮點(diǎn)數(shù),浮點(diǎn)數(shù)丟失數(shù)值的原因?

    1、芯片:STM32F103C8T62、環(huán)境:STM32CubeIDE3、問題:printf浮點(diǎn)數(shù),浮點(diǎn)數(shù)數(shù)值丟失如 float A=25.6666;float B=24.7777;實(shí)際打印
    發(fā)表于 06-03 07:07

    優(yōu)秀的Verilog/FPGA開源項(xiàng)目-浮點(diǎn)運(yùn)算器(FPU)介紹

    浮點(diǎn)運(yùn)算器(英文:floating point unit,簡稱FPU)是計(jì)算機(jī)系統(tǒng)的一部分,它是專門用來進(jìn)行浮點(diǎn)數(shù)運(yùn)算的(CPU中也叫ALU)。
    的頭像 發(fā)表于 04-26 11:27 ?3503次閱讀
    優(yōu)秀的Verilog/FPGA開源項(xiàng)目-<b class='flag-5'>浮點(diǎn)</b>運(yùn)算器(FPU)介紹

    STM32cubemx的和KEIL的有什么不同?

    STM32cubemx軟件有對(duì)應(yīng)各個(gè)系列芯片的. 使用cubemx,生成代碼框架,比如生成MDK的代碼框架。但MDK也有自己的PACK。那么MDK的PACK和CUBEMX的有什么區(qū)別,兩者時(shí)什么關(guān)系呢?
    發(fā)表于 04-11 06:02

    verilog語音實(shí)現(xiàn)浮點(diǎn)運(yùn)算

    Verilog可以通過使用IEEE標(biāo)準(zhǔn)的浮點(diǎn)數(shù)表示來實(shí)現(xiàn)浮點(diǎn)運(yùn)算。下面是一個(gè)基本的Verilog模塊示例,展示了如何進(jìn)行加法、乘法和除法等常見的浮點(diǎn)運(yùn)算操作: module
    發(fā)表于 03-25 21:49

    一文帶你秒懂IEEE 754浮點(diǎn)數(shù)

    一、簡介1、常見的浮點(diǎn)數(shù)表示方式是IEEE754標(biāo)準(zhǔn),它規(guī)定了浮點(diǎn)數(shù)的存儲(chǔ)格式和運(yùn)算規(guī)則,這個(gè)標(biāo)準(zhǔn)定義了兩種浮點(diǎn)數(shù)表示:單精度和雙精度。2、任何一個(gè)浮點(diǎn)數(shù)的二進(jìn)制數(shù)可以寫為:NUM=(
    的頭像 發(fā)表于 03-18 08:09 ?8852次閱讀
    一文帶你秒懂IEEE 754<b class='flag-5'>浮點(diǎn)</b>數(shù)

    使用的H7系列MX編寫代碼后生成靜態(tài)文件,如何更改靜態(tài)工程也使用硬件雙精度浮點(diǎn)數(shù)?

    使用的H7系列MX編寫代碼后生成靜態(tài)文件, 在調(diào)用該靜態(tài)文件并編譯報(bào)錯(cuò),如下: 然后我查看了靜態(tài)工程的編譯指令: 而APP程序的編譯指令: 想請(qǐng)教一下,該如何更改靜態(tài)工程
    發(fā)表于 03-12 06:15

    毫米波雷達(dá)半精度浮點(diǎn)存儲(chǔ)格式分析

    本文介紹了TC3xx單片機(jī)雷達(dá)信號(hào)處理單元SPU支持的半精度浮點(diǎn)格式,將其和32bit整型數(shù)格式進(jìn)行比較,分析了兩者的動(dòng)態(tài)范圍及實(shí)際處理誤差,發(fā)現(xiàn)半精度浮點(diǎn)格式是“性價(jià)比”較高的存儲(chǔ)方式。
    的頭像 發(fā)表于 02-20 08:26 ?383次閱讀
    毫米波雷達(dá)半精度<b class='flag-5'>浮點(diǎn)</b>存儲(chǔ)格式分析