Skip to main content

Блог инженера

Блог о минимализме, инжинерии и программировании.



Целочисленная тригонометрия для Форт.

  | #Форт#программирование

Делаем фукнцию для работы с тригонометрией. Имеем простейшую тригонометрическую библиотеку.

\ trig 99.12.10 11:46 am NAB
: table create does> swap 2* + @ ;
table tsin
0 , 175 , 349 , 523 ,
698 , 872 , 1045 , 1219 ,
1392 , 1564 , 1736 , 1908 ,
2079 , 2250 , 2419 , 2588 ,
2756 , 2924 , 3090 , 3256 ,
3420 , 3584 , 3746 , 3907 ,
4067 , 4226 , 4384 , 4540 ,
4695 , 4848 , 5000 , 5150 ,
5299 , 5446 , 5592 , 5736 ,
5878 , 6018 , 6157 , 6293 ,
6428 , 6561 , 6691 , 6820 ,
6947 , 7071 , 7193 , 7314 ,
7431 , 7547 , 7660 , 7771 ,
7880 , 7986 , 8090 , 8192 ,
8290 , 8387 , 8480 , 8572 ,
8660 , 8746 , 8829 , 8910 ,
8988 , 9063 , 9135 , 9205 ,
9272 , 9336 , 9397 , 9455 ,
9511 , 9563 , 9613 , 9659 ,
9703 , 9744 , 9781 , 9816 ,
9848 , 9877 , 9903 , 9925 ,
9945 , 9962 , 9976 , 9986 ,
9994 , 9998 , 10000 ,
: (sin) ( n -- n')
  dup 90 > if  180 swap -  then
  tsin ;
: sin ( n -- n')
  360 /mod drop dup 0 <
  if  360 +  then  dup 180 >
  if  180 - (sin) negate  else  (sin)  then ;
: cos ( n -- n')  360 /mod drop 90 + sin ;

С ней есть проблема, так делать не следует. Есть слово CREATE, его и следует использовать. Определяем таблицу со значениями синусов для всех целых углов от 0 до 90 градусов.

CREATE tsin  0000 , 0175 , 0349 , 0523 , 0698 , 0872 , 1045 , 1219 , 1392 , 1564 , 1736 , 1908 , 2079 , 2250 , 2419 , 2588 , 2756 , 2924 , 3090 , 3256 , 3420 , 3584 , 3746 , 3907 , 4067 , 4226 , 4384 , 4540 , 4695 , 4848 , 5000 , 5150 , 5299 , 5446 , 5592 , 5736 , 5878 , 6018 , 6157 , 6293 , 6428 , 6561 , 6691 , 6820 , 6947 , 7071 , 7193 , 7314 , 7431 , 7547 , 7660 , 7771 , 7880 , 7986 , 8090 , 8192 , 8290 , 8387 , 8480 , 8572 , 8660 , 8746 , 8829 , 8910 , 8988 , 9063 , 9135 , 9205 , 9272 , 9336 , 9397 , 9455 , 9511 , 9563 , 9613 , 9659 , 9703 , 9744 , 9781 , 9816 , 9848 , 9877 , 9903 , 9925 , 9945 , 9962 , 9976 , 9986 , 9994 , 9998 , 10000 ,

: raw-sin tsin swap cells + @ ;

: (sin) ( n -- n')
  dup 90 > if  180 swap -  then
  raw-sin ;

: sin ( n -- n')
  360 /mod drop dup 0 <
  if  360 +  then  dup 180 >
  if  180 - (sin) negate  else  (sin)  then ;

: cos ( n -- n')  360 /mod drop 90 + sin ;

С create does> было чуть красивее, можно было чуть изменить функцию для случая, когда на стеке хранятся не 16-разрядные числа. Но так, как сделал я - чуть правильней. Данные и функции хранятся раздельно. Возможно, это чуть менее по-фортерски. Пока для меня достаточно, что этот код работает.

About Mikhail Kiselev

Photo of Mikhail Kiselev

Приветствую в моём блоге! 😄 Меня зовут Михаил. Я инженер и программист. Живу в Израиле. Но мой блог связан с работой в Сибири и на Сахалине, путешествую где придётся. Я предпочитаю пост в блог посту в твиттер. Описание полезной технологии или гаджета предпочитаю описанию заката или посиделок в кафе.