247 resolvent analysis - #271
msilvestrec03 wants to merge 26 commits into
Conversation
|
Examples are missing (to be done) |
Example of the jet uploaded |
|
Possible error related to the cholesky factorization. Do not merge |
ArnauMiro
left a comment
ArnauMiro
left a comment
There was a problem hiding this comment.
We should also think some testsuite runs to test the RES module.
| AC_MAT(B,m,ii,ii) = A[ii]; | ||
| } | ||
| else{ | ||
| AC_MAT(B,m,ii,jj) = 0; |
There was a problem hiding this comment.
Preallocate B to zeros with memset or similar and then just loop the diagonal.
| AC_MAT(B,m,ii,ii) = A[ii]; | ||
| } | ||
| else{ | ||
| AC_MAT(B,m,ii,jj) = 0; |
There was a problem hiding this comment.
Preallocate B to zeros with memset or similar and then just loop the diagonal.
| AC_MAT(B,m,ii,ii) = A[ii]; | ||
| } | ||
| else{ | ||
| AC_MAT(B,m,ii,jj) = 0; |
There was a problem hiding this comment.
Preallocate B to zeros with memset or similar and then just loop the diagonal.
| AC_MAT(B,m,ii,ii) = A[ii]; | ||
| } | ||
| else{ | ||
| AC_MAT(B,m,ii,jj) = 0; |
There was a problem hiding this comment.
Preallocate B to zeros with memset or similar and then just loop the diagonal.
|
|
||
| # Import Low Order Models | ||
| from . import POD, PCA, DMD, SPOD, MANIFOLD, GPOD, LAMINE | ||
| from . import POD, PCA, DMD, SPOD, MANIFOLD, GPOD, RES |
| from ..utils import gpu_to_cpu | ||
|
|
||
|
|
||
| def plotEnergy(S,fig=None,ax=None): |
There was a problem hiding this comment.
Isn't this the same as in POD?
There was a problem hiding this comment.
It is the same. I think that the function is also useful here. Should I change anything?
| ''' | ||
| p = cp if type(Phi) is cp.ndarray else np | ||
|
|
||
| # Normalization of modes (?) |
| # for ii in range(len(Matrix[0,:])): | ||
| # Matrix[:,ii] = Matrix[:,ii] / vector_norm(Phi[:,ii]) | ||
|
|
||
| Omega = delta + 1j * freq |
There was a problem hiding this comment.
done (enough?)
| U_res = matmul(Phi, matmul(Fhat_inv, U)) | ||
| V_res = matmul(Phi, matmul(Fhat_inv, V)) | ||
|
|
||
| # # U, S, VT = svd(diag(H)) |
| from libc.string cimport memcpy, memset | ||
| from libc.math cimport sqrt, log, atan2 | ||
| from ..vmmath.cfuncs cimport real, real_complex | ||
| # from ..vmmath.cfuncs cimport |
|
I don't know why is this failing to build intel |
|
@msilvestrec03 it's a timeout on the Intel APT repo. Possibly rerun the action after a while will fix it. |
Okay thanks. I will |
|
Don't merge yet. Wait a meeting with Rafel, as the algorithm could be improved. |
|
@msilvestrec03 @RafelPerello set the merge to draft if you're not ready please |
Updated platform, FFTW, and optimization settings.
…/pyLowOrder into 247-resolvent-analysis
This branch has not been deployed
2 participants
Resolvent analysis for pyLOM