SHRot::SHRot(int numBands) {
numBands_ = numBands;
numCoeffs_ = numBands_*numBands_;
facs_ = std::vector<float>(2*numBands_);
facs_[0] = 1.f;
facs_[1] = 1.f;
for(size_t i = 2; i < facs_.size(); i++)
facs_[i] = i*facs_[i-1];
dySubDiag_ = std::vector<float>();
dydyDiag_ = std::vector<float>();
initDiagonals(numBands_);
}
int main(int argc, char *argv[])
{
int i, n;
long long sum = 1;
n = 1001;
initDiagonals(n);
// method 2
for (i = 3; i <= 1001; i += 2) {
sum += 4 * pow(i, 2) - 6 * i + 6;
}
printf("sum = %lld\n", sum);
return 0;
}
