Может ли значение yerr [] членская переменная класса эволюции в GSL_ODEIV2 (часть библиотеки GSL) быть отрицательным? ⇐ C++

[Anonymous]

Я учусь решать дифференциальные уравнения с помощью библиотеки GSL (gsl_odeiv2) и хотел отслеживать ошибки, допущенные во время вычислений, но, думаю, наткнулся на что-то странное.
Я создал объект эволюции gsl_odeiv2_evolve *evolve = gsl_odeiv2_evolve_alloc(dim);, а затем, внутри цикла while, я начал предпринимать шаги, используя status = gsl_odeiv2_evolve_apply(evolve, control, Step, &sys, &t, t1, &init_time_step, y); метод.
Согласно документации GSL этого конкретного метода yerr[] содержит оценку локальной ошибки, сделанной на одном конкретном шаге. Но в контексте пошагового метода говорится, что yerr[] хранит АБСОЛЮТНУЮ ошибку, что меня немного сбивает с толку. Все ошибки, которые я получаю при использовании gsl_odeiv2_evolve_apply, являются отрицательными, что кажется немного странным.
Кроме того, я не являюсь экспертом в C++ или библиотеке GSL, но я обнаружил, что yerr на самом деле является указателем на двойное значение, но он ведет себя как массив (я знаю, что в C он должен указывать на начало массива), но что меня беспокоит еще больше, так это то, что я могу безостановочно перебирать этот массив, т.е. компилятор не злится на yerr[20], даже несмотря на то, что этот массив должен содержать только одно значение (я проверил его размер с помощью sizeof(evolve->yerr)/sizeof(*(evolve->yerr)) и результат, который я получаю, — ОДИН элемент).
Кроме того, я не эксперт в C++ или библиотеке GSL, но я обнаружил, что yerr на самом деле является указателем для удвоения, но он ведет себя как а массив (я знаю, что в C он должен указывать на начало массива), но меня еще больше беспокоит то, что я могу перебирать этот массив без остановки, т. е. компилятор не злится на yerr[20], даже несмотря на то, что этот массив должен содержать только одно значение (я проверил его размер с помощью sizeof(evolve->yerr)/sizeof(*(evolve->yerr)) и результат, который я получаю, - ОДИН элемент).

Подробнее здесь: https://stackoverflow.com/questions/797 ... iv2-part-o