On 9 December 2014 at 18:24, Carles Bona <carlesb...@gmail.com> wrote: > TS type is alpha >
I'm working on a reimplementation of TSALPHA in PetIGA, https://bitbucket.org/dalcinl/petiga/src/default/src/tsalpha1.c . Eventually, I'll update the code in PETSc with this new implementation. This new implementation features some attempts to implement time-step adaptivity. Additionally, it implements a poor man's attempt to estimate a initial derivative without requiring users to setup and solve a problem involving the mass matrix. We can experiment in this code with some extra APIs to let users specify the initial derivative. What about TSAlphaSetSolution(ts,U,V) and TSAlphaGetSolution(ts,&U,&V), where U and V are the initial solution and derivative vectors? -- Lisandro Dalcin ============ Research Scientist Computer, Electrical and Mathematical Sciences & Engineering (CEMSE) Numerical Porous Media Center (NumPor) King Abdullah University of Science and Technology (KAUST) http://numpor.kaust.edu.sa/ 4700 King Abdullah University of Science and Technology al-Khawarizmi Bldg (Bldg 1), Office # 4332 Thuwal 23955-6900, Kingdom of Saudi Arabia http://www.kaust.edu.sa Office Phone: +966 12 808-0459