I'm working on getting better code coverage on classes and have run into a bit of a snag on the DeepSDP4 class. There are two for loops of similiar structure, on on line 636 and one on line 688. Both of them essentially do:
for (int j = 0; j < SECULAR_INTEGRATION_ORDER; j += 2) {
...
if ((j + 2) < SECULAR_INTEGRATION_ORDER) {
...
}
}
With the current constant SECULAR_INTEGRATION_ORDER set to 2 we never get within the if loop. Is this there for the eventuality that the secular integration order is increased or of variable size or is this exposing an implementation problem?