On 6 July 2017 at 17:10, Alejandro D Otero <[email protected]> wrote: > > # Set rhs function > ts.setRHSFunction(self.evalRHSFunction, self.g) > (where evalRHSFunction has the form f(self, ts, t, Vort, f)) >
Please double check that the RHS function "outputs" the result in the passed-in vector "f", my guess is that you are just returning a new vector that is eventually ignored. -- Lisandro Dalcin ============ Research Scientist Computer, Electrical and Mathematical Sciences & Engineering (CEMSE) Extreme Computing Research Center (ECRC) King Abdullah University of Science and Technology (KAUST) http://ecrc.kaust.edu.sa/ 4700 King Abdullah University of Science and Technology al-Khawarizmi Bldg (Bldg 1), Office # 0109 Thuwal 23955-6900, Kingdom of Saudi Arabia http://www.kaust.edu.sa Office Phone: +966 12 808-0459
