主頁 > 後端開發 > 在多個執行緒中運行時浮點運算是否具有確定性?

在多個執行緒中運行時浮點運算是否具有確定性?

2022-05-16 07:37:39 後端開發

假設我有一個運行計算的函式,例如點積 - 我傳入一個A, B向量陣列和一個浮點陣列C,函式分配: C[i] = dot(A[i], B[i]);

如果我創建并啟動兩個將運行此函式的執行緒,并將相同的三個陣列傳遞給兩個執行緒,那么在什么情況下這種型別的操作(可能使用不同的非亂數學運算等)不能保證給出相同的結果(在同一臺機器上運行相同的應用程式而無需任何重新編譯)?我只對消費類 PC 的環境感興趣。

我知道浮點操作通常是確定性的,但我確實想知道是否可能會發生一些奇怪的事情,并且可能在一個執行緒上計算將使用一個中間 80 位暫存器,但在另一個執行緒中不會。

我假設幾乎可以保證相同的二進制代碼應該在兩個執行緒中運行(有什么方法不會發生這種情況嗎?由于某種原因,該函式被多次編譯,編譯器以某種方式確定它將在多個執行緒中運行,然后編譯出于某種原因,第二個執行緒再次出現?)。但我有點擔心 CPU 內核可能沒有相同的指令集,即使在消費級 PC 上也是如此。

附帶問題 - 在類似情況下 GPU 怎么樣?

//

我假設 x86_64、Windows、c dota.x * b.x a.y * b.y. 無法提供更多資訊 - 使用 Unity IL2CPP,不知道它如何編譯/使用哪些選項。

問題的動機:我正在撰寫一個修改網格的計算幾何程式 - 我將其稱為“幾何網格”。問題是“渲染網格”可能會發生某些幾何位置的多個頂點 - 例如平面著色需要它 - 您有多個具有不同法線的頂點。然而,實際的計算幾何程序僅使用空間位置的純幾何資料。

所以我看到了兩個選擇:

  1. 創建從渲染網格到幾何網格的映射(示例 - 將重復頂點映射到一個唯一頂點),在幾何網格上運行該程序,然后根據結果以某種方式修改渲染網格。
  2. 直接使用渲染網格。由于該程序對所有頂點進行計算,因此效率稍低,但從代碼角度來看要容易得多。但最重要的是,我有點擔心我可能會為實際上具有相同位置的兩個頂點獲得兩個不同的值,這是不應該發生的。僅使用位置,并且兩個這樣的頂點的位置相同。

uj5u.com熱心網友回復:

浮點 (FP) 操作不是關聯的(但它是可交換的)。因此,(x y) z可以給出與 不同的結果x (y z)例如,對于(1e-13 (1 - 1e-13)) == ((1e-13 1) - 1e-13)64 位 IEEE-754 浮點數為 false。C 標準對浮點數的限制不是很嚴格。然而,廣泛使用的IEEE-754標準是。它指定 32 位和 64 位數字運算的精度,包括舍入模式。x86-64 處理器符合 IEEE-754,主流編譯器(例如 GCC、Clang 和 MSVC)默認也符合 IEEE-754。默認情況下,ICC 不兼容,因為它假定 FP 操作是關聯的以提高性能。主流編譯器具有編譯標志來做出這樣的假設,以加快代碼速度。它通常與其他一些相結合,例如假設所有 FP 值都不是 NaN(例如-ffast-math)。此類標志違反了 IEEE-754 合規性,但它們通常用于 3D 或視頻游戲行業,以加快代碼速度。C 標準不需要 IEEE-754,但您可以使用std::numeric_limits<T>::is_iec559.

默認情況下,執行緒可以有不同的舍入模式但是,您可以使用此答案中提供的 C 代碼設定舍入模式另外,請注意,非正規數有時會在某些平臺上被禁用,因為它們的開銷非常高(有關更多資訊,請參閱內容)。

假設沒有破壞 IEEE-754 合規性,舍入模式相同并且執行緒以相同的順序執行操作,那么結果應該至少相同 1 ULP在實踐中,如果它們使用同一個主流編譯器編譯,結果應該是完全一樣的。

問題是使用多個執行緒通常會導致應用的 FP 操作的順序不確定,從而導致結果不確定。更具體地說,對 FP 變數的原子操作通常會導致這樣的問題,因為操作的順序經常在運行時發生變化。如果你想要確定性的結果,你需要使用靜態磁區,避免對 FP 變數的原子操作或更一般的可能導致不同排序的原子操作。同樣的事情也適用于鎖或任何同步機制。

GPU 也是如此。事實上,當開發人員使用原子 FP 操作(例如對值求和)時,此類問題非常常見。他們經常這樣做是因為實作快速歸約很復雜(盡管它更具確定性)并且原子操作在現代 GPU 上非常快(因為它們使用專用的高效單元)。

uj5u.com熱心網友回復:

根據浮點處理器非確定性的公認答案?, C 浮點不是非確定性的。相同的指令序列將給出相同的結果。

不過有幾點需要考慮:

首先,執行 FP 計算的特定 C 源代碼的行為(即結果)可能取決于編譯器和選擇的編譯器選項。例如,它可能取決于編譯器是選擇發出 64 位還是 80 位 FP 指令。但這是確定性的。

其次,類似的 C 源代碼可能會給出不同的結果;例如,由于某些 FP 指令的非關聯行為。這也是確定性的。

默認情況下,確定性不會受到多執行緒的影響。C 編譯器可能不知道代碼是否是多執行緒的。它絕對沒有理由發出不同的 FP 代碼。

誠然,FP 行為取決于所選的舍入模式,并且可以基于每個執行緒進行設定。然而,要做到這一點,某些東西(應用程式代碼)必須為不同的執行緒顯式設定不同的舍入模式。再一次,這是確定性的。(對于應用程式代碼來說,這是一件非常愚蠢的事情,IMO。)


PC 會為不同執行緒使用具有不同行為的不同 FP 硬體的想法對我來說似乎有些牽強。當然,一臺 PC 可以有(比如說)一個 Intel 芯片組和一個 ARM 芯片組,但是同一個 C 應用程式(可執行檔案)的不同執行緒同時在兩個芯片組上運行是不合理的。

對于 GPU 也是如此。事實上,考慮到您需要以與普通(或執行緒)C 完全不同的方式對 GPU 進行編程,我懷疑它們甚至可以共享相同的源代碼。


簡而言之,我認為您正在擔心一個您在現實中不太可能遇到的假設問題……鑒于硬體和 C 編譯器的當前狀態。

轉載請註明出處,本文鏈接:https://www.uj5u.com/houduan/475204.html

標籤:多线程 浮点 确定性的 非确定性的

上一篇:在rust中實作執行緒安全的函式物件

下一篇:在從.exe寫入輸出和讀取輸出之間拆分作業

標籤雲
其他(157675) Python(38076) JavaScript(25376) Java(17977) C(15215) 區塊鏈(8255) C#(7972) AI(7469) 爪哇(7425) MySQL(7132) html(6777) 基礎類(6313) sql(6102) 熊猫(6058) PHP(5869) 数组(5741) R(5409) Linux(5327) 反应(5209) 腳本語言(PerlPython)(5129) 非技術區(4971) Android(4554) 数据框(4311) css(4259) 节点.js(4032) C語言(3288) json(3245) 列表(3129) 扑(3119) C++語言(3117) 安卓(2998) 打字稿(2995) VBA(2789) Java相關(2746) 疑難問題(2699) 细绳(2522) 單片機工控(2479) iOS(2429) ASP.NET(2402) MongoDB(2323) 麻木的(2285) 正则表达式(2254) 字典(2211) 循环(2198) 迅速(2185) 擅长(2169) 镖(2155) 功能(1967) .NET技术(1958) Web開發(1951) python-3.x(1918) HtmlCss(1915) 弹簧靴(1913) C++(1909) xml(1889) PostgreSQL(1872) .NETCore(1853) 谷歌表格(1846) Unity3D(1843) for循环(1842)

熱門瀏覽
  • 【C++】Microsoft C++、C 和匯編程式檔案

    ......

    uj5u.com 2020-09-10 00:57:23 more
  • 例外宣告

    相比于斷言適用于排除邏輯上不可能存在的狀態,例外通常是用于邏輯上可能發生的錯誤。 例外宣告 Item 1:當函式不可能拋出例外或不能接受拋出例外時,使用noexcept 理由 如果不打算拋出例外的話,程式就會認為無法處理這種錯誤,并且應當盡早終止,如此可以有效地阻止例外的傳播與擴散。 示例 //不可 ......

    uj5u.com 2020-09-10 00:57:27 more
  • Codeforces 1400E Clear the Multiset(貪心 + 分治)

    鏈接:https://codeforces.com/problemset/problem/1400/E 來源:Codeforces 思路:給你一個陣列,現在你可以進行兩種操作,操作1:將一段沒有 0 的區間進行減一的操作,操作2:將 i 位置上的元素歸零。最終問:將這個陣列的全部元素歸零后操作的最少 ......

    uj5u.com 2020-09-10 00:57:30 more
  • UVA11610 【Reverse Prime】

    本人看到此題沒有翻譯,就附帶了一個自己的翻譯版本 思考 這一題,它的第一個要求是找出所有 $7$ 位反向質數及其質因數的個數。 我們應該需要質數篩篩選1~$10^{7}$的所有數,這里就不慢慢介紹了。但是,重讀題,我們突然發現反向質數都是 $7$ 位,而將它反過來后的數字卻是 $6$ 位數,這就說明 ......

    uj5u.com 2020-09-10 00:57:36 more
  • 統計區間素數數量

    1 #pragma GCC optimize(2) 2 #include <bits/stdc++.h> 3 using namespace std; 4 bool isprime[1000000010]; 5 vector<int> prime; 6 inline int getlist(int ......

    uj5u.com 2020-09-10 00:57:47 more
  • C/C++編程筆記:C++中的 const 變數詳解,教你正確認識const用法

    1、C中的const 1、區域const變數存放在堆疊區中,會分配記憶體(也就是說可以通過地址間接修改變數的值)。測驗代碼如下: 運行結果: 2、全域const變數存放在只讀資料段(不能通過地址修改,會發生寫入錯誤), 默認為外部聯編,可以給其他源檔案使用(需要用extern關鍵字修飾) 運行結果: ......

    uj5u.com 2020-09-10 00:58:04 more
  • 【C++犯錯記錄】VS2019 MFC添加資源不懂如何修改資源宏ID

    1. 首先在資源視圖中,添加資源 2. 點擊新添加的資源,復制自動生成的ID 3. 在解決方案資源管理器中找到Resource.h檔案,編輯,使用整個專案搜索和替換的方式快速替換 宏宣告 4. Ctrl+Shift+F 全域搜索,點擊查找全部,然后逐個替換 5. 為什么使用搜索替換而不使用屬性視窗直 ......

    uj5u.com 2020-09-10 00:59:11 more
  • 【C++犯錯記錄】VS2019 MFC不懂的批量添加資源

    1. 打開資源頭檔案Resource.h,在其中預先定義好宏 ID(不清楚其實ID值應該設定多少,可以先新建一個相同的資源項,再在這個資源的ID值的基礎上遞增即可) 2. 在資源視圖中選中專案資源,按F7編輯資源檔案,按 ID 型別 相對路徑的形式添加 資源。(別忘了先把檔案拷貝到專案中的res檔案 ......

    uj5u.com 2020-09-10 01:00:19 more
  • C/C++編程筆記:關于C++的參考型別,專供新手入門使用

    今天要講的是C++中我最喜歡的一個用法——參考,也叫別名。 參考就是給一個變數名取一個變數名,方便我們間接地使用這個變數。我們可以給一個變數創建N個參考,這N + 1個變數共享了同一塊記憶體區域。(參考型別的變數會占用記憶體空間,占用的記憶體空間的大小和指標型別的大小是相同的。雖然參考是一個物件的別名,但 ......

    uj5u.com 2020-09-10 01:00:22 more
  • 【C/C++編程筆記】從頭開始學習C ++:初學者完整指南

    眾所周知,C ++的學習曲線陡峭,但是花時間學習這種語言將為您的職業帶來奇跡,并使您與其他開發人員區分開。您會更輕松地學習新語言,形成真正的解決問題的技能,并在編程的基礎上打下堅實的基礎。 C ++將幫助您養成良好的編程習慣(即清晰一致的編碼風格,在撰寫代碼時注釋代碼,并限制類內部的可見性),并且由 ......

    uj5u.com 2020-09-10 01:00:41 more
最新发布
  • Rust中的智能指標:Box<T> Rc<T> Arc<T> Cell<T> RefCell<T> Weak

    Rust中的智能指標是什么 智能指標(smart pointers)是一類資料結構,是擁有資料所有權和額外功能的指標。是指標的進一步發展 指標(pointer)是一個包含記憶體地址的變數的通用概念。這個地址參考,或 ” 指向”(points at)一些其 他資料 。參考以 & 符號為標志并借用了他們所 ......

    uj5u.com 2023-04-20 07:24:10 more
  • Java的值傳遞和參考傳遞

    值傳遞不會改變本身,參考傳遞(如果傳遞的值需要實體化到堆里)如果發生修改了會改變本身。 1.基本資料型別都是值傳遞 package com.example.basic; public class Test { public static void main(String[] args) { int ......

    uj5u.com 2023-04-20 07:24:04 more
  • [2]SpinalHDL教程——Scala簡單入門

    第一個 Scala 程式 shell里面輸入 $ scala scala> 1 + 1 res0: Int = 2 scala> println("Hello World!") Hello World! 檔案形式 object HelloWorld { /* 這是我的第一個 Scala 程式 * 以 ......

    uj5u.com 2023-04-20 07:23:58 more
  • 理解函式指標和回呼函式

    理解 函式指標 指向函式的指標。比如: 理解函式指標的偽代碼 void (*p)(int type, char *data); // 定義一個函式指標p void func(int type, char *data); // 宣告一個函式func p = func; // 將指標p指向函式func ......

    uj5u.com 2023-04-20 07:23:52 more
  • Django筆記二十五之資料庫函式之日期函式

    本文首發于公眾號:Hunter后端 原文鏈接:Django筆記二十五之資料庫函式之日期函式 日期函式主要介紹兩個大類,Extract() 和 Trunc() Extract() 函式作用是提取日期,比如我們可以提取一個日期欄位的年份,月份,日等資料 Trunc() 的作用則是截取,比如 2022-0 ......

    uj5u.com 2023-04-20 07:23:45 more
  • 一天吃透JVM面試八股文

    什么是JVM? JVM,全稱Java Virtual Machine(Java虛擬機),是通過在實際的計算機上仿真模擬各種計算機功能來實作的。由一套位元組碼指令集、一組暫存器、一個堆疊、一個垃圾回收堆和一個存盤方法域等組成。JVM屏蔽了與作業系統平臺相關的資訊,使得Java程式只需要生成在Java虛擬機 ......

    uj5u.com 2023-04-20 07:23:31 more
  • 使用Java接入小程式訂閱訊息!

    更新完微信服務號的模板訊息之后,我又趕緊把微信小程式的訂閱訊息給實作了!之前我一直以為微信小程式也是要企業才能申請,沒想到小程式個人就能申請。 訊息推送平臺🔥推送下發【郵件】【短信】【微信服務號】【微信小程式】【企業微信】【釘釘】等訊息型別。 https://gitee.com/zhongfuch ......

    uj5u.com 2023-04-20 07:22:59 more
  • java -- 緩沖流、轉換流、序列化流

    緩沖流 緩沖流, 也叫高效流, 按照資料型別分類: 位元組緩沖流:BufferedInputStream,BufferedOutputStream 字符緩沖流:BufferedReader,BufferedWriter 緩沖流的基本原理,是在創建流物件時,會創建一個內置的默認大小的緩沖區陣列,通過緩沖 ......

    uj5u.com 2023-04-20 07:22:49 more
  • Java-SpringBoot-Range請求頭設定實作視頻分段傳輸

    老實說,人太懶了,現在基本都不喜歡寫筆記了,但是網上有關Range請求頭的文章都太水了 下面是抄的一段StackOverflow的代碼...自己大修改過的,寫的注釋挺全的,應該直接看得懂,就不解釋了 寫的不好...只是希望能給視頻網站開發的新手一點點幫助吧. 業務場景:視頻分段傳輸、視頻多段傳輸(理 ......

    uj5u.com 2023-04-20 07:22:42 more
  • Windows 10開發教程_編程入門自學教程_菜鳥教程-免費教程分享

    教程簡介 Windows 10開發入門教程 - 從簡單的步驟了解Windows 10開發,從基本到高級概念,包括簡介,UWP,第一個應用程式,商店,XAML控制元件,資料系結,XAML性能,自適應設計,自適應UI,自適應代碼,檔案管理,SQLite資料庫,應用程式到應用程式通信,應用程式本地化,應用程式 ......

    uj5u.com 2023-04-20 07:22:35 more