Commits
1476 1476 | } |
1477 1477 | |
1478 1478 | void |
1479 1479 | print_max_gsl3(gsl_vector *v) |
1480 1480 | { |
1481 1481 | double phi_max = 0.0; |
1482 1482 | double del_max = 0.0; |
1483 1483 | double rat_max = 0.0; |
1484 1484 | |
1485 1485 | const size_t n = v->size; |
1486 - | for (int i=0; i!=n/3; i++) { |
1486 + | for (size_t i=0; i!=n/3; i++) { |
1487 1487 | if (fabs(gsl_vector_get(v, 3*i+0)) > fabs(phi_max)) phi_max = gsl_vector_get(v, 3*i+0); |
1488 1488 | if (fabs(gsl_vector_get(v, 3*i+1)) > fabs(del_max)) del_max = gsl_vector_get(v, 3*i+1); |
1489 1489 | if (fabs(gsl_vector_get(v, 3*i+2)) > fabs(rat_max)) rat_max = gsl_vector_get(v, 3*i+2); |
1490 1490 | } |
1491 1491 | cerr << "phi_max " << phi_max << " del_max " << del_max << " rat_max " << rat_max << endl; |
1492 1492 | } |
1493 1493 | |
1494 1494 | |
1495 1495 | |
1496 1496 | /* |