Eigenvector of a 4x4 matrix

Hello, i’m trying to find a way to compute the eigenvector corresponding to the highest eigenvalue of a 4x4 matrix. On host side i use a VNL function of the ITK toolkit, but an implementation in CUDA “could” be faster :).
since i only use a 4x4 matrix, i hope there is a relatively simple way to compute it. i only know an iterative way for a 3x3 matrix.

has anybody an idea how do do it?

best regards, tobias

Hi, the iterative (Power-iteration - search e.g. wiki) method works the same way for 4×4 than for 3×3.

It will calculate the biggest (in absolut value) eigenvalue and the corresponding eigenvector.

Laszlo Daroczy