Создать функцию на ST

Цитата:

---

Сообщение от zjWlad

Кстати, при последовательном вычислении членов ряда ошибка не возникает. Значение функций SIN_Rad и SIN_Rad2 практически совпадают. По данным fSUB разность 5.96E-8 отмечается лишь на некоторых углах (+/- 50, 245, 255 … - проверял через 5 градусов)). Извиняюсь за грубое редактирование (только для примера).

---

Вас не смущает, что добавив цикл вы только усложнили функцию циклом, ничего не улучшив,
точность не лучше, на 90 и 270 градусах не работает, длительность цикла ПР выросла (по этой причине и отказался от рядов Тейлора). Какой смысл тогда?
Моя доработка хотя бы исправила работу на 90 и 270 градусах.

kondor3000.
Я действительно поспешил и ошибся в тестировании. Похоже в SIN_Rad ошибка спрятана в выражении:

b:=(Pi-a) * (bool_to_real(a>1.5707964))+a * bool_to_real(not(a>=1.5707964)) ; // 1,5707964

Именно оно обращается в ноль при 1,5707964.
Ну а пятикратный вызов простеньких выражений оформленный в виде цикла скорее всего не тяжелее пятикратного вызова pow()
( pow() во все времена была более тяжелой и часто менее точной чем арифметические операции, но возможно в Logic и не так).

kondor3000.
1 Похоже для устранения ошибки в SIN_Rad в выражении для b надо заменить сравнение «больше или равно»
на просто «больше», т.е. как то так:

b:=(Pi-a) * (bool_to_real(a>1.5707964))+a * bool_to_real(not(a>1.5707964)) ; // 1,5707964

2 В исправленной функции (из-за того что сравнение для исключения проводилось с входным аргументом
а не преобразованным значением) были исключены только четыре точки. А их гораздо больше, например,
минус 4,712389, минус 7,853982 и т.п. Можно конечно оговорить область определения, но …

Вы автор – Вам виднее.

Цитата:

---

Сообщение от zjWlad

kondor3000.
1 Похоже для устранения ошибки в SIN_Rad в выражении для b надо заменить сравнение «больше или равно»
на просто «больше», т.е. как то так:

b:=(Pi-a) * (bool_to_real(a>1.5707964))+a * bool_to_real(not(a>1.5707964)) ; // 1,5707964

2 В исправленной функции (из-за того что сравнение для исключения проводилось с входным аргументом
а не преобразованным значением) были исключены только четыре точки. А их гораздо больше, например,
минус 4,712389, минус 7,853982 и т.п. Можно конечно оговорить область определения, но …

Вы автор – Вам виднее.

---

А вот это хорошая мысль, достаточно убрать = и оставить > и всё работает правильно во всех точках.
И я не автор, всего лишь переписал на ST функции и астротаймер с форума.