POV-Ray : Newsgroups : povray.text.scene-files : [Minimum Volume] Bounding Ellipsoid via SVD : Re: [Minimum Volume] Bounding Ellipsoid via SVD Server Time
20 Sep 2021 22:51:21 EDT (-0400)
  Re: [Minimum Volume] Bounding Ellipsoid via SVD  
From: Bald Eagle
Date: 5 Nov 2019 20:10:00
Message: <web.5dc21d59f7b9a3af4eec112d0@news.povray.org>
I copied an pasted / modified the print_matrix () function into all of the other
functions, and puzzled out how to print out a vector ("list").

Got that to compile, and now we can step-by-step what happens, and where it goes
off the rails.   Looks like the trouble starts when trying to create the
Hermitian matrix.

In any event, this is way farther along than I'd thought I'd be by now...   :)

I suppose I need to dump everything to a text file, rather than to std::cout <<

Singular Value Decomposition (SVD):

Enter size of matrix N = (50x50 max): 10

A =
9.0 7.0 3.0 3.0 6.0 5.0 0.0 9.0 3.0 4.0
8.0 4.0 0.0 4.0 7.0 3.0 2.0 6.0 6.0 6.0
1.0 9.0 5.0 3.0 6.0 7.0 5.0 1.0 2.0 8.0
5.0 3.0 7.0 8.0 6.0 6.0 6.0 7.0 5.0 1.0
3.0 5.0 8.0 3.0 0.0 5.0 8.0 2.0 3.0 5.0
0.0 6.0 6.0 6.0 9.0 2.0 5.0 4.0 5.0 9.0
3.0 2.0 3.0 0.0 3.0 1.0 6.0 9.0 8.0 4.0
0.0 2.0 9.0 0.0 7.0 1.0 5.0 8.0 3.0 1.0
3.0 6.0 7.0 1.0 2.0 7.0 5.0 9.0 3.0 0.0
0.0 8.0 2.0 3.0 1.0 7.0 7.0 9.0 6.0 5.0


[matrix]transposed =
9.0 8.0 1.0 5.0 3.0 0.0 3.0 0.0 3.0 0.0
7.0 4.0 9.0 3.0 5.0 6.0 2.0 2.0 6.0 8.0
3.0 0.0 5.0 7.0 8.0 6.0 3.0 9.0 7.0 2.0
3.0 4.0 3.0 8.0 3.0 6.0 0.0 0.0 1.0 3.0
6.0 7.0 6.0 6.0 0.0 9.0 3.0 7.0 2.0 1.0
5.0 3.0 7.0 6.0 5.0 2.0 1.0 1.0 7.0 7.0
0.0 2.0 5.0 6.0 8.0 5.0 6.0 5.0 5.0 7.0
9.0 6.0 1.0 7.0 2.0 4.0 9.0 8.0 9.0 9.0
3.0 6.0 2.0 5.0 3.0 5.0 8.0 3.0 3.0 6.0
4.0 6.0 8.0 1.0 5.0 9.0 4.0 1.0 0.0 5.0


[matrices]multiplied =
315.0 265.0 214.0 259.0 167.0 229.0 194.0 173.0 230.0 231.0
265.0 266.0 195.0 234.0 147.0 235.0 194.0 142.0 169.0 206.0
214.0 195.0 295.0 224.0 220.0 281.0 148.0 159.0 196.0 242.0
259.0 234.0 224.0 330.0 222.0 266.0 209.0 219.0 252.0 250.0
167.0 147.0 220.0 222.0 234.0 214.0 158.0 157.0 200.0 217.0
229.0 235.0 281.0 266.0 214.0 340.0 201.0 212.0 192.0 247.0
194.0 194.0 148.0 209.0 158.0 201.0 229.0 183.0 190.0 223.0
173.0 142.0 159.0 219.0 157.0 212.0 183.0 234.0 202.0 178.0
230.0 169.0 196.0 252.0 200.0 192.0 190.0 202.0 263.0 250.0
231.0 206.0 242.0 250.0 217.0 247.0 223.0 178.0 250.0 318.0


[matrices]multiplied =
198.0 158.0 121.0 114.0 161.0 145.0 108.0 225.0 144.0 124.0
158.0 324.0 245.0 169.0 236.0 267.0 249.0 311.0 211.0 256.0
121.0 245.0 326.0 153.0 232.0 219.0 273.0 301.0 192.0 184.0
114.0 169.0 153.0 153.0 171.0 151.0 151.0 176.0 139.0 152.0
161.0 236.0 232.0 171.0 301.0 178.0 195.0 290.0 204.0 225.0
145.0 267.0 219.0 151.0 178.0 248.0 222.0 273.0 176.0 183.0
108.0 249.0 273.0 151.0 195.0 222.0 289.0 297.0 221.0 207.0
225.0 311.0 301.0 176.0 290.0 273.0 297.0 494.0 303.0 222.0
144.0 211.0 192.0 139.0 204.0 176.0 221.0 303.0 226.0 194.0
124.0 256.0 184.0 152.0 225.0 183.0 207.0 222.0 194.0 265.0


Eigenvectors =
0.7 1.2 1.1 0.8 1.1 1.0 1.1 1.5 1.0 1.0

Hermitian matrix =
1.4 nan nan nan nan nan nan nan nan nan
-1.7 1.0 nan nan nan nan nan nan nan nan
-1.5 nan 1.0 nan nan nan nan nan nan nan
-1.0 nan nan 1.0 nan nan nan nan nan nan
-1.5 nan nan nan 1.0 nan nan nan nan nan
-1.4 nan nan nan nan 1.0 nan nan nan nan
-1.5 nan nan nan nan nan 1.0 nan nan nan
-2.0 nan nan nan nan nan nan 1.0 nan nan
-1.4 nan nan nan nan nan nan nan 1.0 nan
-1.4 nan nan nan nan nan nan nan nan 1.0


[matrices]multiplied =
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan


Hermitian matrix inverse =
0.7 nan nan nan nan nan nan nan nan nan
-1.2 1.0 nan nan nan nan nan nan nan nan
-1.1 nan 1.0 nan nan nan nan nan nan nan
-0.8 nan nan 1.0 nan nan nan nan nan nan
-1.1 nan nan nan 1.0 nan nan nan nan nan
-1.0 nan nan nan nan 1.0 nan nan nan nan
-1.1 nan nan nan nan nan 1.0 nan nan nan
-1.5 nan nan nan nan nan nan 1.0 nan nan
-1.0 nan nan nan nan nan nan nan 1.0 nan
-1.0 nan nan nan nan nan nan nan nan 1.0


[matrices]multiplied =
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan


Reduced matrix =
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan


Eigenvectors =
nan nan nan nan nan nan nan nan 1.0

Hermitian matrix =
nan nan nan nan nan nan nan nan nan
nan 1.0 nan nan nan nan nan nan nan
nan nan 1.0 nan nan nan nan nan nan
nan nan nan 1.0 nan nan nan nan nan
nan nan nan nan 1.0 nan nan nan nan
nan nan nan nan nan 1.0 nan nan nan
nan nan nan nan nan nan 1.0 nan nan
nan nan nan nan nan nan nan 1.0 nan
nan nan nan nan nan nan nan nan 1.0


[matrices]multiplied =
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan


Hermitian matrix inverse =
nan nan nan nan nan nan nan nan nan
nan 1.0 nan nan nan nan nan nan nan
nan nan 1.0 nan nan nan nan nan nan
nan nan nan 1.0 nan nan nan nan nan
nan nan nan nan 1.0 nan nan nan nan
nan nan nan nan nan 1.0 nan nan nan
nan nan nan nan nan nan 1.0 nan nan
nan nan nan nan nan nan nan 1.0 nan
-1.0 nan nan nan nan nan nan nan 1.0


[matrices]multiplied =
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan


Reduced matrix =
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan


Eigenvectors =
nan nan nan nan nan nan nan 1.0

Hermitian matrix =
nan nan nan nan nan nan nan nan
nan 1.0 nan nan nan nan nan nan
nan nan 1.0 nan nan nan nan nan
nan nan nan 1.0 nan nan nan nan
nan nan nan nan 1.0 nan nan nan
nan nan nan nan nan 1.0 nan nan
nan nan nan nan nan nan 1.0 nan
nan nan nan nan nan nan nan 1.0


[matrices]multiplied =
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan


Hermitian matrix inverse =
nan nan nan nan nan nan nan nan
nan 1.0 nan nan nan nan nan nan
nan nan 1.0 nan nan nan nan nan
nan nan nan 1.0 nan nan nan nan
nan nan nan nan 1.0 nan nan nan
nan nan nan nan nan 1.0 nan nan
nan nan nan nan nan nan 1.0 nan
-1.0 nan nan nan nan nan nan 1.0


[matrices]multiplied =
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan


Reduced matrix =
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan


Eigenvectors =
nan nan nan nan nan nan 1.0

Hermitian matrix =
nan nan nan nan nan nan nan
nan 1.0 nan nan nan nan nan
nan nan 1.0 nan nan nan nan
nan nan nan 1.0 nan nan nan
nan nan nan nan 1.0 nan nan
nan nan nan nan nan 1.0 nan
nan nan nan nan nan nan 1.0


[matrices]multiplied =
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan


Hermitian matrix inverse =
nan nan nan nan nan nan nan
nan 1.0 nan nan nan nan nan
nan nan 1.0 nan nan nan nan
nan nan nan 1.0 nan nan nan
nan nan nan nan 1.0 nan nan
nan nan nan nan nan 1.0 nan
-1.0 nan nan nan nan nan 1.0


[matrices]multiplied =
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan
nan nan nan nan nan nan nan


Reduced matrix =
nan nan nan nan nan nan
nan nan nan nan nan nan
nan nan nan nan nan nan
nan nan nan nan nan nan
nan nan nan nan nan nan
nan nan nan nan nan nan


Eigenvectors =
nan nan nan nan nan 1.0

Hermitian matrix =
nan nan nan nan nan nan
nan 1.0 nan nan nan nan
nan nan 1.0 nan nan nan
nan nan nan 1.0 nan nan
nan nan nan nan 1.0 nan
nan nan nan nan nan 1.0


[matrices]multiplied =
nan nan nan nan nan nan
nan nan nan nan nan nan
nan nan nan nan nan nan
nan nan nan nan nan nan
nan nan nan nan nan nan
nan nan nan nan nan nan


Hermitian matrix inverse =
nan nan nan nan nan nan
nan 1.0 nan nan nan nan
nan nan 1.0 nan nan nan
nan nan nan 1.0 nan nan
nan nan nan nan 1.0 nan
-1.0 nan nan nan nan 1.0


[matrices]multiplied =
nan nan nan nan nan nan
nan nan nan nan nan nan
nan nan nan nan nan nan
nan nan nan nan nan nan
nan nan nan nan nan nan
nan nan nan nan nan nan


Reduced matrix =
nan nan nan nan nan
nan nan nan nan nan
nan nan nan nan nan
nan nan nan nan nan
nan nan nan nan nan


Eigenvectors =
nan nan nan nan 1.0

Hermitian matrix =
nan nan nan nan nan
nan 1.0 nan nan nan
nan nan 1.0 nan nan
nan nan nan 1.0 nan
nan nan nan nan 1.0


[matrices]multiplied =
nan nan nan nan nan
nan nan nan nan nan
nan nan nan nan nan
nan nan nan nan nan
nan nan nan nan nan


Hermitian matrix inverse =
nan nan nan nan nan
nan 1.0 nan nan nan
nan nan 1.0 nan nan
nan nan nan 1.0 nan
-1.0 nan nan nan 1.0


[matrices]multiplied =
nan nan nan nan nan
nan nan nan nan nan
nan nan nan nan nan
nan nan nan nan nan
nan nan nan nan nan


Reduced matrix =
nan nan nan nan
nan nan nan nan
nan nan nan nan
nan nan nan nan


Eigenvectors =
nan nan nan 1.0

Hermitian matrix =
nan nan nan nan
nan 1.0 nan nan
nan nan 1.0 nan
nan nan nan 1.0


[matrices]multiplied =
nan nan nan nan
nan nan nan nan
nan nan nan nan
nan nan nan nan


Hermitian matrix inverse =
nan nan nan nan
nan 1.0 nan nan
nan nan 1.0 nan
-1.0 nan nan 1.0


[matrices]multiplied =
nan nan nan nan
nan nan nan nan
nan nan nan nan
nan nan nan nan


Reduced matrix =
nan nan nan
nan nan nan
nan nan nan


Eigenvectors =
nan nan 1.0

Hermitian matrix =
nan nan nan
nan 1.0 nan
nan nan 1.0


[matrices]multiplied =
nan nan nan
nan nan nan
nan nan nan


Hermitian matrix inverse =
nan nan nan
nan 1.0 nan
-1.0 nan 1.0


[matrices]multiplied =
nan nan nan
nan nan nan
nan nan nan


Reduced matrix =
nan nan
nan nan


Eigenvectors =
nan 1.0

Hermitian matrix =
nan nan
nan 1.0


[matrices]multiplied =
nan nan
nan nan


Hermitian matrix inverse =
nan nan
-1.0 1.0


[matrices]multiplied =
nan nan
nan nan


Reduced matrix =
nan


Eigenvectors =
0.7 1.2 1.1 0.8 1.1 1.0 1.1 1.5 1.0 1.0

Eigenvectors =
nan nan nan nan nan nan nan nan nan 1.0

Eigenvectors =
nan nan nan nan nan nan nan nan nan 1.0

Eigenvectors =
nan nan nan nan nan nan nan nan nan 1.0

Eigenvectors =
nan nan nan nan nan nan nan nan nan 1.0

Eigenvectors =
nan nan nan nan nan nan nan nan nan 1.0

Eigenvectors =
nan nan nan nan nan nan nan nan nan 1.0

Eigenvectors =
nan nan nan nan nan nan nan nan nan 1.0

Eigenvectors =
nan nan nan nan nan nan nan nan nan 1.0

Eigenvectors =
nan nan nan nan nan nan nan nan nan 1.0

[matrix]transposed =
0.2 nan nan nan nan nan nan nan nan nan
0.4 nan nan nan nan nan nan nan nan nan
0.3 nan nan nan nan nan nan nan nan nan
0.2 nan nan nan nan nan nan nan nan nan
0.3 nan nan nan nan nan nan nan nan nan
0.3 nan nan nan nan nan nan nan nan nan
0.3 nan nan nan nan nan nan nan nan nan
0.4 nan nan nan nan nan nan nan nan nan
0.3 nan nan nan nan nan nan nan nan nan
0.3 nan nan nan nan nan nan nan nan nan


inverse diagonal matrix =
0.0 nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan


[matrices]multiplied =
15.5 nan nan nan nan nan nan nan nan nan
14.0 nan nan nan nan nan nan nan nan nan
14.8 nan nan nan nan nan nan nan nan nan
16.8 nan nan nan nan nan nan nan nan nan
13.1 nan nan nan nan nan nan nan nan nan
16.5 nan nan nan nan nan nan nan nan nan
13.0 nan nan nan nan nan nan nan nan nan
12.6 nan nan nan nan nan nan nan nan nan
14.5 nan nan nan nan nan nan nan nan nan
16.1 nan nan nan nan nan nan nan nan nan


[matrices]multiplied =
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan


S =
46.7 nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan


U =
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan
nan nan nan nan nan nan nan nan nan nan


V =
0.2 nan nan nan nan nan nan nan nan nan
0.4 nan nan nan nan nan nan nan nan nan
0.3 nan nan nan nan nan nan nan nan nan
0.2 nan nan nan nan nan nan nan nan nan
0.3 nan nan nan nan nan nan nan nan nan
0.3 nan nan nan nan nan nan nan nan nan
0.3 nan nan nan nan nan nan nan nan nan
0.4 nan nan nan nan nan nan nan nan nan
0.3 nan nan nan nan nan nan nan nan nan
0.3 nan nan nan nan nan nan nan nan nan


Post a reply to this message

Copyright 2003-2021 Persistence of Vision Raytracer Pty. Ltd.