i`m doing my thesis comparing CORDIC with polynomial in counting arctan with fixed point. I`m using Q15 format now. I`m using this site CORDIC arctan as a referenced when making with floating point. The problem there`s a lot of error when i try to make it with fixed point.
this is my program -
#include "Unit1.h"
-
#include "math.h"
-
#include "fixed_math.hpp"
-
#define MAXBITS 15
-
static float invGain1;
-
static float atanTable[MAXBITS];
-
static float gain1Cordic();
-
-
-
void initCordic()
-
{
-
/* must call this first to initialise the constants.
-
* of course, here i use the maths library, but the
-
* values would be precomputed.
-
*/
-
float t = 1.0f;
-
int i;
-
for (i = 0; i < MAXBITS; ++i)
-
{
-
atanTable[i] = atan(t);
-
t /= 2;
-
}
-
-
/* set constants */
-
invGain1 = 1/gain1Cordic();
-
}
-
-
/* CORDIC m=1, y-->0 */
-
-
Q15 cordic(Q15 &x0, Q15 &y0, Q15 &z0, Q15 vecmode)
-
{
-
-
short t;
-
Q15 x, y, z;
-
int i;
-
-
t = 1.0f;
-
x = x0; y = y0; z = z0;
-
-
for (i = 0; i < MAXBITS; ++i)
-
{
-
-
double x1;
-
-
if (vecmode >= 0.0f && y < vecmode || vecmode<0.0f && z >= 0.0f)
-
{
-
Q15 x1 = x - y>>t;
-
y = y + x>>t;
-
z = z - atanTable[i];
-
}
-
else
-
{
-
Q15 x1 = x + y>>t;
-
y = y - x>>t;
-
z = z + atanTable[i];
-
}
-
-
x = x1;
-
t /= 2;
-
}
-
-
x0 = x;
-
y0 = y;
-
z0 = z;
-
-
}
-
-
-
static float gain1Cordic()
-
{
-
/* compute gain by evaluating cos(0) without inv gain */
-
float x, y, z;
-
-
x = 1;
-
y = 0;
-
z = 0;
-
cordic(&x, &y, &z, 1.0f);
-
return x;
-
}
-
-
Q15 atanCordic(Q15 a)
-
{
-
/* domain: all a */
-
Q15 x = 1.0f;
-
Q15 z = 0.0f;
-
cordic(&x, &a, &z, 0.0f);
-
return z;
-
}
-
-
void __fastcall TForm1::Button1Click(TObject *Sender)
-
{
-
AnsiString buf;
-
float x;
-
float v,c;
-
-
-
initCordic();
-
-
for (x=0.0;x<=1.0;x=x+0.1)
-
{
-
-
Q15 y = atanCordic(x);
-
v = atan(x);
-
-
-
Memo1->Lines->Add(buf.sprintf("%5.2f\n", x));
-
Memo2->Lines->Add(buf.sprintf("%10.6f\n",Q15ToShort(v)/32768.0));
-
Memo3->Lines->Add(buf.sprintf("%10.6f\n",Q15ToShort(y)/32768.0));
-
-
}
-
-
-
}
-
-
anyone can help?
if you have suggestion about making it with others rather than Q15 format, as long it is fixed point, it is very usefull to
12 9600 sicarie 4,677
Recognized Expert Moderator Specialist
astri-
I'm not sure if I can help with the logic of your program (that's a little over my head!), but it's more than likely that I or someone else browsing through would be able to help you with the error messages that you are getting, if you could post those.
sicarie
RedSon 5,000
Recognized Expert Expert
Can you give us an example of your input and output. What kind of error are you getting?
i`m bolding the font ONLY for error messages(ignore the warning) in my program -
-
#include "Unit1.h"
-
#include "math.h"
-
#include "fixed_math.hpp"
-
#define MAXBITS 15
-
-
static float invGain1;
-
static float atanTable[MAXBITS];
-
static float gain1Cordic();
-
-
-
void initCordic()
-
{
-
/* must call this first to initialise the constants.
-
* of course, here i use the maths library, but the
-
* values would be precomputed.
-
*/
-
float t = 1.0f;
-
int i;
-
for (i = 0; i < MAXBITS; ++i)
-
{
-
atanTable[i] = atan(t);
-
t /= 2;
-
}
-
-
/* set constants */
-
invGain1 = 1/gain1Cordic();
-
}
-
-
/* CORDIC m=1, y-->0 */
-
-
Q15 cordic(Q15 &x0, Q15 &y0, Q15 &z0, Q15 vecmode)
-
{
-
-
short t;
-
Q15 x, y, z;
-
int i;
-
-
t = 1.0f;
-
x = x0; y = y0; z = z0;
-
-
for (i = 0; i < MAXBITS; ++i)
-
{
-
-
double x1;
-
-
if (vecmode >= 0.0f && y < vecmode || vecmode<0.0f && z >= 0.0f)
-
{
-
Q15 x1 = x - y>>t;
-
y = y + x>>t;
-
z = z - atanTable[i];
-
}
-
else
-
{
-
Q15 x1 = x + y>>t;
-
y = y - x>>t;
-
z = z + atanTable[i];
-
}
-
-
x = x1;
-
t /= 2;
-
}
-
-
x0 = x;
-
y0 = y;
-
z0 = z;
-
-
}
-
-
-
static float gain1Cordic()
-
{
-
/* compute gain by evaluating cos(0) without inv gain */
-
float x, y, z;
-
-
x = 1;
-
y = 0;
-
z = 0;
-
cordic(&x, &y, &z, 1.0f);-->THE MAIN PROBLEM,
-
return x; the error messages repeated so much in this line
-
}
-
-
Q15 atanCordic(Q15 a)
-
{
-
/* domain: all a */
-
Q15 x = 1.0f;
-
Q15 z = 0.0f;
-
cordic(&x, &a, &z, 0.0f);
-
return z;
-
}
-
-
void __fastcall TForm1::Button1Click(TObject *Sender)
-
{
-
AnsiString buf;
-
float x;
-
float v,c;
-
-
-
initCordic();
-
-
for (x=0.0;x<=1.0;x=x+0.1)
-
{
-
-
Q15 y = atanCordic(x);
-
v = atan(x);
-
-
Memo1->Lines->Add(buf.sprintf("%5.2f\n", x));
-
Memo2->Lines->Add(buf.sprintf("%10.6f\n",y));
-
Memo3->Lines->Add(buf.sprintf("%10.6f\n",v));
-
-
-
}
-
-
-
}
-
-
this is my error messages
[code]
[C++ errorー] Unit1.cpp(67): E2015 'Q15:: operator -(const Q15 &) const' and 'Q15:: operator -(const Q30 &) const' classification are ambiguous retrieval
[C++ errorー] Unit1.cpp(73): E2015 'Q15:: operator +(const Q15 &) const' and 'Q15:: operator +(const Q30 &) const' classification are ambigous retrieval
[C++ error ー] Unit1.cpp(76): E2015 'Q15::Q15(const short)' and 'Q15::Q15(const float)' classification are ambigous retrieval
[C++ warning] Unit1.cpp(84): W8070 return function value
[C++ warning] Unit1.cpp(95): W8030 'x0' parameter(cordi c(Q15 &,Q15 &,Q15 &,Q15))is used as temporarily variables
[C++ error] Unit1.cpp(95): E2064 'Q15 &' cant initialized as 'float *'
[C++ error] Unit1.cpp(95): E2342 'x0'parameter is define as Q15 & so it cant be change to float *
[C++ warning] Unit1.cpp(95): W8030 'y0' parameter(cordi c(Q15 &,Q15 &,Q15 &,Q15))is used as temporarily variables
[C++ error] Unit1.cpp(95): E2064 'Q15 &' cant initialized as 'float *'
[C++ error] Unit1.cpp(95): E2342 'y0'b parameter is Q15 & type and cant be change to float *
[C++ warning] Unit1.cpp(95): W8030 'z0' parameters(cord ic(Q15 &,Q15 &,Q15 &,Q15))is used as temporarily variables
[C++ error] Unit1.cpp(95): E2064 'Q15 &' cant initialized as 'float *'
[C++ error] Unit1.cpp(95): E2342 'z0' parameter is Q15 & type and cant be change to float *
[C++ warning] Unit1.cpp(97): W8004 'z' subtituted value is not use
[C++ warning] Unit1.cpp(97): W8004 'y' subtituted value is not use [C++ warning] Unit1.cpp(104): W8030 'x0' parameter(cordi c(Q15 &,Q15 &,Q15 &,Q15))is used as temporarily variables
[C++ error] Unit1.cpp(104): E2064 'Q15 &' is 'Q15 *' cant be initialized
[C++ error] Unit1.cpp(104): E2342 'x0' parameters is Q15 & type and cant be change to Q15 *
[C++ warning] Unit1.cpp(104): W8030 'y0' parameters(cord ic(Q15 &,Q15 &,Q15 &,Q15))is used as temporarily variables
[C++ error] Unit1.cpp(104): E2064 'Q15 &' is 'Q15 *' cant be initialized
[C++ error] Unit1.cpp(104): E2342 'y0' parameters is Q15 & type and cant be change to Q15 *
[C++ warning] Unit1.cpp(104): W8030 'z0' パラメータ(cordic(Q1 5 &,Q15 &,Q15 &,Q15))is used as temporarily variables
[C++ warning] Unit1.cpp(104): E2064 'Q15 &' は 'Q15 *' cant be initialized
[C++ error] Unit1.cpp(104): E2342 'z0' parameters is Q15 & type and cant be change to Q15 *
i`m thinking the main problem in my program is because i`m trying to use fixed point (its about the data type). I dont really understand that much.
*ps i`m using c++builder for japanese so the error message output is in japanese , i`m translating by myself to english (so if it`s sound weird please correct it or ask me)
horace1 1,510
Recognized Expert Top Contributor
the function definition of cordic() sepcifies that the formal parameters and return type are Q15 -
Q15 cordic(Q15 &x0, Q15 &y0, Q15 &z0, Q15 vecmode)
-
you call it with actual parameters of type float -
static float gain1Cordic()
-
{
-
/* compute gain by evaluating cos(0) without inv gain */
-
float x, y, z;
-
-
x = 1;
-
y = 0;
-
z = 0;
-
cordic(&x, &y, &z, 1.0f);-->THE MAIN PROBLEM,
-
return x; the error messages repeated so much in this line
-
}
-
where did you get the Q15 fixed point library from?
the function definition of cordic() sepcifies that the formal parameters and return type are Q15 -
Q15 cordic(Q15 &x0, Q15 &y0, Q15 &z0, Q15 vecmode)
-
you call it with actual parameters of type float -
static float gain1Cordic()
-
{
-
/* compute gain by evaluating cos(0) without inv gain */
-
float x, y, z;
-
-
x = 1;
-
y = 0;
-
z = 0;
-
cordic(&x, &y, &z, 1.0f);-->THE MAIN PROBLEM,
-
return x; the error messages repeated so much in this line
-
}
-
where did you get the Q15 fixed point library from?
my professor make it by himself the library, i included it in fixed_math.hpp
RedSon 5,000
Recognized Expert Expert
Can you post the .hpp file with your professors code in it? If it is large consider using the website mediafire.com to upload it and give us the link to it. Also, when you translate error messages, please try to translate word for word. What I mean is don't worry about changing the grammar and word order so that is makes sense in english. We can still get the meaning even though it may not look right, and that way there is less loss in translation. :)
MMcCarthy 14,534
Recognized Expert Moderator MVP
Can you post the .hpp file with your professors code in it? If it is large consider using the website mediafire.com to upload it and give us the link to it. Also, when you translate error messages, please try to translate word for word. What I mean is don't worry about changing the grammar and word order so that is makes sense in english. We can still get the meaning even though it may not look right, and that way there is less loss in translation. :)
Sorry guys I've already told the OP they can't as the code is copyrighted.
RedSon 5,000
Recognized Expert Expert
The professor's code is copyrighted? Or the CORDIC code?
The professor's code is copyrighted? Or the CORDIC code?
sorry my professor code is copyrighted
Sign in to post your reply or Sign up for a free account.
Similar topics |
by: Adrian Parker |
last post by:
I'm using the code below in my project. When I print all of these fixed
length string variables, one per line, they strings in questions do not
properly pad with 0s. strQuantity prints as " 4". Six spaces than the
value of intQuantity. This is correct.
But all the others end up being string objects of only 6 characters long
(with the exception of strTotal). The left most positions of the string
object are being padded with one...
|
by: Roger Leigh |
last post by:
Hello,
I'm writing a fixed-precision floating point class, based on the ideas
in the example fixed_pt class in the "Practical C++ Programming" book
by Steve Oualline (O' Reilly). This uses a long int to store the
value, and the precision (number of decimal points) is variable (it's
a templated class):
template <size_t _decimal_places = 4>
class FixedFloat {
|
by: Shi Mu |
last post by:
any python module to calculate sin, cos, arctan?
|
by: pout |
last post by:
What are the purposes of fixed-point? When should it be used?
I read:
#define Int2Fixed(x) (((long)(short)x) << 16)
and the fixed-point in 16.16 format. Does the 16 in the MACRO refer to
integer or decimal part? For example, if in 8.24, should the macro be:
#define Int2Fixed(x) (((long)(short)x) << 24)?
|
by: Davidlohr Bueso |
last post by:
Hi,
I've been trying to write some trigonometric functions of the math library
(sin, cos, tan, etc, etc.) in C and I've learned that instead of using
series, the best way to go is by using the CORDIC algorithm. The problem
is that I haven't been able to find the algorithm itself, just brief
descriptions of it. Does anyone have any experience with this? Thanks!
Dave
| |
by: Gerald |
last post by:
Recently, my program need to be run in embeded enviroment, and I cann't
use standard library. But I need to use arctan(x), so I implement it
like the following:
inline double pow(double x, size_t n) {
if (n == 0)
return 1;
else if (n % 2 == 0)
return pow(x * x, n >> 1);
else
|
by: riya1012 |
last post by:
hello guys,
I need some help from you. I am doing a DSP project and for that I need
to do some C coding for the conversion of sample data which is in
floating point representation to fixed point representation.
the sample data is in floating point like
0.224128
2.299965
0.448350
-1.779926
|
by: Dave Townsend |
last post by:
Hi,
I'm working on an implementation of the CORDIC
algorithm for a library of functions.
One of the benefits of the CORDIC algorithm is that you
can implement it with simple additions of numbers and
division of (doubles) by powers of two. However, there
doesn't seem to be any mechanism in the C/C++ language
to direct the compiler to do this efficient form of division,
|
by: astri |
last post by:
i`m doing thesis about comparing calculation of arctan by polynomial and CORDIC. I`ve read a lot of journal and books about CORDIC and this is what i understand.
1. make x and y
2. make iterations
3. process with cordic calculation
4. the arctan result is calculating by arctan(y/x).
what i m confused that in cordic equation there`s
Zi=Zo-arctan(2^-i)
|
by: Oralloy |
last post by:
Hello folks,
I am unable to find appropriate documentation on the type promotion of bit-fields when using the generalised comparison operator "<=>".
The problem is that using the GNU compilers, it seems that the internal comparison operator "<=>" tries to promote arguments from unsigned to signed.
This is as boiled down as I can make it.
Here is my compilation command:
g++-12 -std=c++20 -Wnarrowing bit_field.cpp
Here is the code in...
|
by: Hystou |
last post by:
Overview:
Windows 11 and 10 have less user interface control over operating system update behaviour than previous versions of Windows. In Windows 11 and 10, there is no way to turn off the Windows Update option using the Control Panel or Settings app; it automatically checks for updates and installs any it finds, whether you like it or not. For most users, this new feature is actually very convenient. If you want to control the update process,...
| |
by: agi2029 |
last post by:
Let's talk about the concept of autonomous AI software engineers and no-code agents. These AIs are designed to manage the entire lifecycle of a software development project—planning, coding, testing, and deployment—without human intervention. Imagine an AI that can take a project description, break it down, write the code, debug it, and then launch it, all on its own....
Now, this would greatly impact the work of software developers. The idea...
|
by: isladogs |
last post by:
The next Access Europe User Group meeting will be on Wednesday 1 May 2024 starting at 18:00 UK time (6PM UTC+1) and finishing by 19:30 (7.30PM).
In this session, we are pleased to welcome a new presenter, Adolph Dupré who will be discussing some powerful techniques for using class modules.
He will explain when you may want to use classes instead of User Defined Types (UDT). For example, to manage the data in unbound forms.
Adolph will...
|
by: conductexam |
last post by:
I have .net C# application in which I am extracting data from word file and save it in database particularly. To store word all data as it is I am converting the whole word file firstly in HTML and then checking html paragraph one by one.
At the time of converting from word file to html my equations which are in the word document file was convert into image.
Globals.ThisAddIn.Application.ActiveDocument.Select();...
|
by: TSSRALBI |
last post by:
Hello
I'm a network technician in training and I need your help.
I am currently learning how to create and manage the different types of VPNs and I have a question about LAN-to-LAN VPNs.
The last exercise I practiced was to create a LAN-to-LAN VPN between two Pfsense firewalls, by using IPSEC protocols.
I succeeded, with both firewalls in the same network. But I'm wondering if it's possible to do the same thing, with 2 Pfsense firewalls...
|
by: 6302768590 |
last post by:
Hai team
i want code for transfer the data from one system to another through IP address by using C# our system has to for every 5mins then we have to update the data what the data is updated we have to send another system
|
by: muto222 |
last post by:
How can i add a mobile payment intergratation into php mysql website.
| |
by: bsmnconsultancy |
last post by:
In today's digital era, a well-designed website is crucial for businesses looking to succeed. Whether you're a small business owner or a large corporation in Toronto, having a strong online presence can significantly impact your brand's success. BSMN Consultancy, a leader in Website Development in Toronto offers valuable insights into creating effective websites that not only look great but also perform exceptionally well. In this comprehensive...
| |