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; }