数値計算総合

**デフォルトの名無しさん** · 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
ここはプログラム板の数値計算のスレだが、何で"/"記号を普通の算数の割り算だと思ったの？

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

ご教授願います。
Windows7以前の電卓はメモ帳等のテキストエディタに書いた長い計算式をコピペで計算出来たのですが
Windows10になってからは出来ないのでしょうか？
あと、計算式テキストやファイルをコマンドライン引数で渡す方法がありませんでしょうか？
VBSのActiveXで渡してリターンが得られらたら嬉しいんですが
どなたかご存知有りませんか？よろしくお願いします。

**デフォルトの名無しさん** · 2020/06/18(木) 11:08:16.46

>>27
数の値を計らって算出するからじゃないの？
数値計算という言葉には近似計算の歴史も含まれているわけで、
そりゃあ質問スレで質問されたレスを検索してまで読みに行ってまでは無いので目的は知らんけど
intやdoubleのキャストが暗黙とは言えないでしょう？

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

>>29
>>2,4,5のレスを見ても明示されてはいないが、それが整数を対象にした議論で除算も余り切り捨ての(C言語的な)整数の除算であることは、すぐに察しがつくでしょ。
初めに実数を想定して>>26が浮かんだとしても、本人たちが納得しているのだから｢この指摘は正しくないのではないか？ああ、整数だけを想定したものなら成り立つからそういうことなんだろう｣と気づきそうなものだ。
気づかなくて疑問に思うのも別におかしくはないけど、それを>>26のように相手が間違っているはずと決めてかかるような攻撃的な書き方をしなくてもいいんでないの？

**デフォルトの名無しさん** · 2020/06/18(木) 12:14:07.48

>>30
ぁあ「～じゃない」語尾が攻撃的な表現に取られたのね？それは失礼
スレタイに惹かれて来たけど変に絡まれるとはとんだ災難だよ・さようなら

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

そんな事言わずに過疎スレなんだから

**デフォルトの名無しさん** · 2020/06/18(木) 21:54:21.81

C系はintなら暗黙に切り捨てられるし、PythonやJavaScriptなら実数になるよな
この辺は言語によって書き方が変わるし、おそらく>>26の問答はこれから頻出するだろう

**デフォルトの名無しさん** · 2020/06/21(日) 17:26:56.05

物理simをするときに
跳ね返り係数とか摩擦係数とかバネ定数とかって
どこから持ってくるんですか？
リアルの物体で実験して実測値を出すんですか？

**デフォルトの名無しさん** · 2020/06/21(日) 17:30:42.66

>>34
通常は製造メーカーのページかSpecification問い合わせ

**デフォルトの名無しさん** · 2020/06/26(金) 12:17:35.57

数値計算の専門家です
何でも答えます

四則演算までは答えられると思います

**デフォルトの名無しさん** · 2020/06/26(金) 13:12:11.63

0で割ったらダメな理由と
浮動小数点数で 0.0f とか 0.0 を正しく判定する方法を教えてください
float a = 何か;
if(a == 0.0f)
double b = 何か;
if(b == 0.0)

**デフォルトの名無しさん** · 2020/06/26(金) 19:25:52.67

>>37
0で割っても良いです
整数なら例外を発生させる効果があったり
小数なら特殊な値を作る効果があったりします
もちろん環境によって効果は変わります

0.0f 0.0 の比較
何が正しいかは時と場合によります
if (a == 0.) if (a==0.f) で正しい場合も有りますが
そうじゃない場合も有ります

**デフォルトの名無しさん** · 2020/06/27(土) 09:56:05.95

正三角形の高さと重心の高さの比は 3:1 ですが
正四面体の高さと重心の高さの比はなぜ 3:1 ではなくて 4:1 なのですか

**デフォルトの名無しさん** · 2020/06/27(土) 10:18:01.60

>>39
次元の違いスレ違いだから

**デフォルトの名無しさん** · 2020/06/29(月) 01:13:04.00

numerical recipes in Cってまだ現役？

**デフォルトの名無しさん** · 2020/08/02(日) 07:15:28.02

二次元グラフ上で以下の点(A～D)があったとして
;　　　　　Y
;　　　　　|　B
;　　　C　|
Xーーーー+ーーーー
;　A　　　|
;　　　　　|　　　D
ここに入力点としてIを置いた時、
そのIとA～Dのの位置関係を見て各適応割合を算出したいのですが
どうすればできるかな？
例えばIがAとCの間にある時は AC間で割合は変化して、BとDは割合0になる等