Ошибка в переводе десятичных градусов в градусы минуты секунды

В Лиспе операции остатка от деления нету?

(setq d (round n))

Скорее всего это округление к ближайшему целому (0.6 -> 1), а надо с отбрасыванием дробной части (0.6 -> 0).

В общем рабочий код на Обероне:

MODULE A;
	IMPORT Log := StdLog, In := i21sysIn;
	
	PROCEDURE DegToDms*;
		VAR n: REAL; nInt: INTEGER; d, m, s: INTEGER;
	BEGIN
		In.Open; ASSERT(In.done);
		In.Real(n); ASSERT(In.done);
		nInt := SHORT(ENTIER(n*3600 + 0.5));
		d := nInt DIV 3600;
		m := nInt DIV 60 MOD 60;
		s := nInt MOD 60;
		Log.Int(d); Log.String("° "); 
		Log.Int(m); Log.String("' "); 
		Log.Int(s); Log.String("''"); Log.Ln;
	END DegToDms;
	
END A.

A.DegToDms 0.664631042351755
0° 39' 53''
A.DegToDms 1.3349916626548182
1° 20' 6''
A.DegToDms 0.13333333333333333
0° 8' 0''

Google подсказывает что операция остатка от деления в Лиспе называется mod.

Google подсказывает что операция остатка от деления в Лиспе называется mod.

Да, вот и книжка говорит то же самое

mod dividend divisor [Function]

This function returns the value of dividend modulo divisor; in other words, the remainder after division of dividend by divisor, but with the same sign as divisor.

а вот как переписать код с этой функцией у меня не получается.

А я не понимаю, где ты видишь ошибку?

(format «%d° %d' %d"» d m s))

Шаблон %d отсекает дробную часть.

(format "%d" 0.123123123123) -> "0"

А, я понял. Тебя ноль в секундах смущает. Вот тут ты не отбросил дробную часть.

(setq m (truncate (* (- n d) 60)))

И имей в виду, что %d, как я уже написал, не округлит секунды, а отбросит дробную часть.

(defun decimal-deg-to-dms (n)
  (setq d (round n))
  (setq m (* (- n d) 60))
  (setq s (- (* 3600 (- n d)) (* 60 m)))
  (format "%d° %d' %d''" d m s))

Зачем эта императивщина? Для кого let* придумали? Ты же уже здесь давно темы про лиспо строчишь, неужели этого еще не усвоил?

Ты же уже здесь давно темы про лиспо строчишь, неужели этого еще не усвоил?

Если нет, то что?

только

Я тебе про %d два раза написал. Включи мозги. И почему floor? Тебе отбрасывать надо. Делай truncate.

Включил

(defun decimal-deg-to-dms (n)
  (setq d (round (floor n)))
  (setq m (round (floor (* (- n d) 60))))
  (setq s (round (- (* 3600 (- n d)) (* 60 m))))
  (format "%d° %d' %d''" d m s))

(decimal-deg-to-dms 0.664631042351755) => "0° 39' 53''"
(decimal-deg-to-dms 1.3349916626548182) => 1° 20' 6''"
(decimal-deg-to-dms 0.13333333333333333) => "0° 8' 0''"

С другой стороны ерунда какая-то, я тебя игнорирую и уже давно, а ты опять у меня в треде. А я общаюсь с тобой.

С другой стороны ерунда какая-то, я тебя игнорирую и уже давно, а ты опять у меня в треде. А я общаюсь с тобой.

И правильно! А то оскорблю еще религиозные чувства ненароком. :)

А то оскорблю еще религиозные чувства ненароком. :)

Не совсем понял - религиозные чувства, о чём это.

Не совсем понял - религиозные чувства, о чём это.

Ну вот, к примеру. Начинаю оскорблять. :)

(setq d (round (floor n)))
(setq m (round (floor (* (- n d) 60))))
(setq s (round (- (* 3600 (- n d)) (* 60 m))))

(round (floor ())) - это вообще что-то с чем-то. :) Даже если допустить, то попробуй отрицательные углы ввести. Они же допустимы?

(decimal-deg-to-dms 90.664631042351755) -> 90° 39' 53''
(decimal-deg-to-dms -90.664631042351755) -> -91° 20' 7''

Йоу, оскорбил. Говорю же, бери truncate хотя бы (тогда минус будет у всех компонентов DMS). Или знак обрабатывай.

Об этом вам уже не первую тему об этом говорят. setq без нужды – плохой код.

Этот код использую для максимального угла в 4 градуса, другие углы 8'0" , 6'0". Про отрицательные углы даже не задумывался, вот переделанный с truncate.

(defun decimal-deg-to-dms (n)
  (setq d (round (truncate n)))
  (setq m (round (truncate (* (- n d) 60))))
  (setq s (round (- (* 3600 (- n d)) (* 60 m))))
  (format "%d° %d' %d''" d m s))

(decimal-deg-to-dms -90.664631042351755)   "-90° -39' -53''"

Об этом вам уже не первую тему об этом говорят. setq без нужды – плохой код.

Не знаю что ответить.

(decimal-deg-to-dms -90.664631042351755)   "-90° -39' -53''"

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

Покомпонентно складывать углы, вычитать углы (с учетом переполнения, конечно же).

А вот так далеко я даже не думал. Тогда надо widget и оформлять всё как утилу для счёта углов. Наверное так с таким замахом.

(round (truncate

round никакой не нужен. После truncate уже integer.

(truncate ARG &optional DIVISOR)

Truncate a floating point number to an int.
Rounds ARG toward zero.
With optional DIVISOR, truncate ARG/DIVISOR.

(decimal-deg-to-dms -90.664631042351755) "-90° -39' -53""

Только всё равно где-то косячок

(defun dms-to-decimal-deg (d m s)
  (+ d
     (/ m 60.0)
     (/ s 3600.0)))

(dms-to-decimal-deg -90 -39 -53)    -90.66472222222222

Потому что ты (оскорбляю чувства) округлял секунды и восстанавливаешь угол из округленного значения фактически. Чего ты ожидаешь?

Тупо переписал

[code] (defun decimal-deg-to-dms (n) (setq d (truncate n)) (setq m (truncate (* (- n d) 60))) (setq s (truncate (- (* 3600 (- n d)) (* 60 m)))) (format «%d° %d’ %d''» d m s))

(decimal-deg-to-dms -90.664631042351755) «-90° -39’ -52''»

должно же быть «-90° -39’ -53''» [/code]

мне надо сидеть с мануалом и вспоминать, что там да как.

должно же быть «-90° -39’ -53"

round у секунд оставь. round перед truncate был лишним, а не этот.

-90.664631042351755 не равно -90.66472222222222

Господи. :( Ошибка в переводе десятичных градусов в градусы минуты секунды (комментарий)

тут округляется тоже

(dms-to-decimal-deg 0 8 1) 0.1336111111111111

(decimal-deg-to-dms 0.1336111111111111) «0° 8’ 1''»

Ну и [code] (decimal-deg-to-dms -90.664631042351755) «-90° -39’ -53''»

(decimal-deg-to-dms -90.66472222222222) «-90° -39’ -53''» [/code]

(decimal-deg-to-dms 0.1336111111111111) «0° 8’ 1"»

У тебя фундаментальное непонимание. Я не знаю, что с этим делать. Все здесь правильно. Число очень хорошо раскладывается. Округление секунд тут на уровне самых младших значащих разрядов мантиссы получается. Число такое везучее. В примере с "-90° -39’ -53"" совсем все не так.

https://en.wikipedia.org/wiki/Double-precision_floating-point_format

UPD Сделай так: (decimal-deg-to-dms 0.1336111111111112) -> «0° 8’ 1"» и вся идиллия развалится.

Число такое везучее.

нет, просто сначала

(dms-to-decimal-deg 0 8 1) 0.1336111111111111

(decimal-deg-to-dms 0.1336111111111111) «0° 8’ 1''»

В примере с "-90° -39’ -53"" совсем все не так.

Знаю, просто забуксовал.

нет, просто сначала

Не нет, а да. Число везучее, так как округления практически не получается. Единица у тебя получается из оставшегося 0.9999999999999600.

В общем, так. Если хочешь поточнее, то можешь оставлять секунды неокругленными.

https://en.wikipedia.org/wiki/Decimal_degrees#Example

Note that with this formula only D can be negative and only S may have a fractional value.

Ну тебя прёт. Конечно программирование требует таких знаний.

Вот приведу аналогию. По ГОСТ 3325-85 рассчитал допуск соосности вала, для подшипника нормального класса точности, получил 0.0056 мм, длина вала 1000 мм. Показал знакомому технологу, он сказал что это теоретическая цифра в реалиях не достижимая на обычных заводах, а вот про практическую цифру он промолчал, хотя все обычные заводы делают валы. У всех размеров подшипников нормального класса точности допустимый перекос колец 8 минут. И допуск соосности вала будет лежать в тысячных.