Commit 02d58c08 authored by Ding, Lijie's avatar Ding, Lijie
Browse files

adjust param range based on prec results

parent c09617a8

analyze/ML_analyze.py

+2 −2
Original line number Diff line number Diff line
@@ -109,8 +109,8 @@ def calc_svd(folder, parameters): 


    #F = [np.concatenate((all_gq_r[i], all_gq_phi[i])) for i in range(len(all_gq_r))]

    F = all_gq_phi

    F = all_gq_r

    qphi = qr

    #F = all_gq_r

    #qphi = qr



    print("all_feature shape:", np.array(all_feature).shape)

    svd = np.linalg.svd(F)

analyze/main_ML_analyze.py

+2 −2
Original line number Diff line number Diff line
@@ -10,7 +10,7 @@ import time 
def main():



    print("analyzing data using ML model")

    folder = "../data/20241202_rand"

    folder = "../data/20241202_andes_rand"

    #folder = "../data/20241125_rand"

    rand_num = 5000

    rand_max = 2000
@@ -24,7 +24,7 @@ def main(): 
    print("parameters", parameters)

    print("total number of parameters", len(parameters))



    calc_svd(folder, parameters)

    #calc_svd(folder, parameters)



    #return 0

    random.shuffle(parameters)

code/suspension

(94.9 KiB)

File changed.

No diff preview for this file type.

View original file
View changed file

code/suspension.cpp

+2 −2
Original line number Diff line number Diff line
@@ -35,9 +35,9 @@ suspension::suspension(double R0_, double Rmu_, int n_, double sigma_, double sq 
        Rmu = 1.0 * R0;

        n = 100 + 100 * rand_uni(gen);

        // n = 100;

        sigma = 0.0 + 0.4 * rand_uni(gen);

        sigma = 0.0 + 0.3 * rand_uni(gen);

        sqrtD = (0 + 3 * rand_uni(gen)) * R0;

        gxy = (0 + 10 * rand_uni(gen)) * R0;

        gxy = (0 + 30 * rand_uni(gen)) * R0;

    }

    else

    {