Отключение плавучки и KGP Опции

[Ash_snz]

Господа. имеется тема для размышления.
Есть процессор мипсовый с сопроцессором для плавучки. 32-х разрядные вычисления делаются на отлично.
Цель попробовать отказаться от 32-х разрядов и все плавучие вычисления реализовать на 16 разрядах, соответственно не используя дополнительное ядро.

Я так понимаю один из вариантов -использование фиксированной запятой. Сам я с этим не сталкивался, кто знаком подскажите где можно изучить вопрос?

Второй вариант - перекинуть это на сильные плечи компилятора. Пробовал с помощью KGP и библиотек сбилденных на soft-float реализовать ту же программу. Компилятор ругается на многие андефайны, относящиеся к плавучке, хотя в libm их видно. Как итог - не собирается.

Кто, что может подсказать по этому вопросу?

[Ash_snz]

math.h на месте.
вот так ругается:

Код

(.text+OxleO): undefined reference to "_floatsisf'
main.о:  In function 'rs_transmit_digit':
(.text+Oxlf8):  undefined reference to "_floatsisf'
main, о:  In function    'rs_transmit_digit' :
(.text+0x204):  undefined reference to '_mulsf3'
main.о:  In function 'rs_transmit_digit1:
(.text+0x224):  undefined reference to '_fixsfsi'
main, о :  In function    rs_transmit_digit' :
(.text+0x808):  undefined reference to "_subsf31
main.o:  In function 'cycle_2_datchika':
(.text+0xac8):  undefined reference to "_addsf31
main.o:  In function "cycle_2_datchika1:
(.text+0xb2c): undefined reference to "_extendsfdf2'
main.o:  In function ''cycle_2_datchika' :
(.text+0xb58): undefined reference to "_divdf31
main.o:  In function 'cycle_2_datchika1:
(.text+Oxcdc): undefined reference to "_truncdfsf2'
main.o:  In function "cycle_2_datchika1:
( . text+0xce4) : undefined reference to   _truncdf sf 2 1
main.o:  In function vcycle_2_datchika1:
(.text+0xe48):  undefined reference to "_subdf3'
main.o:   In function "cycle_2_datchika1:
(.text+0xe68):  undefined reference to "_truncdfsf21
с:\kgp_mips_elf\mips-kgp-elf\lib\libm.a(lib_a-w_sqrt.o): In function "sqrf: (.text.sgrt+0x54):  undefined reference to -_ltdf21
с:\kgp_mips_elf\mips-kgp-elf\lib\libm.a(lib_a-w_sqrt.o): In function "sqrt": (.text.sqrt+0xa8):  undefined reference to "_divdf3'
c:\kgp_mips_elf\mips-kgp-elf\lib\libm.adib_a-e_sqrt.o):  In function "_ieee754__sqrt' :
(.text._ieee754_sqrt+0x34):  undefined reference to "_muldf3'
с:\kgp_mips_elf\mips-kgp-elf\lib\libm.a(lib_a-e_sqrt.o):  In function _ieee754_sqrt1 :
(.text._ieee754_sqrt+0x48):  undefined reference to "_adddf3'
с:\kgp_mips_elf\mips-kgp-elf\lib\libm.a(lib_a-e_sqrt.o):  In function "_ieee754_sqrt1:
(.text._ieee754_sqrt+0x84): undefined reference to *_subdf3'
с:\kgp_mips_elf\mips-kgp-elf\lib\libm.a(lib_a-e_sqrt.o):  In function "_ieee754_sqrt1:
(.text._ieee754_sqrt+0x98): undefined reference to "_divdf31
cs-make:  ***   [main.elf]   Error 1

Прям и не знаю что делать...

я думаю решим вопрос, только зачем надевать штаны через голову если есть FPU, сдается что Вы не пытаетесь решить задачу а пытаетесь решить задачу так как Вам хочется.

рас уж поднялась тема про мипсы то можно по ходу и преесобрать все. какое ядро?

Ядро многострадального R3000 (комдив, вестимо ). Дело даже не в том - как решить задачу. Это попытка научиться и попробовать работать без плавучки, а также выяснить разницу в производительности.
Вот только не понимаю, что можно пересобрать? у нас есть сборка под наш проц и lib.m как с плавучкой так и без, с плавучкой работает, а вот без - нет.

Чувствую где-то пробел в моих знаниях

[alx2]

math.h на месте.
вот так ругается:

Код

(.text+OxleO): undefined reference to "_floatsisf'
main.о:  In function 'rs_transmit_digit':

Так это у Вас вовсе не компилятор ругается, а линкер!
Следовательно, разбираться надо не с инклудами, а с библиотеками. Судя по именам символов, у Вас не подключена библиотека libgcc.

[Ash_snz]

Так это у Вас вовсе не компилятор ругается, а линкер!
Следовательно, разбираться надо не с инклудами, а с библиотеками. Судя по именам символов, у Вас не подключена библиотека libgcc.

блин, как догадался? пойду проверю.

и почему я так не могу догадаться?!

[alx2]

блин, как догадался? пойду проверю.

Вы будете смеяться. Я тупо погрепал имя символа по своим библиотекам. Сначала тут:

Код

alx% grep -rl _floatsisf /usr/local/arm-elf/lib
/usr/local/arm-elf/lib/thumb/libm.a
/usr/local/arm-elf/lib/libm.a

Но здесь этот символ оказался Undefined:

Код

alx% arm-elf-nm /usr/local/arm-elf/lib/libm.a|grep _floatsisf
         U __floatsisf
         U __floatsisf
         U __floatsisf
         U __floatsisf
         U __floatsisf
         U __floatsisf
         U __floatsisf
         U __floatsisf
         U __floatsisf
         U __floatsisf
         U __floatsisf
         U __floatsisf
         U __floatsisf
         U __floatsisf
         U __floatsisf
         U __floatsisf
         U __floatsisf

А вот тут он нашелся:

Код

alx% grep -rl _floatsisf /usr/local/lib/gcc/arm-elf
/usr/local/lib/gcc/arm-elf/4.5.2/plugin/include/insn-flags.h
/usr/local/lib/gcc/arm-elf/4.5.2/libgcc.a
/usr/local/lib/gcc/arm-elf/4.5.2/thumb/libgcc.a
/usr/local/lib/gcc/arm-elf/4.5.2/fpu/libgcc.a
alx% arm-elf-nm /usr/local/lib/gcc/arm-elf/4.5.2/libgcc.a |grep _floatsisf
000001d0 T __floatsisf

[Ash_snz]

Вы будете смеяться. Я тупо погрепал имя символа по своим библиотекам.

Благодарствую за подсказку!
В итоге добавил libgcc и libg и все собралось. Скорость выполнения уменьшилась в 50 раз, точность (на удивление) осталась на том же уровне! /тестим дальше/

Остался последний и главный вопрос.
Существуют ли способы реализации плавучки, а точнее даже дробных вычислений не с 32-мя разрядами, а с 16-ю?
Понимаю, что для перехода придется изменять программу. вопрос в том как? Куда двигаться, кто знает?
может кто реализовывал вычисления на 16-ти разрядных системах без сопроцессоров?

Сообщение отредактировал Ash_snz - Nov 13 2012, 16:18

[maksimp]

Скорость выполнения уменьшилась в 50 раз,

Ядро ARM тратит около 50 тактов на сложение или умножение с 32 битной плавающей точкой при отсутствии аппаратной поддержки. То есть у вас аппаратно операции за 1 такт выполняются.

точность (на удивление) осталась на том же уровне!

Если реализация аппаратного и программного вычисления соответствуют стандарту, то операциии + - * / дают чётко определённый результат. Поэтому так и должно быть.

Существуют ли способы реализации плавучки, а точнее даже дробных вычислений не с 32-мя разрядами, а с 16-ю?

Есть такое. Ищите в сети "half precision floating point".
Есть даже процессоры с аппаратной поддержкой такого.

[Ash_snz]

Ядро ARM тратит около 50 тактов на сложение или умножение с 32 битной плавающей точкой при отсутствии аппаратной поддержки. То есть у вас аппаратно операции за 1 такт выполняются.

Наше ядро тратит до 16 тактов на деление, на умножение до 8 кажется тактов.

Если реализация аппаратного и программного вычисления соответствуют стандарту, то операциии + - * / дают чётко определённый результат. Поэтому так и должно быть.

Есть такое. Ищите в сети "half precision floating point".
Есть даже процессоры с аппаратной поддержкой такого.

Спасибо, уже ищу.

скорее увеличилась .

С точностью смотрите окуратно: у меня количество знаков после запятой было критично и вначале точность поплыла, увеличивал количество знаков после запятой, чтоб проходить по точности.

Как меняли количество знаков? что за методы?

[cant]

Как меняли количество знаков? что за методы?

В общем ничего особенного: вначале из лонг два байта было на дробь, т.е точность 1/65536. потом добавил пол байта за счет уменьшения целой части - дипозон чисел позволял.

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