Take a look at the literature that does similar things for facial recognition -- search for the term "eigenface", for example.
The point to make in this context is that the information you are looking for does not actually require you to consider high-resolution images. You may have $10,000\times 10,000$ pixels, for which any non-trivial ...
The eigenvalues given correspond to the flux jacobians in the Euler equations. Physically, they represent the speed at which waves of information can travel in a space-time domain.
The goal is to compute the eigenvalues at each face based on the values in each neighboring cell so there is no 'left' or 'right' face in this equation, but instead a left and ...
Let us proceed systematically:
numerical precision of data (you said from medical imaging)
number of operations required for standard methods (as from libraries)
possible out-of-core computation (i.e. not the whole matrix at all times in memory).
In all cases, I am afraid, you would have to be prepared to suffer. Incidentally, out-of-core methods are very ...
As noted above by Thijs Steel a randomized svd is a solution but the number 78800000 is out of our computers computation ability.So you can proceed to the rsvd algorithm by :
import numpy as np
n = 788
mu = 0
sigma = 1
A = np.random.normal(mu, sigma, (n,n))
Omega = np.random.normal(mu, sigma, (n,n))
def rsvd(A, Omega):
Y = A @ Omega
Q, _ = np....