C++相談室 part145

**デフォルトの名無しさん** · 2019/09/13(金) 17:13:24.60

C++に関する質問やら話題やらはこちらへどうぞ。
ただし質問の前にはFAQに一通り目を通してください。
IDE (VC++など)などの使い方の質問はその開発環境のスレにお願いします。

前スレ
C++相談室 part144
https://mevius.5ch.net/test/read.cgi/tech/1563769115/

このスレもよろしくね。
【初心者歓迎】C/C++室 Ver.105【環境依存OK】
https://mevius.5ch.net/test/read.cgi/tech/1556142878/

■長いソースを貼るときはここへ。■
　http://codepad.org/
　https://ideone.com/

[C++ FAQ]
https://isocpp.org/wiki/faq/
http://www.bohyoh.com/CandCPP/FAQ/ (日本語)

----- テンプレここまで -----

**デフォルトの名無しさん** · 2019/10/03(木) 21:59:37.34

>>294 エラーメッセージとか出てないの？

**デフォルトの名無しさん** · 2019/10/03(木) 22:06:09.67

>>294
clangさんが教えてくれた
std::string hogehoge() { return A<_T3>::template get<std::string>();}

**デフォルトの名無しさん** · 2019/10/03(木) 22:13:42.78

ありがとうございました。
通りました。

つか、そんなtemplate文の使い方初めて知った・・

**デフォルトの名無しさん** · 2019/10/03(木) 22:31:15.90

この場合Aがテンプレートなので、Aの先のgetが何者かコンパイラは判断つかないので
"<" が比較演算なのかテンプレートのかっこなのか判断つかないから
らしい
clangは教えてくれるんだから判断ついてるわけだが・・・
まったく醜い文法を持った言語だよ

**デフォルトの名無しさん** · 2019/10/03(木) 22:49:50.59

この場合のエラーメッセージ、GCCはほんとに何言ってるのか分からないのを出すのでclang様様

**デフォルトの名無しさん** · 2019/10/04(金) 01:11:35.70

C++テンプレートのプログラミングでのエラーを読み解く訓練をしているうちに
実世界での初心者が投げてくる質問へのエスパー力がやしなわれる
これ豆

**デフォルトの名無しさん** · 2019/10/04(金) 02:30:32.92

そのうちエラーメッセージを見ずにエラーを修正できるようになる

**デフォルトの名無しさん** · 2019/10/04(金) 02:54:08.89

ちなVC2015か2017か忘れたけど、一時期までのVCは
その場合のtemplateを書くと逆にエラーにされてた

**デフォルトの名無しさん** · 2019/10/04(金) 03:18:41.39

goで変数用いて
*******
*****
***
*
***
*****
*******
って出力するにはどうすりゃいいのっと

**デフォルトの名無しさん** · 2019/10/04(金) 03:19:28.61

あ、↑のは砂時計みたいな形にしたい

**蟻人間** ◆T6xkBnTXz7B0 · 2019/10/04(金) 04:12:53.14

>>303
二重のループか仮想画面

**デフォルトの名無しさん** · 2019/10/04(金) 04:19:18.14

>>305
申し訳ないが初心者ゆえ、数字の二重ループはできるけどそれを文字列にする仕方がわからず

**デフォルトの名無しさん** · 2019/10/04(金) 05:45:45.52

goスレ行け

**デフォルトの名無しさん** · 2019/10/04(金) 08:20:48.47

>>298
なぜか偶然ですが、このスレを見る前から昨日からその辺に関係するパーサーを
作っていて悩んでいたところなんですが、
多くの言語では、+,-,*,/,<,> などの演算子が混ざったいわゆる「数式」は文脈に依存せずに
上位のコンパイル作業に入る前に、いったん「数式ツリー」に直すことが出来ます。
ここで「文脈に依存せず」というのは、変数名、関数名、型名などの宣言情報を
調べなくても済むという（いわゆる意味解析の結果の情報を必要としないと言う
）ことです。
ところが、template機能が入った後のC++では、A<B と書いたとき、
< がテンプレートの記号なのか比較演算子なのかは A がテンプレート名であるか
どうかによって決まります。
Bが型名であることが明確な場合には文脈によらず < はテンプレート記号だと断定
できます。ところが、Bに 1, 456, 3.1415 のような数値や、"Hello" などの文字列
が来ている場合には、（その場では）判断が付きません。
文脈に依存せずに判断を付けたいなら、もっと先まで呼んで、< に対応する >
があるかどうかを調べるといける可能性があります。

1. "<"に対応する">"が存在するかどうか調べれば、これまでの宣言情報を使わなくても
　トークン・ツリーを作ることが出来る。トークン・ツリーと言ったのは、数式とは
　限らないからです。例えば、A<B> x[100]; などとすることがありますが、A<B>
　は型名であり、数式では有りません。ただし、他の多くの言語では数式ツリーは
　作ることが多いのですが、トークン・ツリーはまず作りません。
2. 1の方法を使わないなら、宣言情報を調べて、Aがtemplate名であるかどうかを調べる
　必要がある。ところが、Q::A<B> などのような場合があるので、そんなに
　単純な作業でく、Q::A が型名なのか、変数名なのかはちゃんと調べる必要がある。

**デフォルトの名無しさん** · 2019/10/04(金) 08:25:56.64

>>308
補足
「Bが型名であることが明確な場合には文脈によらず < はテンプレート記号だと断定
できます。」
の部分ですが、B<CPerson(実引数列) などとした場合、CPersonは型名ですが、
CPerson(実引数列) 全体は、型名ではなくCPersonクラスの引数つき
コンストラクタを呼び出した結果の「一時オブジェクト」です。なので、
Bが変数名の場合は、< は比較演算子の可能性があります。

**デフォルトの名無しさん** · 2019/10/04(金) 08:33:37.62

>>308
誤：なぜか偶然ですが、このスレを見る前から昨日からその辺に関係するパーサーを
正：なぜか偶然ですが、このスレを見る前の昨日からその辺に関係するパーサーを

誤：文脈に依存せずに判断を付けたいなら、もっと先まで呼んで、< に対応する >
正：文脈に依存せずに判断を付けたいなら、もっと先まで読んで、< に対応する >

誤：単純な作業でく、Q::A が型名なのか、変数名なのかはちゃんと調べる必要がある。
正：単純な作業ではなく、Q::A がtemplate名なのか、変数名なのかはちゃんと調べる必要がある。

**デフォルトの名無しさん** · 2019/10/04(金) 08:37:14.54

>>225
>>224=228や>>226=219みたいなどうみても初心者なヤツにプロが敬語使う必要ないと思うよ
C++スゲー=俺スゲーみたいな明らかにプログラミングに向いてない精神構造のクソバカが多く居るからねここ

**デフォルトの名無しさん** · 2019/10/04(金) 08:41:34.48

>>306
一般論として、文字列の長さを len として、文字列型の変数を str としたとき、
for ( len 回のループ ) {
　str = str + "*"
}
で行けます。

**デフォルトの名無しさん** · 2019/10/04(金) 09:33:34.23

>>311
ようクソバカ

**デフォルトの名無しさん** · 2019/10/04(金) 09:41:54.97

>>313
すまんな図星だったか？
邪魔だから出てけよ

**デフォルトの名無しさん** · 2019/10/04(金) 11:43:51.24

>>308
C++言語をパースする方法として、結論的には、2.ではダメで、1.でなくてはならないようです。例として、A<実引数列>::B の場合を考えると、template は実体化する
際の実引数によって展開される内容は大幅に変わりますので、
A<int>::B とした場合には、B は template 名で、
A<string>::B とした場合には、B は変数名であることもありえます。すると、
A<int>::B<C> の場合の B<C> は template 実体化ですが、
A<string>::B < C の場合は、変数B と変数 Cに対しての比較演算となります。
なので、もし、>>308 の 2.のやり方でパース使用とすると、Aがtemplate名であると
分かっただけでは B が「何なのか」を結論付けるには不十分で、A<実引数列>の
実引数列をほとんど完全に調べ上げてから、その方名に応じた template A の
定義を探して template の overload 解決をしてから、A<実引数列> 部分(のclass)を
完全に「特定」する必要があります。ここまでやるとなると、コンパイラの層分けが
上手く出来ないことになります。

**デフォルトの名無しさん** · 2019/10/04(金) 12:03:58.19

>>315
A::B くらいの場合なら、意味解析の情報だけを元に、Bが何なのかを特定をすることは
可能です。ところが、A<1+2,decltype(x*y)>::B の様な場合、x, y の意味論的な型
を調べた後、x*yに適応可能な operator*() 関数が定義されて言うかどうかも調べ、
定義されていれば、その戻り値の型を調べてx*y の型として、A<・・・> の部分が
結局何なのかを特定しなくてはなりません。今述べた x*y の部分の処理は、
トークン解析層(?)より上位のコンパイル層で行う処理です。なのでこのやり方
だと、コンパイラ内部の「層」の切り分けが難しくなります（技術的に不可能な
わけではありませんが。）。

**デフォルトの名無しさん** · 2019/10/04(金) 12:07:27.45

>>315
誤：なので、もし、>>308 の 2.のやり方でパース使用とすると、Aがtemplate名であると
正：なので、もし、>>308 の 2.のやり方でパースしようとすると、Aがtemplate名であると

誤：実引数列をほとんど完全に調べ上げてから、その方名に応じた template A の
正：実引数列をほとんど完全に調べ上げてから、その型名に応じた template A の

>>316
誤：を調べた後、x*yに適応可能な operator*() 関数が定義されて言うかどうかも調べ、
正：を調べた後、x*yに適応可能な operator*() 関数が定義されているかどうかも調べ、

**デフォルトの名無しさん** · 2019/10/04(金) 12:36:35.66

>>312
ありがとうございます！
やってみます

**デフォルトの名無しさん** · 2019/10/04(金) 12:41:17.15

c++って文法に曖昧さがないことどうやって保証してんだろ
形式的なアプローチ無理だと思うんだが

**デフォルトの名無しさん** · 2019/10/04(金) 13:05:22.64

いや、そんな保証は端っから諦めてるだろ
引数なしの関数宣言とデフォコン呼び出しの曖昧さなんか
長年放置してたのをC++11でようやく対応したし
関数ポインタにキャストするoperatorはtypedef使えとか
言語としての完全性なんか重視してない

**デフォルトの名無しさん** · 2019/10/04(金) 13:28:40.80

>>319
曖昧さはあるといわれています。例えば、CPersonというクラスがあったとき、関数内で
CPerson person();
と書いた場合、person という名前で、戻り値の型が CPerson の関数のプロトタイプ宣言なのか、
CPerson クラスの person という変数名(オブジェクト名)の定義なのかの曖昧さがあります。
後者は、CPerson person("name", 25, MALE); などと同じ系統で、実引数が全く無い場合
に相当します。
また、template class の場合に、
A<B<・・・>>
と書くと、>> が右シフト演算子に解釈されてしまうので、回避するために空白を1つ入れて
A<B<・・・> >
としなければならなかった(過去形)時代も有ります。
しかし、パーサーを作る側からすれば、">>" を ">" 2個だと解釈するのはかなり大変な
労力が必要でです。

また、x < y を比較演算のつもりで A< x < y >　などと書くと、A<・・・>
の中に x という template class 名に対しての x<y> が入っていると解釈
されてしまい、> が足りないというエラーになるかもしれません。
回避するには、
A< (x < y) >
と書くと良いと思われます。

**デフォルトの名無しさん** · 2019/10/04(金) 14:48:30.02

C++というより３Dプログラミングの質問で申し訳ないのだけど、
行列ライブラリglmの以下のコードがbugってる気がするので詳しい人がいたら確認して欲しい。

template<typename T, qualifier Q>
GLM_FUNC_QUALIFIER qua<T, Q> rotate(qua<T, Q> const& q, T const& angle, vec<3, T, Q> const& v)
{
vec<3, T, Q> Tmp = v;

// Axis of rotation must be normalised
T len = glm::length(Tmp);
if(abs(len - static_cast<T>(1)) > static_cast<T>(0.001))
{
T oneOverLen = static_cast<T>(1) / len;
Tmp.x *= oneOverLen;
Tmp.y *= oneOverLen;
Tmp.z *= oneOverLen;
}

T const AngleRad(angle);
T const Sin = sin(AngleRad * static_cast<T>(0.5));

return q * qua<T, Q>(cos(AngleRad * static_cast<T>(0.5)), Tmp.x * Sin, Tmp.y * Sin, Tmp.z * Sin);
}

クオータニオンqをベクタ軸ｖ周りにAngleRad回転させる関数のソースらしいんだが、
最後の行掛け算の順番が逆のような気がする。
qua<T,Q>(...) * qが正解だと思うんだ。

glmはファイル構成がバージョン間で結構派手に変更されてるんで、
最近のバージョンではext/quaternion_transform.inlに上のコードが含まれてる。
古いバージョンではgtx/quaternion.inlにあったような気がする。

**デフォルトの名無しさん** · 2019/10/04(金) 14:55:11.92

>>322
あるライブラリでのq1 * q2 という書き方が、q2 * q1
がどういう意味に解釈されるかは、

**デフォルトの名無しさん** · 2019/10/04(金) 15:01:05.28

>>322
(>>323 は誤投です。)
あるライブラリでの q1 * q2 という書き方が、別のライブラリでは
q2 * q1 と等価なことが有り得て、それはベクトルや行列計算全般に
言える記法の違いなのでこのソースだけではなんとも言えない。
ベクトルの場合、成分表示を1列の縦行列と考えるか1行の横行列
と考えるかによって、賭け全の左右が逆さまになる。
その結果、成分表示で書かない表記も、左右どちらかで書くかは
好みの問題となる。
特に Quaternion の場合も同様。

**デフォルトの名無しさん** · 2019/10/04(金) 15:02:02.33

>>324
誤：と考えるかによって、賭け全の左右が逆さまになる。
正：と考えるかによって、掛け算の左右が逆さまになる。

**デフォルトの名無しさん** · 2019/10/04(金) 15:45:42.88

OpenGL系は縦ベクトルだね
D3Dを前提にしたライブラリとは逆になる(行列も転置

**デフォルトの名無しさん** · 2019/10/04(金) 15:58:00.91

>>326
数学や物理学では縦ベクトルが多いが、D3Dはそれに従ってない。

**デフォルトの名無しさん** · 2019/10/04(金) 16:08:07.51

テンプレが糞だと思う瞬間はいくつもあるが
その一つはネストしてるときの閉じ括弧
>> じゃだめで > > ってわざわざ分けないといかんうざい

**デフォルトの名無しさん** · 2019/10/04(金) 16:08:35.30

＞数学や物理学では縦ベクトルが多いが、D3Dはそれに従ってない。
みたいだね
D3Dからそれ以外に移るとみんなここで引っかかるけど
実は行列周りは(大文字は行列として)
OpenGLに合わせるとA*B*C*xとなるのを
D3Dではx*C*B*Aと書けた方が、左からの計算が全部結果がベクトルなので計算量が減るというメリットはある
MSはC++での書きやすさと効率を考えたのかもしれない

**デフォルトの名無しさん** · 2019/10/04(金) 16:09:02.58

>>326-327
縦横ってより右から掛けるか左から掛けるかだと思ってた
意味は一緒なんだろうけど

**デフォルトの名無しさん** · 2019/10/04(金) 16:11:58.82

>>328
C++11からは連続して書けるようになったで

**デフォルトの名無しさん** · 2019/10/04(金) 16:14:31.61

C++03いやC++98？のまま更新止まってるやつ
俺が想像するより遙かに多いのかもしかして

もうC++11でさえ要注意な旧規格になってるんだが

**デフォルトの名無しさん** · 2019/10/04(金) 16:16:30.52

>>329
メモリ上の順序も変わるぞ
SIMDとかで計算されるときの効率にも影響する
必ずしもD3D(横ベクトル方式)が高速化に都合が良いとは言えないんじゃないか

行列の各要素のメモリ上の並びって
glの行列AとD3Dの行列A'はいちいち転置しなくても結局同じになるんだっけ

**デフォルトの名無しさん** · 2019/10/04(金) 16:22:07.03

ああ確かにメモリの中身は同じだった
格納方式も違うから結果的に同じになるんだよね
そこが余計ややこしい

**デフォルトの名無しさん** · 2019/10/04(金) 16:29:39.34

>>329
行列計算は、縦ベクトル方式、横ベクトル方式でも、（それぞれで）「結合則」があるので、
A*B*C*x = A*(B*(C*x))
でもあり、
A*B*C*x = ((A*B)*C)*x)
でもあり、どっちで計算しても結果は同じになる。

また、3D計算の場合、点の座標がいくつあっても、行列を一度計算しておくと、
使いまわしできるので、先に行列の方を計算しておいて、最後に座標の変換
を行うと、点が2個以上あるときには効率が高い。

**デフォルトの名無しさん** · 2019/10/04(金) 16:31:49.50

>>332
C++98くらいから、実際問題上使いたくない機能ばかりが追加されるような
傾向があったのに、Upgrade 版は適用外になったりしたりして買う機を逃した。

**デフォルトの名無しさん** · 2019/10/04(金) 16:38:53.10

>>329
>MSはC++での書きやすさと効率を考えたのかもしれない
>>335 に既に書いたけど、OpenGL 形式でも、右結合でも、
左結合でも、好きな方で計算は出来る。だから、行列部分を
先に計算しても良いし、座標部分から計算しても良い。
座標部分から計算していくと一見計算量が減るように思える
かも知れないが、点の個数が多くなると、一度だけ行列部分の
積を予め計算しておいて、その結果を全ての点の変換部分に
遣うと効率が劇的に向上するので、むしろ、効率が良いのは逆。

**デフォルトの名無しさん** · 2019/10/04(金) 16:40:53.55

>>335
確かに(OpenGLでも上で書いたような式が使えるライブラリ使うとして)括弧でくくったりして優先順位変えればいいし
言う通り行列はまとめればいいんだけど
どちらにもメリットデメリットあると言いたかった

**デフォルトの名無しさん** · 2019/10/04(金) 16:51:20.07

>>338
実は、アメリカで有名な 3D Graphics の本が、横ベクトル方式を採用していた
事が関係しているかもしれない。要はその本が数学や物理の標準とは何故か
逆の記法を採用していたということ。

**デフォルトの名無しさん** · 2019/10/04(金) 17:00:49.77

>>337
GL
(A*B*C)*(x0(縦),x1(縦),x2(縦),x3(縦), ... , xN(縦))(横)
D3D
(x0(横),x1(横),x2(横),x3(横), ... , xN(横))(縦)*(C*B*A)
が最速やね

◆QZaw55cn4c · 2019/10/04(金) 19:54:25.76

>>335
行列って結合則は約束されてましたっけ？

**デフォルトの名無しさん** · 2019/10/04(金) 20:17:27.71

当たり前だろ。
そんなレベルの知識しかないからフォントライブラリ書けないんだよ

◆QZaw55cn4c · 2019/10/04(金) 20:27:59.81

>>342
では証明してください（キリッ）

**デフォルトの名無しさん** · 2019/10/04(金) 21:24:32.75

>>343
ベクトル部分はおいておくとして、まず正方行列の部分だけに限定すれば
A(BC)=(AB)C　　・・・(1)
が証明できればどんな場合でも結合側が成り立つことが数学的帰納法で
証明できます。なので、A,B,C が正方行列の時に(1)を証明すれば全体の証明が
ほぼ終わります。
【(1)の証明】
行列の積の定義により
　(AB)ij=Σ_{k=1}^n A_{ik} B_{kj}　　・・・(2)
です。なので、
{(AB)C}mn = Σ_l(AB)_{ml}C_{ln}
　　　　　= Σ_k Σ_l A_{mk} B_{kl} C_{ln}　・・・(3)
です。全く同様に、
{A(BC)}mn = Σ_k A_{mk} (BC)_{kn}
　　　　　= Σ_k A_{mk} Σ_l B_{kl} C_{ln}
　　　　　= Σ_k Σ_l A_{mk} B_{kl} C_{ln}　・・・(4)
となり、(3), (4) が一致することから、
　{A(BC)}mn = {(AB)C}mn
が言えます。これが (1) に他なりません。　　　　(Q.E.D)

次に、このように3つの行列の場合ではなく、一般の個数の行列の場合、
最初に書いたように数学的帰納法を使います。それは、あなた自身で
お考えください。考える力の練習になります。

**デフォルトの名無しさん** · 2019/10/04(金) 22:10:52.87

doubleの値を++aした時にaの値は+1と規定されてる？

**デフォルトの名無しさん** · 2019/10/04(金) 23:19:26.96

>>345
boolen型以外では、++x は正確に x+=1 と等価です。
x+=1 は、正確に x = x + 1 と等価です。x がポインタ型の場合は、
BYTE 単位で見ると object size 単位で増加しますが、double/float/int/short/char
などの「算術型」の場合は、単純に 1 が足されるだけです。

https://en.cppreference.com/w/cpp/language/operator_incdec
For non-boolean operands, the expression ++x is exactly equivalent to x += 1, and the expression --x is exactly equivalent to x -= 1, that is, the prefix increment or decrement is an lvalue expression that identifies the modified operand.

**デフォルトの名無しさん** · 2019/10/04(金) 23:36:52.32

>>346
ありがとうございます。
という事はキャストの時におかしくなっているのかな？

aの値は0.0～100.0程度
for(double i = a; i < 1000.0; ++i){
std::cout << (int)i << "," << std::endl;
}
期待している結果: 1,2,3,4,5...と連番になる
たまに起こる結果: 1,2,3,5,6...と歯抜けになる

**デフォルトの名無しさん** · 2019/10/04(金) 23:42:33.92

丸められただけでは

**デフォルトの名無しさん** · 2019/10/04(金) 23:55:23.10

>>341-344
また5ch名物のしょーもない煽り煽られの流れかとおもったら
ちゃんと勉強になる流れでワロタ

**デフォルトの名無しさん** · 2019/10/05(土) 00:03:15.79

行列の結合則なんて高校で習っただろ
線形代数わからないプログラマって生き残れないと思うわ

**デフォルトの名無しさん** · 2019/10/05(土) 00:05:26.41

>>347
・そのコードは、完全に上記の通りですか？

例えば、ループ回数が1000回より遥かに多くなると、
誤差の都合でそうなることがあるかもしれませんが、1000回くらい
だと、正しい処理系ではその現象は起きないはずです。

**デフォルトの名無しさん** · 2019/10/05(土) 01:17:38.13

キャストする前にlround呼んだらええんちゃうん。

**デフォルトの名無しさん** · 2019/10/05(土) 01:44:04.24

(0.2+0.3)+0.4 == 0.2+(0.3+0.4)の結果がfalseになるのは正しいの？

**デフォルトの名無しさん** · 2019/10/05(土) 01:54:00.97

浮動小数点を勉強しておいで

**デフォルトの名無しさん** · 2019/10/05(土) 02:07:17.64

有効数字1桁の計算で誤差は生じない

**デフォルトの名無しさん** · 2019/10/05(土) 07:04:22.03

10進の0.1を2進表現すると循環小数定期

>>355
バカ乙

**デフォルトの名無しさん** · 2019/10/05(土) 08:03:13.32

>>336
本当、人って色々だな
俺はC++98は未完成感があまりに強くて続きはまだかと待ちかねていた
具体的には例えば右辺値参照だ
一時オブジェクトにいちいちconstがつくのがイヤでイヤで待望のやつがやっと来た

**デフォルトの名無しさん** · 2019/10/05(土) 08:28:45.86

>>337, >>340
ちょっと誤解招くのでやっぱ突っ込んでおくけど、
多数のベクトルを変換するのに行列をまとめておかないと無駄になるのはどっちの方式でも同じだよ
俺が言ったのは(少なくとも同一の行列が)一度限りしか出てこない場合の話
まぁそれを最適化する必要があるのか、またそういう処理を大量にループするなら最適化の方法は他にあるだろうけど
>>333
GL方式でも列優先だから結局d3dと同じ無駄はあるんやで
最近のは知らんけど

**デフォルトの名無しさん** · 2019/10/05(土) 08:49:33.11

>>353
計算誤差です。
(0.2+0.3)+0.4 == 0.2+(0.3+0.4) が偽になっても
(0.2+0.3)+0.4 == (0.2+0.3)+0.4 や
0.2+(0.3+0.4) == 0.2+(0.3+0.4)
は必ず真になることが保障されています。

**デフォルトの名無しさん** · 2019/10/05(土) 08:55:38.66

>>359
コンピュータにおける浮動小数点は、内部表現は10進数ではなく2進数
で表現されており、有効桁数は 10 進数で 8 桁や 15桁などではなく
2進数で xx BIT という風になっています。
0.2, 0.3, 0.4 は、10進数だと、有効桁数が1桁でも完全に区切れが
よく表現できていますが、2進数の表現だと厳密には無限小数になってしまい、
どんなに有効桁数を長くしても厳密には表現できません。そのため計算誤差
が生じるのです。0.5 や 0.25 や 0.125 は 2 進数でも区切れ良く表現できるため、
　(0.125+0.25)+0.5 == 0.125+(0.25+0.5)
は誤差が生じることがないため、必ず真になるはずです。

**デフォルトの名無しさん** · 2019/10/05(土) 10:06:32.18

>>360
浮動小数の計算誤差を考慮して比較する時にstd::abs(a-b) < std::numeric_limit<double>::epsilon()というふうにするのはどう思う？

**デフォルトの名無しさん** · 2019/10/05(土) 11:34:20.93

>>361
誤差の量が計算の順序や書き方によって変わってくるので、
そのようなヘッダに書かれているような誤差定数を使うことは
余り意味がありません。

**デフォルトの名無しさん** · 2019/10/05(土) 11:37:49.72

実用数学

**デフォルトの名無しさん** · 2019/10/05(土) 11:43:08.23

>>361
機械イプシロンはそういうふうに使うものじゃない

**デフォルトの名無しさん** · 2019/10/05(土) 12:09:22.39

cpprefjpが>>361と全く同じ間違いしてたわ馬鹿が湧くのはこいつのせいか
cppreference.comは正しかったからそっち見れ

**デフォルトの名無しさん** · 2019/10/05(土) 12:34:57.85

>>365
x-yの元になっているxやyがどのような計算過程を経てきたかわからないのに
　std::abs(x-y) <= std::numeric_limits<T>::epsilon() * std::abs(x+y) * ulp
でxとyがalmost equalsと言い切ってよい根拠とは、

ちゅか常識的に考えて、x、yがそれぞれ±e、±gの誤差を有するなら
　std::abs(x-y) <= std::numeric_limits<T>::epsilon() * (std::abs(e) + std::abs(g))
あたりの判定に落ち着かねばおかしい（これで正解とは言っていない

なんでstd::abs(x+y)みたいな場合によっては莫大な値になりえる係数を掛けねばならないの？

**デフォルトの名無しさん** · 2019/10/05(土) 13:04:02.75

イプシロンに誤差を掛けても意味ない。
つか、許容誤差の絶対値が与えられているならイプシロン使う必要もない。
cppreference.comのは許容誤差をULP単位で与える場合のやりかたで、
abs(x+y)はxとyの指数のうち大きい方を意味するものだったはず。

**デフォルトの名無しさん** · 2019/10/05(土) 13:12:59.71

>>366
なんか色々根本的に勘違いしてる
cppreference.comのalmost_equal関数の仕事はあくまで与えられたxとyがalmost equalかどうかの判定だよ
> x、yがそれぞれ±e、±gの誤差を有するなら
なにこれ？「誤差」って何に対する何の誤差？

**デフォルトの名無しさん** · 2019/10/05(土) 13:28:00.03

>>367
>イプシロンに誤差を掛けても意味ない。
>つか、許容誤差の絶対値が与えられているならイプシロン使う必要もない。
左様ですなｽﾏﾝｶｯﾀorz

>>368
>なにこれ？「誤差」って何に対する何の誤差？
真の値に対する誤差じゃわﾊﾟｵｰﾝ
真の値に対する誤差に依存しないalmost equal判定とか、使いどころはきわめて限定されるはず、

**デフォルトの名無しさん** · 2019/10/05(土) 14:05:34.77

>abs(x+y)はxとyの指数のうち大きい方を意味するものだったはず。
xが (xの符号)×1.bbbb...×2^m
yが (yの符号)×1.cccc...×2^n
だとしたときに、m>nだとすると
　x + y =
　{ (xの符号)×1.bbbb... + ((yの符号)×1.cccc...×2^-(m-n)) } ×2^m
として計算されるのでだいたい2^m（x≒yなら2^(m+1））という意味か左様か、

**デフォルトの名無しさん** · 2019/10/05(土) 14:08:05.45

んまーulpをxとyの来歴に応じて調整すべき量とするなら
cppreference.comのalmost_equal関数は理解してやらないでもない

**デフォルトの名無しさん** · 2019/10/05(土) 16:03:17.50

だれか、俺に分数を教えてくれないか？
足し算と引き算がマジ解らなくなってるんだ。
https://ideone.com/IfCcGz

一番のネックはネガフラグの扱い。
マイナスとプラスの反転演算でどうかくとスマートだろうか。
だれか教えてプリーズ。

**デフォルトの名無しさん** · 2019/10/05(土) 16:21:25.21

分子に符号もたせりゃいいだろ

**デフォルトの名無しさん** · 2019/10/05(土) 16:25:36.94

>>373
え、それでよかったっけ？

**デフォルトの名無しさん** · 2019/10/05(土) 16:34:31.54

>>350
今習わないよ
俺の時でも選択

**デフォルトの名無しさん** · 2019/10/05(土) 17:18:56.47

一応以下のような感じになった。Thx!
http://coliru.stacked-crooked.com/

**デフォルトの名無しさん** · 2019/10/05(土) 17:19:54.63

http://coliru.stacked-crooked.com/a/5fa6ab838f9ff15a
こっちだった。Orz

んじゃ。

**デフォルトの名無しさん** · 2019/10/05(土) 17:28:26.58

>>377
1. constexprにしてもいいんじゃない？
2. 出力関数作成しよう

template<typename Char>
std::basic_ostream<Char>& operator<<(std::basic_ostream<Char>& os, const Fraction& frac){
return os << frac.GetNumerator() <<'/'<<frac.GetDenominator();
}

**デフォルトの名無しさん** · 2019/10/05(土) 19:08:43.36

>>378
1.できなくはない。機械的に接頭語つけるだけでいいよね。
2.暇だから作るわ。

行ってきまーす。

**デフォルトの名無しさん** · 2019/10/05(土) 19:46:54.47

https://ideone.com/uzd7EF
ただいまー。イデオンで動かせるようにした。Gcdがバグってないといいなぁ。
あとvc2019では意味不明なエラーで動かない。

**デフォルトの名無しさん** · 2019/10/05(土) 19:58:28.51

フィードバックハブに投げておいた。

**デフォルトの名無しさん** · 2019/10/06(日) 10:28:58.03

>>380
マルチ禁止

**デフォルトの名無しさん** · 2019/10/06(日) 13:42:53.69

>>382
ソーリー。

**デフォルトの名無しさん** · 2019/10/07(月) 19:17:43.20

Ｃ＋＋ってなんでヘッダファイルいるんですか
最近マシン早いんだし
全部ソースに書いてインポートすればいいのに

**デフォルトの名無しさん** · 2019/10/07(月) 19:38:58.27

やれば？(クレヨンしんちゃん風に)

**デフォルトの名無しさん** · 2019/10/07(月) 19:52:54.51

来年まで待ってね

**デフォルトの名無しさん** · 2019/10/07(月) 20:58:42.59

>>384
むしろソースがいらなくて全部ヘッダに書けばいい
inline変数とかもできるようになったことだし

**デフォルトの名無しさん** · 2019/10/07(月) 21:05:07.29

>>387
たまに趣味で何でもかんでもtemplateでゴリゴリ実装してるとガチでそれになるww

**デフォルトの名無しさん** · 2019/10/07(月) 21:26:10.24

>>385
ヤレばデキるは魔法の言葉

**デフォルトの名無しさん** · 2019/10/08(火) 09:06:39.07

プリコンパイルヘッダーって使ってます？

使っているプロジェクト見たことないですわ。

**デフォルトの名無しさん** · 2019/10/08(火) 20:03:58.69

>>387
こういう馬鹿野郎が出てくるからメタプロとかテンプレ好き連中ってのは嫌いなんだよ。

**デフォルトの名無しさん** · 2019/10/08(火) 21:04:41.31

C++なんて趣味でやってるんだからテンプレ好きでもええやん

**デフォルトの名無しさん** · 2019/10/08(火) 21:16:41.69

出てくることの何が嫌なんだろうか
はっきり言って病気だ

**デフォルトの名無しさん** · 2019/10/08(火) 21:42:59.95

headerオンリーってのは可搬性に優れた良いものだろ
特にテンプレート多用している奴は、一部関数だけソースファイル分離する意味が薄い

**デフォルトの名無しさん** · 2019/10/08(火) 21:43:38.32

>>391
どんなものにだって馬鹿は出てくるだろうに