473,405 Members | 2,210 Online
Bytes | Software Development & Data Engineering Community
Post Job

Home Posts Topics Members FAQ

Join Bytes to post your question to a community of 473,405 software developers and data experts.

Complex number prints wierd

I'm trying to print a complex number and i get this as output: Phi:
(1.618034,0.000000).
What does it mean? How do I get it to print normally? Here's the code:

----------fib.hpp-------------
#ifndef FIB_HPP
#define FIB_HPP
#include <iostream>
#include <iomanip>
#include <cstdlib>
#include <complex>
using namespace std;

namespace
{
template <long N>
class Fib
{
public:
static const long double val=Fib<N-2>::val+Fib<N-1>::val;
};

template<>
class Fib<2>
{
public:
static const long double val=1;
};

template<>
class Fib<1>
{
public:
static const long double val=1;
};

const complex<long double> phi= Fib<32>::val/Fib<31>::val;
}
#endif
--------------------------

------Main.cpp------
#include "fib.hpp"

int main()
{
cout << "Phi: " << fixed << phi << endl;
system("PAUSE");
return 0;
}
------------------------

Can you help me?

Nov 8 '05 #1
9 1812
Protoman wrote:
I'm trying to print a complex number and i get this as output: Phi:
(1.618034,0.000000).
What does it mean?
It means that your complex number has a real part of 1.618034 and an
imaginary part of 0.0.

How do I get it to print normally?


This is the way complex numbers print *normally*. What did you expect?
[code snipped]
Best

Kai-Uwe Bux
Nov 8 '05 #2
OK. Is this a fast way to get fibonacci numbers?

Nov 8 '05 #3
You are calculating the fibonacci of 32 and 31 at compile time. So,
when the program runs, it doesn't have to calculate the sequence.
From the executable's standpoint: yes, this is fast.


Nov 8 '05 #4
Is there anyway I can make it faster? Do you know if there's a name for
what I'm doing for Fib?

Nov 8 '05 #5

Chris Goller wrote:
You are calculating the fibonacci of 32 and 31 at compile time. So,
when the program runs, it doesn't have to calculate the sequence.


There is no such thing as the Fibonacci of a number.

Granted, there is 32nd and 31st Fibonacci number, but this program is
certainly not calculating either of them. Nor is it successfully
calculating the golden mean.

This program is simply dividing the sum of all integers from 1 to 32 by
the sum of all integers from 1 to 31.

Greg

Nov 8 '05 #6
Greg wrote:

Chris Goller wrote:
You are calculating the fibonacci of 32 and 31 at compile time. So,
when the program runs, it doesn't have to calculate the sequence.
There is no such thing as the Fibonacci of a number.

Granted, there is 32nd and 31st Fibonacci number, but this program is
certainly not calculating either of them. Nor is it successfully
calculating the golden mean.


Your assesment is a little off, at least if you are referring to the program
in the original post: if the program compiled, it would compute Fibonacci
numbers and the golden mean; the problem with the code is that it uses
in-place initialization for non-integral types. Here is a version that
compiles. The algorithm is completely unchanged.

#include <iostream>
#include <iomanip>
#include <cstdlib>
#include <complex>
using namespace std;

template <long N>
class Fib
{
public:
static const unsigned long val=Fib<N-2>::val+Fib<N-1>::val;
};

template<>
class Fib<2>
{
public:
static const unsigned long val=1;
};

template<>
class Fib<1>
{
public:
static const unsigned long val=1;
};

const double phi= double(Fib<32>::val)/double(Fib<31>::val);

int main()
{
cout << "Phi: " << fixed << phi << endl;
return 0;
}

Guess what it prints:

Phi: 1.618034

That passes for an approximation of the golden mean. Moreover, you can print
the first few terms of the sequence:

#include <iostream>
#include <iomanip>
#include <cstdlib>
#include <complex>
using namespace std;

template <long N>
class Fib
{
public:
static const unsigned long val=Fib<N-2>::val+Fib<N-1>::val;
};

template<>
class Fib<2>
{
public:
static const unsigned long val=1;
};

template<>
class Fib<1>
{
public:
static const unsigned long val=1;
};

const double phi= double(Fib<32>::val)/double(Fib<31>::val);

int main()
{
cout << Fib<1>::val << " "
<< Fib<2>::val << " "
<< Fib<3>::val << " "
<< Fib<4>::val << " "
<< Fib<5>::val << " "
<< Fib<6>::val << " "
<< Fib<7>::val << "\n";
return 0;
}
And that prints:

1 1 2 3 5 8 13

Looks like Fibonacci numbers to me.
Also, it is actually clear from the code that the program computes the
Fibonacci sequence.

This program is simply dividing the sum of all integers from 1 to 32 by
the sum of all integers from 1 to 31.


No it does not. The sum of all integers from 1 to 32 is 16*31. The sum of
all integers from 1 to 31 is 15*31. The quotient is 16/15.

Best

Kai-Uwe Bux
Nov 8 '05 #7
So, how does the compiler execute the metaprogram?

Nov 8 '05 #8
Protoman wrote:

So, how does the compiler execute the metaprogram?


Well. If *you* were the compiler. What would *you* need to
do in order to create an executable program.
(Actually: A surprisingly large number of so called 'clever'
programming is just an adoption of what people do in real
life. So if you direct your thinking into 'What would I do
if all I have is paper and pencil?' is surprisingly often
an excellent tool in understanding what is going on.)

The compiler comes to compiling:

const complex<long double> phi= Fib<32>::val/Fib<31>::val;

This is a declaration of a complex<double> variable called phi.
The variable is initialized with the value of
Fib<32>::val/Fib<31>::val

For this the compiler needs to know
Fib<32>::val
and Fib<31>::val

Fib<32> is the instantiation of a template. Thus the compiler
looks up the Fib template and substitutes 32 for N

template <long N>
class Fib
{
public:
static const long double val=Fib<N-2>::val+Fib<N-1>::val;
};

becomes

class Fib
{
public:
static const long double val=Fib<30>::val+Fib<31>::val;
};

In order to make this compilable the compiler has to come up with
the initial value for val. For this it has to evaluate the initialization
part:

Fib<30>::val+Fib<31>::val

Again: The compiler is looking for an instantiation of Fib<30> and since
there is none it creates one, substituting 30 for N

class Fib
{
public:
static const long double val=Fib<28>::val+Fib<29>::val;
};

Part of compiling that one, makes the compiler look for an instatiation
of Fib<28>. Since there is none, the compiler creates one, using 28 for N

class Fib
{
public:
static const long double val=Fib<26>::val+Fib<27>::val;
};

And so on, and so on.
Finally the compiler will have the request to intialize one of the
generated classes with:

class Fib
{
public:
static const long double val=Fib<1>::val+Fib<2>::val;
};

Looking for Fib<1> the compiler figures out, that the programmer specified
that one:

template<>
class Fib<1>
{
public:
static const long double val=1;
};

so it doesn't need to create one on its own. Same for Fib<2>.

Since those 2 classes are 'complete' the compiler now can use
them to continue working on

class Fib
{
public:
static const long double val=Fib<1>::val+Fib<2>::val;
};

(which was the template instatiation for Fib<3>).
Fib<1>::val is known, it equals 1.
Fib<2>::val is known, it equals 1
Thus the initialization value for Fib<3>::val thus must be 2

Having this value, the compiler can continue on finishing
Fib<4> and so on, and so on, until finally Fib<32>::val can
be calculated by the compiler.
Now the whole story starts again for Fib<31>.

If both values Fib<32>::val and Fib<31>::val are known to the
copmiler, it can easily calculate Fib<32>::val / Fib<31>::val
and assign that value to phi.
However why one would want to do all of that with complex numbers
is bejond my imagination.

--
Karl Heinz Buchegger
kb******@gascad.at
Nov 8 '05 #9
So that's how metaprogramming works. Thanks!!!

Nov 9 '05 #10

This thread has been closed and replies have been disabled. Please start a new discussion.

Similar topics

13
by: Chris Cioffi | last post by:
Hello all! I'm trying to subclass int such that once it reaches a certain value, it flips back to the starting count. This is for use as an X12 control number. Subclassing int and getting the...
116
by: Mike MacSween | last post by:
S**t for brains strikes again! Why did I do that? When I met the clients and at some point they vaguely asked whether eventually would it be possible to have some people who could read the data...
6
by: Herrcho | last post by:
in K&R Chapter 6.3 it mentions two methods to calculate NKEYS. and points out the first one which is to terminate the list of initializers with a null pointer, then loop along keytab until the...
3
by: Russ | last post by:
I'd like to get output formatting for my own classes that mimics the built-in output formatting. For example, >>> x = 4.54 >>> print "%4.2f" % x 4.54 In other words, if I substitute a class...
12
by: vj | last post by:
Hi! I have a piece of code (shown below) involving complex numbers. The code is not running and giving error ("Invalid floating point operation" and "SQRT:Domain error"). I would be very...
2
by: wudoug119 | last post by:
This is my code and it will take any number that I input and say it is a prime number. Please help me... int Prime(int prime) //declares isPrime as a function using integers { ...
4
by: SweetLeftFoot | last post by:
Hello, i have designed some code that works out the first 250 prime numbers and prints them to the screen. However i need to implement 2 functions, one of which returns a 1 if the number is a prime...
7
by: Tina | last post by:
Dear all, I'm looking for a routine which calculates the Gamma Function for a complex valued variable. I'm using #include <complex> to work with complex numbers. I define a complex number...
25
by: jacob navia | last post by:
The C99 standard forgot to define the printf equivalent for complex numbers Since I am revising the lcc-win implementation of complex numbers I decided to fill this hole with "Z" for...
0
by: Charles Arthur | last post by:
How do i turn on java script on a villaon, callus and itel keypad mobile phone
0
BarryA
by: BarryA | last post by:
What are the essential steps and strategies outlined in the Data Structures and Algorithms (DSA) roadmap for aspiring data scientists? How can individuals effectively utilize this roadmap to progress...
1
by: Sonnysonu | last post by:
This is the data of csv file 1 2 3 1 2 3 1 2 3 1 2 3 2 3 2 3 3 the lengths should be different i have to store the data by column-wise with in the specific length. suppose the i have to...
0
by: Hystou | last post by:
There are some requirements for setting up RAID: 1. The motherboard and BIOS support RAID configuration. 2. The motherboard has 2 or more available SATA protocol SSD/HDD slots (including MSATA, M.2...
0
by: Hystou | last post by:
Most computers default to English, but sometimes we require a different language, especially when relocating. Forgot to request a specific language before your computer shipped? No problem! You can...
0
Oralloy
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,...
0
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...
0
agi2029
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,...
0
isladogs
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...

By using Bytes.com and it's services, you agree to our Privacy Policy and Terms of Use.

To disable or enable advertisements and analytics tracking please visit the manage ads & tracking page.