diff --git a/run_secorder.cpp b/run_secorder.cpp index 9da3dd3bfde2a646e898191cfcd73dd15a26ee26..02bde284732fcc78d721a610d8e34e2e868e01c3 100644 --- a/run_secorder.cpp +++ b/run_secorder.cpp @@ -86,7 +86,7 @@ int main(int argc, char *argv[]) { for(int i=0; i1.0); } fprintf(fhnd,"\n"); } diff --git a/secorder_rec_1p.cpp b/secorder_rec_1p.cpp index f18aaa34923f878465b5c48906081334e0cee87d..0a921ce41804bd7163f327998633b4de97f313a4 100644 --- a/secorder_rec_1p.cpp +++ b/secorder_rec_1p.cpp @@ -207,30 +207,7 @@ gsl_matrix* secorder_rec_1p(int N_nodes, double p, } - //////////////////////////////////////////////////////////// - // Step 5: Calculate Bernoulli matrix - // Simply make the Bernoulli variable be 1 - // if the Gaussian variable is greater than 1 - //////////////////////////////////////////////////////////// - - cout << "Generating Bernoulli matrix..."; - cout.flush(); - // calculate bernoulli matrix - gsl_matrix *W_ber = gsl_matrix_alloc(N_nodes, N_nodes); - for(int i=0; i1.0); - } - } - - // free Gaussian matrix - gsl_matrix_free(W_gaus); - - cout << "done\n"; - cout.flush(); - - // return Bernoulli matrix - return W_ber; + return W_gaus; }