Commit 270bba63 authored by Emmanuel Branlard's avatar Emmanuel Branlard

Using scipy for speedup of cholesky with float32 precision

parent 8093b151
Pipeline #13619 passed with stage
in 1 minute and 13 seconds
......@@ -9,6 +9,7 @@ import warnings
import numpy as np
import pandas as pd
import scipy
from pyconturb.coherence import get_coh_mat
from pyconturb.core import TimeConstraint
......@@ -164,7 +165,7 @@ def gen_turb(spat_df, T=600, dt=1, con_tc=None, coh_model='iec',
all_mags[i_f, :]) * coh_mat
# get cholesky decomposition of sigma matrix
cor_mat = np.linalg.cholesky(sigma)
cor_mat = scipy.linalg.cholesky(sigma,overwrite_a=True, check_finite=False, lower=True)
# if constraints, assign data unc_pha
if constrained:
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment