VPF::Странная реализация деления в std::complex

borisbn

Дата 22.1.2014, 18:33 (ссылка)

(нет голосов)

Загрузка ...

Эксперт

Профиль
Группа: Завсегдатай
Сообщений: 4875
Регистрация: 6.2.2010
Где: Ростов-на-Дону

Репутация: 22
Всего: 135

Здравствуйте.

Сегодня рылся в исходниках STL в реализации Builder'а и нашёл такую штуку:

Код

 inline complex<double>& 
  complex<double>::operator/= (double rhs)  
  {
    double denom = rhs*rhs;
    double re    = (__re*rhs)/denom;
    double im    = (rhs*__im)/denom;
    __re         = re;
    __im         = im;
    return *this;
  }

Вопрос(ы):
1) Зачем понадобилось сначала получать квадрат rhs, а затем делить на него, умножив на rhs ? Почему не сделать просто (как в gcc)

Код

  template<typename _Tp>
    complex<_Tp>&
    complex<_Tp>::operator/=(const _Tp& __t)
    {
      _M_real /= __t;
      _M_imag /= __t;
      return *this;
    }

?
2) Почему с re умножение идёт в порядке __re * rhs, а с im - наоборот ?

Спасибо.

--------------------

Женщины отличаются от программистов тем, что у них чары состоят из стрингов

tzirechnoy

Дата 22.1.2014, 19:03 (ссылка)

(нет голосов)

Загрузка ...

Эксперт

Профиль
Группа: Завсегдатай
Сообщений: 1173
Регистрация: 30.1.2009

Репутация: 2
Всего: 16

Скопировали из комплексного деления, видимо. Приравнять сомножытели с rhs.im к нулю -- догадались, сократить на rhs -- нет.

Можэт, правда, хотели bug-to-bug compatibility с делением на complex(rhs, 0) -- но это вряд ли. Странное это дело -- спецыально вкрячивать ошыбку.

akizelokro

Дата 23.1.2014, 02:39 (ссылка)

(нет голосов)

Загрузка ...

Крокодил

Профиль
Группа: Участник
Сообщений: 761
Регистрация: 30.7.2007

Репутация: 1
Всего: 5

Стандартная формула так записывается (из учебника, где просто мнимая часть делителя = 0)

Почему не сделать как в gcc? Потому что есть формула расчёта комплексных чисел (да и вообще желательна чёткая типизация, раз понятно, что такое комплексные числа и какие на сегодняшний день есть варианты работы с ними). Если кто-нибудь придумает третий вид чисел и пойдёт их шпарить в таком делении, будут интересные эффекты.

--------------------

a = a + b; b = a - b; a = a - b;

borisbn

Дата 23.1.2014, 08:25 (ссылка)

(нет голосов)

Загрузка ...

Эксперт

Профиль
Группа: Завсегдатай
Сообщений: 4875
Регистрация: 6.2.2010
Где: Ростов-на-Дону

Репутация: 22
Всего: 135

Цитата(tzirechnoy @ 22.1.2014, 19:03

)

Скопировали из комплексного деления, видимо. Приравнять сомножытели с rhs.im к нулю -- догадались, сократить на rhs -- нет.

Блин... А точно...
Вот говорили мне: "Не ищи зебр, если видишь следы копыт"...
Я-то думал, что это нужно для какого-нибудь приведения степеней или ещё для чего-то подобного.
Спасибо за ответы.
Закрываю.

--------------------

Женщины отличаются от программистов тем, что у них чары состоят из стрингов

0 Пользователей читают эту тему (0 Гостей и 0 Скрытых Пользователей)
0 Пользователей:
« Предыдущая тема \| C/C++: Общие вопросы \| Следующая тема »