sin関数の計算結果がGR-SAKURAと違う？

ド素人ですｗ。以前にGR-SAKURAのフォーラムで「8×8×8のLED-CUBEを作ったけど、うまく動かない」で助けて頂いた者です。

https://japan.renesasrulz.com/gr_user_forum_japanese/f/gr-sakura/3717/thread

同じLED-CUBEをとある事情で、SAKURA－＞秋月で買ったRX621ボードへの移植をしてます。（せっかくかわいいSAKURAボードは、他の用事に使おうかとｗ）

e2studioで作ってます。なんとかほぼ動くようになったのですが、sinを使ってる所でSAKURAとLEDの表示が違います…。

SAKURAの方のプログラムの該当部分は、

1 : void ripples(int times, int repeat)
2 :  {
3 :    int i,j,k,z;
4 :    float distance,ztemp;
5 :    clearData();

6 :    for(k=0;k<=repeat;++k) {
7 :    for(i=0;i<=3;++i) {
8 :    for(j=0;j<=3;++j) {
9 :    distance=sqrt((3.5-(float)i)*(3.5-(float)i)+(3.5-(float)j)*(3.5-(float)j));
10:    ztemp=3.5+sin(distance/1.3+(float)k/3.14159)*3.5;
11:    z=(int)(ztemp+0.5);
12:    bitSet(dispData[z][i],j);
13:    bitSet(dispData[z][7-i],j);
14:    bitSet(dispData[z][i],7-j);
15:    bitSet(dispData[z][7-i],7-j);
16:    }
17:    }
18:    LEDDisp(times);
19:    clearData();
20:    }
21:  }

LED表示用配列の辺は見ないで頂いて（移植に伴いbitsetー＞ビットシフトに変えてます）、9～11行目の式で、がじぇるねのwebコンパイラでビルドしたSAKURAと、e2studioでビルドしたRX621ボードとで、計算結果が異なる原因が何かあるでしょうか…？

sqrtとsinをsqrtfとsinfに変えてみたりとか、型指定を外してみたりとかしても、変わらない様です。

（printfで計算結果を確認しようとしたのですがうまくいかないので、ほんとに計算結果が違うかどうかは確認できていません…）

どなたかご助言頂けると助かります。

Parents

fujita nozomu over 3 years ago

sin や sqrt 等数学関数は math.h をインクルードしないと正常動作しませんがされてますか
Cancel
Up 0 Down

Reply

Cancel
kohei over 3 years ago in reply to fujita nozomu

もちろん、エラーが出ずにビルド・書き込み・実物実行してるので、インクルードしてます。
（今見たら…、sin・sqrtとsinf・sqrtfと、どちらを使うべきなのか分からないので、math.hとmathf.hと両方インクルードしてたつもりだったのですが、LibraryGeneratorの「構成」の画面（先ほどのrepに張り付けてる）でmathf.hにチェック入れてるけど、インクルードはmath.hしか記述してませんでした。）
Cancel
Up 0 Down

Reply

Cancel
fujita nozomu over 3 years ago in reply to kohei

> もちろん、エラーが出ずにビルド・書き込み・実物実行してるので、

関数プロトタイプが含まれるヘッダファイルのインクルードを忘れてもエラーとはならずに警告が出るだけだと思いますが一つの警告も出ていないでしょうか?
Cancel
Up 0 Down

Reply

Cancel
kohei over 3 years ago in reply to fujita nozomu

インクルードはしてるつもりなのですが…。警告は「W0561017:The evaluation period has expired」というのが１つ出てるようですが…。
Cancel
Up 0 Down

Reply

Cancel

Reply

kohei over 3 years ago in reply to fujita nozomu

インクルードはしてるつもりなのですが…。警告は「W0561017:The evaluation period has expired」というのが１つ出てるようですが…。
Cancel
Up 0 Down

Reply

Cancel

Children

fujita nozomu over 3 years ago in reply to kohei

> 警告は「W0561017:The evaluation period has expired」というのが１つ

「お試し期間終わったから金払え」ということですがとりあえず無視して良いです。
Cancel
Up 0 Down

Reply

Cancel
kohei over 3 years ago in reply to fujita nozomu

あ、そのヤツか～ｗｗ。プログラムのサイズ制限があるんだったかしら？（辞書引いても分からなかったｗ）
ありがとです～。（「変数やユーザー定義関数の評価をせー」とかいう意味かと思ったｗ。）
Cancel
Up 0 Down

Reply

Cancel

Okamiya Yuuki over 3 years ago in reply to kohei

がじぇるね岡宮です。

＃すみません、コードの中で、for文のkの演算が足りてないですね。つねに0で計算しているため、不十分です。

CCRXと、Webコンパイラで使っているGCCで違いが出るか、GR-SAKURAで試してみましたが、結果としては一緒でした。

for文2個で4回ずつ計16回ループしているため、かなり限定的な演算の比較なので、どれくらいsinの結果が一致してるかは未検証です。できれば、どういった差分が出るのか教えていただけるとありがたいです。

■Renesas CCRX V3.02.00

（新規プロジェクト作成で、FITモジュールのbspのみ適用、LibGenでmath.hを追加）

コード

#include <stdio.h>
#include <math.h>
void main(void)
{
 int i,j,k,z;
 float distance,ztemp;
 for(i=0;i<=3;++i) {
 for(j=0;j<=3;++j) {
 distance=sqrt((3.5-(float)i)*(3.5-(float)i)+(3.5-(float)j)*(3.5-(float)j));
 ztemp=3.5+sin(distance/1.3+(float)k/3.14159)*3.5;
 z=(int)(ztemp+0.5);
 printf("%d",z);printf(",");printf("%d",i);printf(",");printf("%d",j);
 printf("\r\n");
 printf("%d",z);printf(",");printf("%d",7-i);printf(",");printf("%d",j);
 printf("\r\n");
 printf("%d",z);printf(",");printf("%d",i);printf(",");printf("%d",7-j);
 printf("\r\n");
 printf("%d",z);printf(",");printf("%d",7-i);printf(",");printf("%d",7-j);
 printf("\r\n");
 printf("\r\n");
 }
 }
 while(1) {

 // TODO: add application code here
 }

}

結果

1,0,0
1,7,0
1,0,7
1,7,7

3,0,1
3,7,1
3,0,6
3,7,6

4,0,2
4,7,2
4,0,5
4,7,5

5,0,3
5,7,3
5,0,4
5,7,4

3,1,0
3,6,0
3,1,7
3,6,7

5,1,1
5,6,1
5,1,6
5,6,6

6,1,2
6,6,2
6,1,5
6,6,5

7,1,3
7,6,3
7,1,4
7,6,4

4,2,0
4,5,0
4,2,7
4,5,7

6,2,1
6,5,1
6,2,6
6,5,6

7,2,2
7,5,2
7,2,5
7,5,5

7,2,3
7,5,3
7,2,4
7,5,4

5,3,0
5,4,0
5,3,7
5,4,7

7,3,1
7,4,1
7,3,6
7,4,6

7,3,2
7,4,2
7,3,5
7,4,5

5,3,3
5,4,3
5,3,4
5,4,4

■Webコンパイラ（KPIT GNURX-ELF Toolchain 14.03.. ちょっと古くて恐縮ですが、Mac版IDE for GRと合わせるためバージョン変更してません）

（mathはNewlib prebuild使用）

コード

void setup()
{
 Serial.begin(9600);
 delay(1000);
 int i,j,k,z;
 float distance,ztemp;
 for(i=0;i<=3;++i) {
 for(j=0;j<=3;++j) {
 distance=sqrt((3.5-(float)i)*(3.5-(float)i)+(3.5-(float)j)*(3.5-(float)j));
 ztemp=3.5+sin(distance/1.3+(float)k/3.14159)*3.5;
 z=(int)(ztemp+0.5);
 Serial.print(z);Serial.print(",");Serial.print(i);Serial.print(",");Serial.print(j);
 Serial.println();
 Serial.print(z);Serial.print(",");Serial.print(7-i);Serial.print(",");Serial.print(j);
 Serial.println();
 Serial.print(z);Serial.print(",");Serial.print(i);Serial.print(",");Serial.print(7-j);
 Serial.println();
 Serial.print(z);Serial.print(",");Serial.print(7-i);Serial.print(",");Serial.print(7-j);
 Serial.println();
 Serial.println();
 delay(100);
 }
 }

}

void loop(){
}

結果

1,0,0
1,7,0
1,0,7
1,7,7

3,0,1
3,7,1
3,0,6
3,7,6

4,0,2
4,7,2
4,0,5
4,7,5

5,0,3
5,7,3
5,0,4
5,7,4

3,1,0
3,6,0
3,1,7
3,6,7

5,1,1
5,6,1
5,1,6
5,6,6

6,1,2
6,6,2
6,1,5
6,6,5

7,1,3
7,6,3
7,1,4
7,6,4

4,2,0
4,5,0
4,2,7
4,5,7

6,2,1
6,5,1
6,2,6
6,5,6

7,2,2
7,5,2
7,2,5
7,5,5

7,2,3
7,5,3
7,2,4
7,5,4

5,3,0
5,4,0
5,3,7
5,4,7

7,3,1
7,4,1
7,3,6
7,4,6

7,3,2
7,4,2
7,3,5
7,4,5

5,3,3
5,4,3
5,3,4
5,4,4

Okamiya Yuuki over 3 years ago in reply to Okamiya Yuuki

がじぇるね岡宮です。

試しにk=3として、同じ演算をしてみましたが、結果としては同じでした。RX621は持っていないため、やれてません。あまり解決の糸口にならずすみませんが、とりあえずの実験結果でした。

CCRX

GCC

0,0,0
0,7,0
0,0,7
0,7,7

0,0,1
0,7,1
0,0,6
0,7,6

1,0,2
1,7,2
1,0,5
1,7,5

2,0,3
2,7,3
2,0,4
2,7,4

0,1,0
0,6,0
0,1,7
0,6,7

2,1,1
2,6,1
2,1,6
2,6,6

3,1,2
3,6,2
3,1,5
3,6,5

4,1,3
4,6,3
4,1,4
4,6,4

1,2,0
1,5,0
1,2,7
1,5,7

3,2,1
3,5,1
3,2,6
3,5,6

5,2,2
5,5,2
5,2,5
5,5,5

6,2,3
6,5,3
6,2,4
6,5,4

2,3,0
2,4,0
2,3,7
2,4,7

4,3,1
4,4,1
4,3,6
4,4,6

6,3,2
6,4,2
6,3,5
6,4,5

7,3,3
7,4,3
7,3,4
7,4,4