collapse

Welcome!



Welcome to Robot Rebels, an online robot maker community.

Register today to post your projects, ask questions, share knowledge and meet like-minded people from around the world.


The RobotRebel.org Community

Author Topic: CORDIC algorithms  (Read 1127 times)

deshipu

  • Member
  • *
  • Posts: 113
  • Radomir Dopieralski
    • Homepage
CORDIC algorithms
« on: December 21, 2016, 04:46:58 PM »
We often need to use trigonometric and inverse trigonometric functions in our microcontroller projects. Whether it's calculating inverse kinematics, figuring out angles and distances, stabilizing a quadcopter or otherwise dealing with angles -- those functions are indispensable. Personally I'm using the atan2 and acos functions in my robots for calculating inverse kinematics of the legs. I'm using the standard floating-point functions from math.h, but that is suboptimal.

First of all, the implementations of those functions are rather large. Just adding them to my AVR code makes it swell by some 3000 bytes. That is partially because they probably have some lookup tables and complicated logic, but also because they are floating-point functions, and pull in all the floating point multiplication, division and other operations, which are quite costly on a poor microcontroller without a co-processor.

Second, they take and return values that are scaled rather arbitrarily (-1:1, -PI:PI, etc.), so afterwards I have to do some more floating point operations to get them to the 0-1800 range that I need for my servos.

So I decided to try and find integer replacements for those functions. I couldn't find any that would be suitable. They were either floating-point, using some strange fixed-point types, using huge lookup tables, etc. However, during my search I discovered a family of algorithms called CORDIC, which allows me to calculate trigonometric and inverse trigonometric functions with minimal code, using only a lookup into a tiny table, additions and bit shifts. If you are interested about how those algorithms work, I wrote a short article about it at https://hackaday.io/project/6050-tote/log/50664-cordic

Anyways, I still couldn't find a decent implementation of those algorithms for the AVR, so I set out to write my own. I started by prototyping them in Python (because it's easy to debug for me and fast to write), and when I got them to work with a reasonable accuracy (I don't need much, 0.1° is enough, my servos won't move at better accuracy anyways), I re-wrote them in AVR C. The code is at https://bitbucket.org/thesheep/imath/src/tip/imath.c and it currently takes about 514 bytes when compiled with my avr-gcc. I suspect a couple of bytes could be shaved here and there from it, but I'm pretty happy with it already. Oh, I also added a square root function in there, for completeness.

Enjoy!

Protowrxs

  • Member
  • *
  • Posts: 139
    • The Protowrxs World
Re: CORDIC algorithms
« Reply #1 on: December 21, 2016, 10:27:17 PM »
Although I really don't understand it all it has always been a pain trying to do many math functions on an Arduino even for dead reckoning, etc for me. I understand the concept of IK but not the method. I followed up on some of your links and it does help me some, just have much more to read to try to fully understand.

Thanks for sharing your info. Can this stuff be added as an Arduino library in any manner?
(Sorry, more stuff I do not understand)

I'd really like to try again to build a four legged walker but do it right this time. Maybe someday.

Stephen
The past is done, the future is not here, live today...
My Bots: http://www.Protowrxs.com/Robotics

deshipu

  • Member
  • *
  • Posts: 113
  • Radomir Dopieralski
    • Homepage
Re: CORDIC algorithms
« Reply #2 on: December 22, 2016, 04:40:14 AM »
Can this stuff be added as an Arduino library in any manner?

Not sure what you mean. If you are asking if you can use this in your Arduino programs, then yes, the easiest way is to just copy the imath.h and imath.c files into the directory of your sketch, and do #include "imath.h" in your code. If you are asking whether I plan to make a pull request to the Arduino core to have this code included by default, then the answer is no.

 

* Search


* Recent Topics

Just scary to me... by jinx
[November 16, 2017, 03:05:30 AM]


"Artie" version 3.0 (RTV3) by jinx
[November 14, 2017, 03:11:04 AM]


Anycubic Kossel (pulley) by Bajdi
[November 13, 2017, 04:48:04 PM]


grbl-LPC by jinx
[November 12, 2017, 02:48:46 PM]


mks ftf2.4 simple bezel by jinx
[November 11, 2017, 04:18:13 PM]


Controllers by jinx
[November 11, 2017, 06:48:21 AM]


Magnetometer vs. Electric Motor by ZeroMax
[November 08, 2017, 02:22:47 PM]


Raspberry Pi / Python eBook - free for the next 11 hours! by MEgg
[November 02, 2017, 06:10:53 PM]


who's using atom by jinx
[November 02, 2017, 03:49:38 AM]


3d review of a delta by jinx
[October 27, 2017, 11:30:38 PM]


Coding a Nano i2c Slave by BaldwinK
[October 25, 2017, 05:40:27 AM]


Drawdio by mogul
[October 24, 2017, 02:35:00 PM]


My scope focuser upgrade by jscottb
[October 22, 2017, 03:13:23 PM]


thermoplastic bed surface by Bajdi
[October 22, 2017, 03:12:37 PM]


Kossel Build by Deity
[October 20, 2017, 11:23:27 AM]

* Recent Posts

Re: Just scary to me... by jinx
[November 16, 2017, 03:05:30 AM]


Just scary to me... by Protowrxs
[November 15, 2017, 05:31:10 PM]


Re: "Artie" version 3.0 (RTV3) by jinx
[November 14, 2017, 03:11:04 AM]


Re: Anycubic Kossel (pulley) by Bajdi
[November 13, 2017, 04:48:04 PM]


Re: Anycubic Kossel (pulley) by Deity
[November 13, 2017, 04:08:05 PM]


Re: "Artie" version 3.0 (RTV3) by lukeyes2
[November 12, 2017, 06:29:44 PM]


grbl-LPC by jinx
[November 12, 2017, 02:48:46 PM]


mks ftf2.4 simple bezel by jinx
[November 11, 2017, 04:18:13 PM]


Controllers by jinx
[November 11, 2017, 06:48:21 AM]


Magnetometer vs. Electric Motor by ZeroMax
[November 08, 2017, 02:22:47 PM]


Re: Raspberry Pi / Python eBook - free for the next 11 hours! by MEgg
[November 02, 2017, 06:10:53 PM]


Re: "Artie" version 3.0 (RTV3) by lukeyes2
[November 02, 2017, 10:47:29 AM]


Re: who's using atom by jinx
[November 02, 2017, 03:49:38 AM]


Re: "Artie" version 3.0 (RTV3) by Bajdi
[November 01, 2017, 04:11:20 PM]


Re: who's using atom by Bajdi
[November 01, 2017, 04:08:47 PM]