std::vector<double>::const_iterator maxIter = std::max_element( array.begin(), array.end() );

        std::vector<double>::const_iterator it;

        // find first element on left of maximum lower than threshold
        it = std::lower_bound (maxIter, array.begin(), threshold);
        size_t leftindex = std::distance(array.begin(), it);
        double leftValue = axis[leftindex];

leider liefert mir das genau it = maxIter;

Kann es sein, das nicht von maxIter nach array.begin() gesucht wird weil begin kleiner ist als maxIter?

std::vector<double>::const_reverse_iterator rev_maxIter = maxIter;
it = std::lower_bound (rev_maxIter, array.rend(), threshold);

Ich hoffe, lower_bound kommt damit klar...

Alternativ halt die Schleife gerade selber schreiben.