cblas_sscal について

Question

前提

open blasライブラリについての質問です。
単精度浮動小数点数型に対するベクタのスカラ倍をする演算ライブラリで，
cblas_sscalについてですが，計算結果がおかしいので，使い方に問題があるのかを教えてほしく質問します。

openblas のデストレイビューションは複数あると思いますが私が使っているのは以下に記載したバージョンです。

発生している問題・エラーメッセージ

コードを実行すると，以下のように32番目から演算をしていないように思います。
nを63 以下にすると正確に演算するようです。

0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045
0.05 0.055 0.06 0.065 0.07 0.075 0.08 0.085 0.09 0.095 
0.1 0.105 0.11 0.115 0.12 0.125 0.13 0.135 0.14 0.145 
0.15 0.155 32 33 34 35 36 37 38 39 
40 41 42 43 44 45 46 47 48 49 
50 51 52 53 54 55 56 57 58 59 
60 61 62 63 

該当のソースコード

#include<numeric>
#include<vector>

int main(){

    int n = 64;
    std::vector<float> v(n);
    std::iota(v.begin(), v.end(), 0);
    cblas_sscal(n, 0.005f, v.data(), 1 );
       
    for(auto& i:v)
        std::cout << i << ' ';
    
    std::cout <<'\n';
}

補足情報（FW/ツールのバージョンなど）

https://www.openblas.net/
openblas 0.3.21
clang
C++17

Answer 1

ここの https://github.com/xianyi/OpenBLAS/releases OpenBLAS-0.3.21-x64.zipをつかって、Windows / VC2019でコードをコンパイルして実行しましたが、再現しませんでした。

0 0.005 0.01 0.015 0.02 0.025 0.03 0.035 0.04 0.045 0.05 0.055 0.06 0.065 0.07 0.075 0.08 0.085 0.09 0.095 0.1 0.105 0.11 0.115 0.12 0.125 0.13 0.135 0.14 0.145 0.15 0.155 0.16 0.165 0.17 0.175 0.18 0.185 0.19 0.195 0.2 0.205 0.21 0.215 0.22 0.225 0.23 0.235 0.24 0.245 0.25 0.255 0.26 0.265 0.27 0.275 0.28 0.285 0.29 0.295 0.3 0.305 0.31 0.315

Comments

[jbe00214]

再現の検証ありがとうございました。私の使い方は正しいのであれば，PC内のライブラリの問題かもしれませんので改めて点検してみます。どうもありがとうございました。

Answer 2

原因がわかりました。
私の現の環境は以下の通りです。

xcodeの 'Build settings'において，'Linking' オプションの'Dead Code Stripping' オプションをYesに設定すると，不正な結果となり，Noに設定すると，正しい結果となりました。xcode 13までのバージョンでは発生していなかった問題ですが，参考までにお知らせします。

なお，Xcode 13（たぶん）までは，Dead Code Stripping オプションはデフォルトでYesだったのですが，Xcode14（たぶん）からはNoに変更されています。

環境
macos 13.0.1
Xcode 14.1
clang
c++17
openblas 0.3.21
staticlibraly