 // dimension of the characteristics, chi
 static const int p = 3;
 static __constant__ int c_p = 3;
+static const int dim_cholfactor = 6; // = p * p - ((p - 1) * p) / 2;
+static const int c_dim_cholfactor = 6;
 
 // dimension of the chi population parameter, theta
 static const int dim_theta = 6;

        CholUpdateR1(cholfactor, scaled_proposal, c_p, downdate);
        
        // copy cholesky factor for this data point back to global memory
        cholfact_index = 0
        for (int j=0; j<dim_cholfactor; j++) {
            prop_cholfact[n * j + i] = cholfactor[j];
        }

    // Construct initial cholesky factor for the covariance matrix of the theta proposal distribution
    int dim_cholfactor_theta = dim_theta * dim_theta - ((dim_theta - 1) * dim_theta) / 2;
    thrust::host_vector<double> h_theta_cholfact(dim_cholfactor_theta);
    int diag_index = 0;
    for (int j=0; j<dim_theta; j++) {
        h_theta_cholfact[diag_index] = 1e-2;
        diag_index += j + 2;
    }
    
