static void test_inc_x(size_t i, const double *dx, const size_t *jdx, size_t ndx)
{
size_t j, k;
if (jdx) {
for (k = 0; k < ndx; k++) {
X[i * P + jdx[k]] += dx[k];
}
} else if (ndx) {
assert(ndx == P);
for (j = 0; j < P; j++) {
X[i * P + j] += dx[j];
}
}
recompute();
mlogit_set_x(&MLOGIT, i, 0, P, X + i * P);
//print_error("\tx...");
test_x();
//print_error("ok\n\tmean...");
test_mean();
//print_error("ok\n\tcov...");
test_cov();
//print_error("ok\n");
}
int main(){
int ind=0;
switch(ind){
case 0:
test_accuracy();
break;
case 1:
test_cov();
break;
case 2:
test_corner();
break;
case 3:
test_part();
break;
case 4:
test_stfun();
break;
case 5:
test_psd();
break;
case 6:
test_cxx();
}
}
