数値計算総合

**デフォルトの名無しさん** · 2020/03/18(水) 04:44:36.36

プログラミングで使う数値計算のスレです
正規化, 微分積分, 行列式, etc

**デフォルトの名無しさん** · 2020/03/18(水) 18:04:20.14

質問スレで0~255を0~3の範囲に変換したいと質問した
回答は x * 3 / 255 を得られた
これについて自分なりに考えてみる

まず x / 255 でxが0~255の範囲であった時に、値が0~1の範囲に正規化される
つまりxを3倍すれば正規化の範囲も3倍になるということになる

回答を頂く前に (x / (255 * 3)) * 10 という式を自分で考えていたが、回答の式のほうがスマートに実装できた

x / 255という式を最初に発想できていれば自力で解決できたかもしれないのが悔やまれる

**デフォルトの名無しさん** · 2020/03/20(金) 05:19:59.16

数値計算の、既にライブラリとかが整備されてそうなやつ(2次方程式の解とか)を、
自前で書かないといけない場合って、
ライブラリの自作やアルゴリズムとかの勉強を除くと残ってるのだろうか？

ライセンス絡みで、仕事では使えないから自分で書くとか、
仕様が自分の用途にはoverkillだからサブセットを書くとかは
あるかも知れないけど。

**デフォルトの名無しさん** · 2020/03/20(金) 09:46:28.25

>>2
>回答は x * 3 / 255 を得られた
x * 4 / 256 (= x / 64)の間違いだよね？

**デフォルトの名無しさん** · 2020/03/20(金) 10:02:06.25

>>4
ごめ、>>2は間違い
質問スレで貰った回答は x * 4 / 256 で合ってる

**デフォルトの名無しさん** · 2020/03/20(金) 13:14:46.24

xが0なら2を1なら－ 1を返すにはどうやる?

**デフォルトの名無しさん** · 2020/03/20(金) 21:26:13.31

2-3*x

**デフォルトの名無しさん** · 2020/03/20(金) 21:28:20.69

xが0または1しか取り柄ないなら
x==0?2:-1

**デフォルトの名無しさん** · 2020/03/22(日) 14:38:45.70

スパコン界隈はまだFortranが現役なんだっけ？
最高齢かつ現役ってスゲーな

**デフォルトの名無しさん** · 2020/03/23(月) 18:53:04.48

2点の場合は線形近似でぴったりいきますね

**デフォルトの名無しさん** · 2020/04/01(水) 03:27:42.30

-1 0 1 0 -1 0 1 ...

みたいな数列が欲しいんだけど、どうすればいいんだろう

-1 + time % 3

だと -1 0 1 -1 0 1 ... という循環になるけど

**デフォルトの名無しさん** · 2020/04/01(水) 06:41:46.41

>>11

if文使うとこうなるけど数学的にスマートな方法がわからない

if x % 2 == 0:
　y = 0
elif x % 4 == 1:
　y = 1
elif x % 4 == 3:
　y = -1

**デフォルトの名無しさん** · 2020/04/01(水) 07:10:56.36

>>11
数学的にきれいなのは三角関数でπ/2ずつ増やすことだろうけど無駄に計算量が増えるだけだな。
以下のでどう？
int x[4]={-1,0,1,0};
int a=x[i%4];

**デフォルトの名無しさん** · 2020/04/01(水) 07:27:36.92

>>13
おお、なるほど
三角関数使うのは俺も考えたんだけど、こっちのほうがスマートだね

**デフォルトの名無しさん** · 2020/04/02(木) 17:45:51.29

>>9
C系はポインタのエイリアス問題があるからね。
restrictができたけど

**デフォルトの名無しさん** · 2020/04/02(木) 18:39:18.20

密度が均一だと仮定した場合の地下H(m)における重力

**デフォルトの名無しさん** · 2020/04/02(木) 19:13:11.13

1-abs(i%4-2)

**デフォルトの名無しさん** · 2020/04/02(木) 19:17:27.60

>>16
質問したいのか何なのか、もう少し日本語を頑張れ

**デフォルトの名無しさん** · 2020/04/02(木) 19:39:13.42

これでもできる
(i%4-1)(i%4+1)(3-i%4)/3

**デフォルトの名無しさん** · 2020/04/02(木) 20:05:36.09

>>17,19
>>17これいいね
計算量も少なそう

**デフォルトの名無しさん** · 2020/05/07(木) 09:29:58.63

x=f(x)で、f(x)が複雑でx=定数の形にするのが大変なんで、数値計算でxの値をだしたいんですが、いい方法は有りますでしょうか

**デフォルトの名無しさん** · 2020/05/08(金) 09:42:46.66

>>21
g(x)=f(x)-xとしてg(x)=0を解けばいいじゃん
g(x)=0を解くライブラリなんてその辺にいくらでも転がってるでしょ
gslとか

**デフォルトの名無しさん** · 2020/05/08(金) 14:22:46.68

>>21
計算式や条件がわからんと

**デフォルトの名無しさん** · 2020/05/11(月) 19:50:17.40

良スレに育てよ

**デフォルトの名無しさん** · 2020/06/18(木) 05:19:02.39

イージングの公式
値 += (目標値 - 現在の値) * 減速値

**デフォルトの名無しさん** · 2020/06/18(木) 10:45:06.58

>>4->>5
なんで？0~255を0~3の範囲に変換したいのに
x=255を代入したらさ
255×3÷255=3は良いとして
255×4÷256= 3.984375だよ？3以下じゃない
同じく255÷64= 3.984375だよ？
0~255を0~3の範囲に変換出来てないじゃない

**デフォルトの名無しさん** · 2020/06/18(木) 10:51:10.14

>>26
ここはプログラム板の数値計算のスレだが、何で"/"記号を普通の算数の割り算だと思ったの？