1

以下のテンプレートの関数を、

complex<double> c = u16todbm<complex<double> >(10, 32768.0, 0.0);

と呼ぶと問題なくコンパイルできますが、

complex<short> c = u16todbm<complex<short> >(10, 32768.0, 0.0);

に変更すると、*+を使用している行で、

error: no match for ‘operator*’
error: no match for ‘operator+’

が発生してコンパイルが出来ません。

・呼ばれる側

template<typename samp_type>
samp_type u16todbm(samp_type a, double add, double mul)
{
  // scale offset
  a = a * mul;
  a = a + add;

  return a;
}

・呼ぶ側

int main(int argc, char* args[])
{
  complex<short> c = u16todbm<complex<short> >(10, 32768.0, 0.0); // コンパイルエラー
  complex<double> c = u16todbm<complex<double> >(10, 32768.0, 0.0); // コンパイル通る
}

complexとdoubleの*や+が定義されていないのか?と思いますが、実装方法がわかりません。
お教え願えないでしょうか。

3
  • 2
    std::complex<> で浮動小数点型 (float, double, long double) でない型を使った場合の動作は未定義らしいですが、お使いの処理系は大丈夫でしょうか? 参考: stackoverflow.com/q/11108743/4937930
    – yaegashi
    Commented 2016年11月14日 4:47
  • 1
    unspecified の訳は未規定ですね。未定義とは違います。とりあえずウチの [email protected] でも operator* と operator + が無いのエラーになりました。
    – 774RR
    Commented 2016年11月14日 5:43
  • コメントありがとうございます。こちらはg++5.3.1-1ubuntuです。再現できて良かったです。
    – cpphoge
    Commented 2016年11月14日 9:13

2 件の回答 2

1
int a = 1.0;

とはかけない言語仕様(浮動小数点数を整数へ変換する際に切り捨てや丸めなどいくつかの方法があるため一意に変換することができず暗黙変換が定義されない)なのでcomplex<T>でTが整数型の場合は四則演算可能な型はTへの暗黙変換可能な型に対してしか定義されていない(定義すべきでないから)と思います。complexの定義を見ていないので推測ですが厳密にはTの値を実部と仮定するTからcomplex<T>への暗黙変換が定義されているように思います。
以上からu16todbmは以下のような使い方が適切だと思います。

u16todbm<complex<short>>(10, 32768, 0);

あえて浮動小数点数を指定したいのであれば整数への変換を明示的に実装したu16todbm相当の別のテンプレート関数を定義すべきと思いますが、実装してみると若干不格好なものになりますね。

template<typename samp_type, typename cast_type>
samp_type u16todbm_cast(samp_type a, double add, double mul)
{
  a = a * (cast_type)mul;
  a = a + (cast_type)add;
  return a;
}

手元にC++環境がなく確認できてませんので間違いなどあった場合はご容赦ください

1
  • ありがとうございます。教えて頂いたu16todbm_castでコンパイルは通りました。ありがとうございます!
    – cpphoge
    Commented 2016年11月14日 9:17
1

short に (16bit short であるとして)
- 32768 を加えると符号付整数オーバーフローは未定義動作であるとか
- 0.5 を乗じても精度落ちするとか
- abs()sqrt() を使っているところで精度落ちしないよう工夫しているコードが
 浮動小数点数しか考慮してないとか

いろいろあって std::complexfloat double long double しか規定していません。

それら以外の型については unspecified つまり未規定とされています。
std::complex<short> の動作は未規定、つまり
- std::complex<short> を使うこと自体は間違っていないが
- 得られる結果は誰も定義していない、すなわち
- 期待通りの結果が得られるかもしれない
- 期待はずれの結果が得られるかもしれない

std::complex<short> を演算するとコンパイルエラーになるように
gcc はわざとそういう operator を提供していないんでしょう。

mul/add である double は複素数の実部だけを与えて虚部は 0 でよいということなら

template<typename T>
T u16todbm(T a, double add, double mul) {
    a = a * T(mul);
    a = a + T(add);
    return a;
}

でいけます。
これを std::complex<short> に対して使って意味があるかは別問題。
途中であふれたり精度落ちして期待通りの結果になるとは思えません。

int main() {
    std::complex<short> c(1, -2);
    std::cout << u16todbm(c, 1, 2) << std::endl; // おそらく期待通り
    std::cout << u16todbm(c, 0.5, -0.5) << std::endl; // たぶん期待外れ
}

joke code でこんなの書いたことはあります。何とか期待通り。

int main() {
    std::complex<std::string> a("Hello", "Alice");
    std::complex<std::string> b("World", "Bob");
    std::cout << a+b << std::endl;
}
1
  • ありがとうございます。確かにshort型にした場合、演算途中で値が0になることが分かりました。途中でdouble精度で計算をしないとだめでした。精度を念頭にコードを書いてみます。ありがとうございました。
    – cpphoge
    Commented 2016年11月14日 9:18

この質問に回答するには、ログインする必要があります。

求めていた回答ではありませんか? のタグが付いた他の質問を参照する。