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 1211 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

LEGO by MEgg
[December 11, 2017, 02:59:41 PM]


Your workshop space by viswesh
[December 07, 2017, 09:32:47 PM]


printed CNC engraver by tinhead
[December 05, 2017, 05:09:44 PM]


Albert Michelson's Harmonic Analyzer by jinx
[December 05, 2017, 01:05:48 PM]


Raspberry Pi / Python eBook - free for the next 11 hours! by MEgg
[December 03, 2017, 04:20:50 PM]


Fun Fact by Impala
[December 03, 2017, 11:26:45 AM]


CZ-1 2.0: problems, calibration and printing, part2 by MEgg
[December 02, 2017, 12:15:34 PM]


Commissioning a Nano Leg Module by BaldwinK
[November 28, 2017, 11:42:31 AM]


T-HR3 by jinx
[November 21, 2017, 03:41:11 PM]


metal chassis desktop bot by jinx
[November 21, 2017, 01:28:34 PM]


MKS Sbase by jinx
[November 20, 2017, 03:53:21 PM]


Skull Line Follower by Ittiporng
[November 20, 2017, 03:23:18 AM]


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]

* Recent Posts

Re: Autodownload items (was: Re: LEGO) by MEgg
[December 11, 2017, 02:59:41 PM]


Re: Your workshop space by viswesh
[December 07, 2017, 09:32:47 PM]


Re: Your workshop space by tinhead
[December 07, 2017, 08:50:59 PM]


Re: Your workshop space by MEgg
[December 07, 2017, 03:13:51 PM]


Re: Your workshop space by mogul
[December 07, 2017, 02:43:22 PM]


Re: printed CNC engraver by tinhead
[December 05, 2017, 05:09:44 PM]


printed CNC engraver by jinx
[December 05, 2017, 04:59:18 PM]


Albert Michelson's Harmonic Analyzer by jinx
[December 05, 2017, 01:05:48 PM]


Re: Your workshop space by craighissett
[December 04, 2017, 09:32:32 AM]


Re: Your workshop space by tinhead
[December 04, 2017, 09:26:32 AM]


Re: Raspberry Pi / Python eBook - free for the next 11 hours! by MEgg
[December 03, 2017, 04:20:50 PM]


Re: Fun Fact by Impala
[December 03, 2017, 11:26:45 AM]


Changed calibration after long time by MEgg
[December 02, 2017, 12:15:34 PM]


Commissioning a Nano Leg Module by BaldwinK
[November 28, 2017, 11:42:31 AM]


T-HR3 by jinx
[November 21, 2017, 03:41:11 PM]