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 668 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: 138
    • 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

Encoders Really? by DangerousThing
[Today at 05:37:59 AM]


Blockchains, Ethereum, and a Safer SkyNet by jinx
[June 26, 2017, 01:17:32 PM]


Z motor mount by jinx
[June 26, 2017, 04:05:25 AM]


Spinny Flowers by Impala
[June 25, 2017, 06:24:33 PM]


Automobile drag coefficient by Impala
[June 23, 2017, 01:18:46 AM]


Arduino fork by ossipee
[June 20, 2017, 03:20:37 PM]


ghetto compressor revisited by 1 what
[June 19, 2017, 04:07:09 AM]


Don Quixote Re-loaded by Gareth
[June 18, 2017, 07:10:39 AM]


Spammers by MEgg
[June 16, 2017, 11:11:26 AM]


5 min meeting timer by mogul
[June 16, 2017, 04:55:08 AM]


Picture Topic by Impala
[June 14, 2017, 05:17:43 PM]


Bluetooth Low Energy in JavaScript by Ladvien
[June 10, 2017, 05:10:21 PM]


Gait Creation Using Excel by BaldwinK
[June 08, 2017, 01:41:11 PM]


Which micro motor solution would be the strongest ? by programer786
[June 07, 2017, 02:50:31 PM]


Fun Fact by jinx
[June 07, 2017, 03:08:01 AM]

* Recent Posts

Re: Encoders Really? by DangerousThing
[Today at 05:37:59 AM]


Re: Blockchains, Ethereum, and a Safer SkyNet by jinx
[June 26, 2017, 01:17:32 PM]


Z motor mount by jinx
[June 26, 2017, 04:05:25 AM]


Spinny Flowers by Impala
[June 25, 2017, 06:24:33 PM]


Re: Blockchains, Ethereum, and a Safer SkyNet by mogul
[June 23, 2017, 01:53:15 AM]


Re: Automobile drag coefficient by Impala
[June 23, 2017, 01:18:46 AM]


Blockchains, Ethereum, and a Safer SkyNet by mtriplett
[June 22, 2017, 11:47:08 AM]


Re: Automobile drag coefficient by Impala
[June 21, 2017, 09:45:12 PM]


Re: Arduino fork by ossipee
[June 20, 2017, 03:20:37 PM]


Re: Arduino fork by mogul
[June 20, 2017, 02:28:32 PM]


Arduino fork by ossipee
[June 20, 2017, 01:35:33 PM]


Re: ghetto compressor revisited by 1 what
[June 19, 2017, 04:07:09 AM]


Re: Don Quixote Re-loaded by Gareth
[June 18, 2017, 07:10:39 AM]


Re: Spammers by MEgg
[June 16, 2017, 11:11:26 AM]


Re: 5 min meeting timer by mogul
[June 16, 2017, 04:55:08 AM]