spectrum_calc.cpp

// Hami->FreqsFiltered  is a *double
// Hami->IntensFiltered  is a *double
// TheoreticalSpec.Points is a *double
double LW = 2, Offset = 2000, FreqStep = 1 , tmp = 0 , tmp1 = 0, sqLW = LW * LW, CurrFreq=0;

for (int i = 1; i <= 131072; i++) {
        CurrFreq = Offset - FreqStep * (i - 1);
        tmp1 = 0;
        for (int j = 1; j <= Hami->nFreqsFiltered; j++) {
            tmp = CurrFreq - Hami->FreqsFiltered[j];
            tmp *= tmp;
            tmp1 += Hami->IntensFiltered[j] * LW / (tmp + sqLW);
        }
        TheoreticalSpec.Points[i] = tmp1;
    }