473,395 Members | 1,694 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,395 software developers and data experts.

PDL::Slatec Problem of eigenvalue of Matrix

I have installed the PDL module and use it to calculate eigenvalue of
a
matrix. The PDL do gives me some eigenvalue and eigenvector. But the
problem is that these values are different from the eigenvalue and
eigenvector calculated in Matlab.

A simple example is:

use PDL::LiteF;
my $mat = pdl [ 3.4603 , 1.7529 , 1.6990 , 2.1181],
[ 2.5459 , 2.5615 , 2.2525 , 2.3606],
[ 2.3606 , 2.2515 , 2.5628 , 2.5459],
[ 2.1181 , 1.6982 , 1.7537 , 3.4603];
($eigvals,$eigvecs) = eigsys($mat);
print $eigvecs,$eigvals,"\n";
The results are:

For PDL module is:
[
[0.0541307 -0.554549 0.789897 -0.256142]
[0.0412208 0.635852 0.202431 -0.743649]
[ 0.87661 -0.252953 -0.320313 -0.254888]
[-0.476366 -0.473479 -0.482167 -0.562502]
] and eigenvalue
[0.272291 0.63145 1.7178 9.42336]
But for matlab, the results are

-0.4788 -0.6962 0.4294 -0.0371
-0.5203 -0.1238 -0.5616 0.7063
-0.5203 0.1239 -0.5619 -0.7059
-0.4788 0.6961 0.4295 0.0370
and the eigenvalue is
9.3294 0 0 0
0 1.3519 0 0
0 0 1.0631 0
0 0 0 0.3004
As you can see, the results from PDL is greatly different from Matlab.
So is my installation not correct? or something wrong with my perl
script?

Could somebody who has the PDL installed try this simple example, and
tell me
what is you eigenvalue and eigenvector?

Thanks
Jul 19 '05 #1
0 2291

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

Similar topics

2
by: B McInnes | last post by:
Hello, I am trying to store approximately 50- 100 million variable length strings into an array and of course this is not working due to memory constraints. I have converted the strings to unique...
0
by: ryjfgjl | last post by:
In our work, we often receive Excel tables with data in the same format. If we want to analyze these data, it can be difficult to analyze them because the data is spread across multiple Excel files...
0
by: emmanuelkatto | last post by:
Hi All, I am Emmanuel katto from Uganda. I want to ask what challenges you've faced while migrating a website to cloud. Please let me know. Thanks! Emmanuel
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...
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
marktang
by: marktang | last post by:
ONU (Optical Network Unit) is one of the key components for providing high-speed Internet services. Its primary function is to act as an endpoint device located at the user's premises. However,...
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
jinu1996
by: jinu1996 | last post by:
In today's digital age, having a compelling online presence is paramount for businesses aiming to thrive in a competitive landscape. At the heart of this digital strategy lies an intricately woven...
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
tracyyun
by: tracyyun | last post by:
Dear forum friends, With the development of smart home technology, a variety of wireless communication protocols have appeared on the market, such as Zigbee, Z-Wave, Wi-Fi, Bluetooth, etc. Each...

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.