コード例 #1
0
ファイル: rossler.cpp プロジェクト: vspetrov/nldyn
void Rossler::jac(const_it_t &state,
                  state_t & out, double time) {
    assert(out.size() == dim*dim);
    auto it = out.begin();
    it[0] = 0;        it[1] = -1;  it[2] = -1; it += dim;
    it[0] = 1;        it[1] = a;   it[2] = 0; it += dim;
    it[0] = state[2]; it[1] = 0;   it[2] = state[0]-c;
}