/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c

Bug Summary

File:	ts/impls/glee/glee.c
Warning:	line 850, column 10 Access to field 's' results in a dereference of a null pointer (loaded from variable 'tab')

Annotated Source Code

[?] Use j/k keys for keyboard navigation

Code for time stepping with the General Linear with Error Estimation method

Notes:

The general system is written as

Udot = F(t,U)

#include <petsc/private/tsimpl.h> /*I "petscts.h" I*/

#include <petscdm.h>

static PetscBool cited = PETSC_FALSE;

static const char citation[] =

"@ARTICLE{Constantinescu_TR2016b,\n"

" author = {Constantinescu, E.M.},\n"

" title = {Estimating Global Errors in Time Stepping},\n"

" journal = {ArXiv e-prints},\n"

" year = 2016,\n"

" adsurl = {http://adsabs.harvard.edu/abs/2015arXiv150305166C}\n}\n";

static TSGLEEType TSGLEEDefaultType = TSGLEE35"35";

static PetscBool TSGLEERegisterAllCalled;

static PetscBool TSGLEEPackageInitialized;

static PetscInt explicit_stage_time_id;

typedef struct _GLEETableau *GLEETableau;

struct _GLEETableau {

char *name;

PetscInt order; /* Classical approximation order of the method i*/

PetscInt s; /* Number of stages */

PetscInt r; /* Number of steps */

PetscReal gamma; /* LTE ratio */

PetscReal *A,*B,*U,*V,*S,*F,*c; /* Tableau */

PetscReal *Fembed; /* Embedded final method coefficients */

PetscReal *Ferror; /* Coefficients for computing error */

PetscReal *Serror; /* Coefficients for initializing the error */

PetscInt pinterp; /* Interpolation order */

PetscReal *binterp; /* Interpolation coefficients */

PetscReal ccfl; /* Placeholder for CFL coefficient relative to forward Euler */

};

typedef struct _GLEETableauLink *GLEETableauLink;

struct _GLEETableauLink {

struct _GLEETableau tab;

GLEETableauLink next;

};

static GLEETableauLink GLEETableauList;

typedef struct {

GLEETableau tableau;

Vec *Y; /* Solution vector (along with auxiliary solution y~ or eps) */

Vec *X; /* Temporary solution vector */

Vec *YStage; /* Stage values */

Vec *YdotStage; /* Stage right hand side */

Vec W; /* Right-hand-side for implicit stage solve */

Vec Ydot; /* Work vector holding Ydot during residual evaluation */

Vec yGErr; /* Vector holding the global error after a step is completed */

PetscScalar *swork; /* Scalar work (size of the number of stages)*/

PetscScalar *rwork; /* Scalar work (size of the number of steps)*/

PetscReal scoeff; /* shift = scoeff/dt */

PetscReal stage_time;

TSStepStatus status;

} TS_GLEE;

/*MC

TSGLEE23 - Second order three stage GLEE method

This method has three stages.

s = 3, r = 2

Level: advanced

.seealso: TSGLEE

M*/

/*MC

TSGLEE24 - Second order four stage GLEE method

This method has four stages.

s = 4, r = 2

Level: advanced

.seealso: TSGLEE

M*/

/*MC

TSGLEE25i - Second order five stage GLEE method

This method has five stages.

s = 5, r = 2

Level: advanced

.seealso: TSGLEE

M*/

/*MC

TSGLEE35 - Third order five stage GLEE method

This method has five stages.

100

s = 5, r = 2

101

102

Level: advanced

103

104

.seealso: TSGLEE

105

M*/

106

/*MC

107

TSGLEEEXRK2A - Second order six stage GLEE method

108

109

This method has six stages.

110

s = 6, r = 2

111

112

Level: advanced

113

114

.seealso: TSGLEE

115

M*/

116

/*MC

117

TSGLEERK32G1 - Third order eight stage GLEE method

118

119

This method has eight stages.

120

s = 8, r = 2

121

122

Level: advanced

123

124

.seealso: TSGLEE

125

M*/

126

/*MC

127

TSGLEERK285EX - Second order nine stage GLEE method

128

129

This method has nine stages.

130

s = 9, r = 2

131

132

Level: advanced

133

134

.seealso: TSGLEE

135

M*/

136

137

/*@C

138

TSGLEERegisterAll - Registers all of the General Linear with Error Estimation methods in TSGLEE

139

140

Not Collective, but should be called by all processes which will need the schemes to be registered

141

142

Level: advanced

143

144

.seealso: TSGLEERegisterDestroy()

145

@*/

146

PetscErrorCode TSGLEERegisterAll(void)

147

{

148

PetscErrorCode ierr;

149

150

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 150; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

151

if (TSGLEERegisterAllCalled) PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);

152

TSGLEERegisterAllCalled = PETSC_TRUE;

153

154

{

155

#define GAMMA0.25 0.5

156

/* y-eps form */

157

const PetscInt

158

p = 1,

159

s = 3,

160

r = 2;

161

const PetscReal

162

A[3][3] = {{1.0,0,0},{0,0.5,0},{0,0.5,0.5}},

163

B[2][3] = {{1.0,0,0},{-2.0,1.0,1.0}},

164

U[3][2] = {{1.0,0},{1.0,0.5},{1.0,0.5}},

165

V[2][2] = {{1,0},{0,1}},

166

S[2] = {1,0},

167

F[2] = {1,0},

168

Fembed[2] = {1,1-GAMMA0.25},

169

Ferror[2] = {0,1},

170

Serror[2] = {1,0};

171

ierr = TSGLEERegister(TSGLEEi1"BE1",p,s,r,GAMMA0.25,&A[0][0],&B[0][0],&U[0][0],&V[0][0],S,F,NULL((void*)0),Fembed,Ferror,Serror,0,NULL((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),171,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

172

}

173

{

174

#undef GAMMA0.25

175

#define GAMMA0.25 0.0

176

/* y-eps form */

177

const PetscInt

178

p = 2,

179

s = 3,

180

r = 2;

181

const PetscReal

182

A[3][3] = {{0,0,0},{1,0,0},{0.25,0.25,0}},

183

B[2][3] = {{1.0/12.0,1.0/12.0,5.0/6.0},{1.0/12.0,1.0/12.0,-1.0/6.0}},

184

U[3][2] = {{1,0},{1,10},{1,-1}},

185

V[2][2] = {{1,0},{0,1}},

186

S[2] = {1,0},

187

F[2] = {1,0},

188

Fembed[2] = {1,1-GAMMA0.25},

189

Ferror[2] = {0,1},

190

Serror[2] = {1,0};

191

ierr = TSGLEERegister(TSGLEE23"23",p,s,r,GAMMA0.25,&A[0][0],&B[0][0],&U[0][0],&V[0][0],S,F,NULL((void*)0),Fembed,Ferror,Serror,0,NULL((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),191,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

192

}

193

{

194

#undef GAMMA0.25

195

#define GAMMA0.25 0.0

196

/* y-y~ form */

197

const PetscInt

198

p = 2,

199

s = 4,

200

r = 2;

201

const PetscReal

202

A[4][4] = {{0,0,0,0},{0.75,0,0,0},{0.25,29.0/60.0,0,0},{-21.0/44.0,145.0/44.0,-20.0/11.0,0}},

203

B[2][4] = {{109.0/275.0,58.0/75.0,-37.0/110.0,1.0/6.0},{3.0/11.0,0,75.0/88.0,-1.0/8.0}},

204

U[4][2] = {{0,1},{75.0/58.0,-17.0/58.0},{0,1},{0,1}},

205

V[2][2] = {{1,0},{0,1}},

206

S[2] = {1,1},

207

F[2] = {1,0},

208

Fembed[2] = {0,1},

209

Ferror[2] = {-1.0/(1.0-GAMMA0.25),1.0/(1.0-GAMMA0.25)},

210

Serror[2] = {1.0-GAMMA0.25,1.0};

211

ierr = TSGLEERegister(TSGLEE24"24",p,s,r,GAMMA0.25,&A[0][0],&B[0][0],&U[0][0],&V[0][0],S,F,NULL((void*)0),Fembed,Ferror,Serror,0,NULL((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),211,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

212

}

213

{

214

#undef GAMMA0.25

215

#define GAMMA0.25 0.0

216

/* y-y~ form */

217

const PetscInt

218

p = 2,

219

s = 5,

220

r = 2;

221

const PetscReal

222

A[5][5] = {{0,0,0,0,0},

223

{-0.94079244066783383269,0,0,0,0},

224

{0.64228187778301907108,0.10915356933958500042,0,0,0},

225

{-0.51764297742287450812,0.74414270351096040738,-0.71404164927824538121,0,0},

226

{-0.44696561556825969206,-0.76768425657590196518,0.20111608138142987881,0.93828186737840469796,0}},

227

B[2][5] = {{-0.029309178948150356153,-0.49671981884013874923,0.34275801517650053274,0.32941112623949194988,0.85385985637229662276},

228

{0.78133219686062535272,0.074238691892675897635,0.57957363498384957966,-0.24638502829674959968,-0.18875949544040123033}},

229

U[5][2] = {{0.16911424754448327735,0.83088575245551672265},

230

{0.53638465733199574340,0.46361534266800425660},

231

{0.39901579167169582526,0.60098420832830417474},

232

{0.87689005530618575480,0.12310994469381424520},

233

{0.99056100455550913009,0.0094389954444908699092}},

234

V[2][2] = {{1,0},{0,1}},

235

S[2] = {1,1},

236

F[2] = {1,0},

237

Fembed[2] = {0,1},

238

Ferror[2] = {-1.0/(1.0-GAMMA0.25),1.0/(1.0-GAMMA0.25)},

239

Serror[2] = {1.0-GAMMA0.25,1.0};

240

ierr = TSGLEERegister(TSGLEE25I"25i",p,s,r,GAMMA0.25,&A[0][0],&B[0][0],&U[0][0],&V[0][0],S,F,NULL((void*)0),Fembed,Ferror,Serror,0,NULL((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),240,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

241

}

242

{

243

#undef GAMMA0.25

244

#define GAMMA0.25 0.0

245

/* y-y~ form */

246

const PetscInt

247

p = 3,

248

s = 5,

249

r = 2;

250

const PetscReal

251

A[5][5] = {{0,0,0,0,0},

252

{- 2169604947363702313.0 / 24313474998937147335.0,0,0,0,0},

253

{46526746497697123895.0 / 94116917485856474137.0,-10297879244026594958.0 / 49199457603717988219.0,0,0,0},

254

{23364788935845982499.0 / 87425311444725389446.0,-79205144337496116638.0 / 148994349441340815519.0,40051189859317443782.0 / 36487615018004984309.0,0,0},

255

{42089522664062539205.0 / 124911313006412840286.0,-15074384760342762939.0 / 137927286865289746282.0,-62274678522253371016.0 / 125918573676298591413.0,13755475729852471739.0 / 79257927066651693390.0,0}},

256

B[2][5] = {{61546696837458703723.0 / 56982519523786160813.0,-55810892792806293355.0 / 206957624151308356511.0,24061048952676379087.0 / 158739347956038723465.0,3577972206874351339.0 / 7599733370677197135.0,-59449832954780563947.0 / 137360038685338563670.0},

257

{- 9738262186984159168.0 / 99299082461487742983.0,-32797097931948613195.0 / 61521565616362163366.0,42895514606418420631.0 / 71714201188501437336.0,22608567633166065068.0 / 55371917805607957003.0,94655809487476459565.0 / 151517167160302729021.0}},

258

U[5][2] = {{70820309139834661559.0 / 80863923579509469826.0,10043614439674808267.0 / 80863923579509469826.0},

259

{161694774978034105510.0 / 106187653640211060371.0,-55507121337823045139.0 / 106187653640211060371.0},

260

{78486094644566264568.0 / 88171030896733822981.0,9684936252167558413.0 / 88171030896733822981.0},

261

{65394922146334854435.0 / 84570853840405479554.0,19175931694070625119.0 / 84570853840405479554.0},

262

{8607282770183754108.0 / 108658046436496925911.0,100050763666313171803.0 / 108658046436496925911.0}},

263

V[2][2] = {{1,0},{0,1}},

264

S[2] = {1,1},

265

F[2] = {1,0},

266

Fembed[2] = {0,1},

267

Ferror[2] = {-1.0/(1.0-GAMMA0.25),1.0/(1.0-GAMMA0.25)},

268

Serror[2] = {1.0-GAMMA0.25,1.0};

269

ierr = TSGLEERegister(TSGLEE35"35",p,s,r,GAMMA0.25,&A[0][0],&B[0][0],&U[0][0],&V[0][0],S,F,NULL((void*)0),Fembed,Ferror,Serror,0,NULL((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),269,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

270

}

271

{

272

#undef GAMMA0.25

273

#define GAMMA0.25 0.25

274

/* y-eps form */

275

const PetscInt

276

p = 2,

277

s = 6,

278

r = 2;

279

const PetscReal

280

A[6][6] = {{0,0,0,0,0,0},

281

{1,0,0,0,0,0},

282

{0,0,0,0,0,0},

283

{0,0,0.5,0,0,0},

284

{0,0,0.25,0.25,0,0},

285

{0,0,0.25,0.25,0.5,0}},

286

B[2][6] = {{0.5,0.5,0,0,0,0},

287

{-2.0/3.0,-2.0/3.0,1.0/3.0,1.0/3.0,1.0/3.0,1.0/3.0}},

288

U[6][2] = {{1,0},{1,0},{1,0.75},{1,0.75},{1,0.75},{1,0.75}},

289

V[2][2] = {{1,0},{0,1}},

290

S[2] = {1,0},

291

F[2] = {1,0},

292

Fembed[2] = {1,1-GAMMA0.25},

293

Ferror[2] = {0,1},

294

Serror[2] = {1,0};

295

ierr = TSGLEERegister(TSGLEEEXRK2A"exrk2a",p,s,r,GAMMA0.25,&A[0][0],&B[0][0],&U[0][0],&V[0][0],S,F,NULL((void*)0),Fembed,Ferror,Serror,0,NULL((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),295,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

296

}

297

{

298

#undef GAMMA0.25

299

#define GAMMA0.25 0.0

300

/* y-eps form */

301

const PetscInt

302

p = 3,

303

s = 8,

304

r = 2;

305

const PetscReal

306

A[8][8] = {{0,0,0,0,0,0,0,0},

307

{0.5,0,0,0,0,0,0,0},

308

{-1,2,0,0,0,0,0,0},

309

{1.0/6.0,2.0/3.0,1.0/6.0,0,0,0,0,0},

310

{0,0,0,0,0,0,0,0},

311

{-7.0/24.0,1.0/3.0,1.0/12.0,-1.0/8.0,0.5,0,0,0},

312

{7.0/6.0,-4.0/3.0,-1.0/3.0,0.5,-1.0,2.0,0,0},

313

{0,0,0,0,1.0/6.0,2.0/3.0,1.0/6.0,0}},

314

B[2][8] = {{1.0/6.0,2.0/3.0,1.0/6.0,0,0,0,0,0},

315

{-1.0/6.0,-2.0/3.0,-1.0/6.0,0,1.0/6.0,2.0/3.0,1.0/6.0,0}},

316

U[8][2] = {{1,0},{1,0},{1,0},{1,0},{1,1},{1,1},{1,1},{1,1}},

317

V[2][2] = {{1,0},{0,1}},

318

S[2] = {1,0},

319

F[2] = {1,0},

320

Fembed[2] = {1,1-GAMMA0.25},

321

Ferror[2] = {0,1},

322

Serror[2] = {1,0};

323

ierr = TSGLEERegister(TSGLEERK32G1"rk32g1",p,s,r,GAMMA0.25,&A[0][0],&B[0][0],&U[0][0],&V[0][0],S,F,NULL((void*)0),Fembed,Ferror,Serror,0,NULL((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),323,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

324

}

325

{

326

#undef GAMMA0.25

327

#define GAMMA0.25 0.25

328

/* y-eps form */

329

const PetscInt

330

p = 2,

331

s = 9,

332

r = 2;

333

const PetscReal

334

A[9][9] = {{0,0,0,0,0,0,0,0,0},

335

{0.585786437626904966,0,0,0,0,0,0,0,0},

336

{0.149999999999999994,0.849999999999999978,0,0,0,0,0,0,0},

337

{0,0,0,0,0,0,0,0,0},

338

{0,0,0,0.292893218813452483,0,0,0,0,0},

339

{0,0,0,0.0749999999999999972,0.424999999999999989,0,0,0,0},

340

{0,0,0,0.176776695296636893,0.176776695296636893,0.146446609406726241,0,0,0},

341

{0,0,0,0.176776695296636893,0.176776695296636893,0.146446609406726241,0.292893218813452483,0,0},

342

{0,0,0,0.176776695296636893,0.176776695296636893,0.146446609406726241,0.0749999999999999972,0.424999999999999989,0}},

343

B[2][9] = {{0.353553390593273786,0.353553390593273786,0.292893218813452483,0,0,0,0,0,0},

344

{-0.471404520791031678,-0.471404520791031678,-0.390524291751269959,0.235702260395515839,0.235702260395515839,0.195262145875634979,0.235702260395515839,0.235702260395515839,0.195262145875634979}},

345

U[9][2] = {{1,0},{1,0},{1,0},{1,0.75},{1,0.75},{1,0.75},{1,0.75},{1,0.75},{1,0.75}},

346

V[2][2] = {{1,0},{0,1}},

347

S[2] = {1,0},

348

F[2] = {1,0},

349

Fembed[2] = {1,1-GAMMA0.25},

350

Ferror[2] = {0,1},

351

Serror[2] = {1,0};

352

ierr = TSGLEERegister(TSGLEERK285EX"rk285ex",p,s,r,GAMMA0.25,&A[0][0],&B[0][0],&U[0][0],&V[0][0],S,F,NULL((void*)0),Fembed,Ferror,Serror,0,NULL((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),352,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

353

}

354

355

PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);

356

}

357

358

/*@C

359

TSGLEERegisterDestroy - Frees the list of schemes that were registered by TSGLEERegister().

360

361

Not Collective

362

363

Level: advanced

364

365

.seealso: TSGLEERegister(), TSGLEERegisterAll()

366

@*/

367

PetscErrorCode TSGLEERegisterDestroy(void)

368

{

369

PetscErrorCode ierr;

370

GLEETableauLink link;

371

372

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 372; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

373

while ((link = GLEETableauList)) {

374

GLEETableau t = &link->tab;

375

GLEETableauList = link->next;

376

ierr = PetscFree5(t->A,t->B,t->U,t->V,t->c)PetscFreeA(5,376,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,&(t->A),&(t->B),&(t->U),&(t->V),
&(t->c));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),376,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

377

ierr = PetscFree2(t->S,t->F)PetscFreeA(2,377,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,&(t->S),&(t->F)); CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),377,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

378

ierr = PetscFree (t->Fembed)((*PetscTrFree)((void*)(t->Fembed),378,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
) || ((t->Fembed) = 0,0)); CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),378,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

379

ierr = PetscFree (t->Ferror)((*PetscTrFree)((void*)(t->Ferror),379,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
) || ((t->Ferror) = 0,0)); CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),379,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

380

ierr = PetscFree (t->Serror)((*PetscTrFree)((void*)(t->Serror),380,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
) || ((t->Serror) = 0,0)); CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),380,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

381

ierr = PetscFree (t->binterp)((*PetscTrFree)((void*)(t->binterp),381,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
) || ((t->binterp) = 0,0)); CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),381,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

382

ierr = PetscFree (t->name)((*PetscTrFree)((void*)(t->name),382,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
) || ((t->name) = 0,0)); CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),382,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

383

ierr = PetscFree (link)((*PetscTrFree)((void*)(link),383,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
) || ((link) = 0,0)); CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),383,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

384

}

385

TSGLEERegisterAllCalled = PETSC_FALSE;

386

387

}

388

389

/*@C

390

TSGLEEInitializePackage - This function initializes everything in the TSGLEE package. It is called

391

from TSInitializePackage().

392

393

Level: developer

394

395

.seealso: PetscInitialize()

396

@*/

397

PetscErrorCode TSGLEEInitializePackage(void)

398

{

399

PetscErrorCode ierr;

400

401

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 401; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

402

if (TSGLEEPackageInitialized) PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);

403

TSGLEEPackageInitialized = PETSC_TRUE;

404

ierr = TSGLEERegisterAll();CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),404,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

405

ierr = PetscObjectComposedDataRegister(&explicit_stage_time_id);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),405,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

406

ierr = PetscRegisterFinalize(TSGLEEFinalizePackage);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),406,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

407

408

}

409

410

/*@C

411

TSGLEEFinalizePackage - This function destroys everything in the TSGLEE package. It is

412

called from PetscFinalize().

413

414

Level: developer

415

416

.seealso: PetscFinalize()

417

@*/

418

PetscErrorCode TSGLEEFinalizePackage(void)

419

{

420

PetscErrorCode ierr;

421

422

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 422; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

423

TSGLEEPackageInitialized = PETSC_FALSE;

424

ierr = TSGLEERegisterDestroy();CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),424,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

425

426

}

427

428

/*@C

429

TSGLEERegister - register an GLEE scheme by providing the entries in the Butcher tableau

430

431

Not Collective, but the same schemes should be registered on all processes on which they will be used

432

433

Input Parameters:

434

+ name - identifier for method

435

. order - order of method

436

. s - number of stages

437

. r - number of steps

438

. gamma - LTE ratio

439

. A - stage coefficients (dimension s*s, row-major)

440

. B - step completion coefficients (dimension r*s, row-major)

441

. U - method coefficients (dimension s*r, row-major)

442

. V - method coefficients (dimension r*r, row-major)

443

. S - starting coefficients

444

. F - finishing coefficients

445

. c - abscissa (dimension s; NULL to use row sums of A)

446

. Fembed - step completion coefficients for embedded method

447

. Ferror - error computation coefficients

448

. Serror - error initialization coefficients

449

. pinterp - order of interpolation (0 if unavailable)

450

- binterp - array of interpolation coefficients (NULL if unavailable)

451

452

Notes:

453

Several GLEE methods are provided, this function is only needed to create new methods.

454

455

Level: advanced

456

457

.seealso: TSGLEE

458

@*/

459

PetscErrorCode TSGLEERegister(TSGLEEType name,PetscInt order,PetscInt s, PetscInt r,

460

PetscReal gamma,

461

const PetscReal A[],const PetscReal B[],

462

const PetscReal U[],const PetscReal V[],

463

const PetscReal S[],const PetscReal F[],

464

const PetscReal c[],

465

const PetscReal Fembed[],const PetscReal Ferror[],

466

const PetscReal Serror[],

467

PetscInt pinterp, const PetscReal binterp[])

468

{

469

PetscErrorCode ierr;

470

GLEETableauLink link;

471

GLEETableau t;

472

PetscInt i,j;

473

474

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 474; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

475

ierr = TSGLEEInitializePackage();CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),475,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

476

ierr = PetscNew(&link)PetscMallocA(1,PETSC_TRUE,476,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,(size_t)(1)*sizeof(**((&link))),((&link)));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),476,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

477

t = &link->tab;

478

ierr = PetscStrallocpy(name,&t->name);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),478,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

479

t->order = order;

480

t->s = s;

481

t->r = r;

482

t->gamma = gamma;

483

ierr = PetscMalloc5(s*s,&t->A,r*r,&t->V,s,&t->c,r*s,&t->B,s*r,&t->U)PetscMallocA(5,PETSC_FALSE,483,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,(size_t)(s*s)*sizeof(**(&t->A)),(&t->A),(size_t
)(r*r)*sizeof(**(&t->V)),(&t->V),(size_t)(s)*sizeof
(**(&t->c)),(&t->c),(size_t)(r*s)*sizeof(**(&
t->B)),(&t->B),(size_t)(s*r)*sizeof(**(&t->U
)),(&t->U));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),483,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

484

ierr = PetscMalloc2(r,&t->S,r,&t->F)PetscMallocA(2,PETSC_FALSE,484,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,(size_t)(r)*sizeof(**(&t->S)),(&t->S),(size_t)
(r)*sizeof(**(&t->F)),(&t->F));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),484,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

485

ierr = PetscArraycpy(t->A,A,s*s)((sizeof(*(t->A)) != sizeof(*(A))) || PetscMemcpy(t->A,
A,(s*s)*sizeof(*(t->A))));;CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),485,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

486

ierr = PetscArraycpy(t->B,B,r*s)((sizeof(*(t->B)) != sizeof(*(B))) || PetscMemcpy(t->B,
B,(r*s)*sizeof(*(t->B))));;CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),486,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

487

ierr = PetscArraycpy(t->U,U,s*r)((sizeof(*(t->U)) != sizeof(*(U))) || PetscMemcpy(t->U,
U,(s*r)*sizeof(*(t->U))));;CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),487,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

488

ierr = PetscArraycpy(t->V,V,r*r)((sizeof(*(t->V)) != sizeof(*(V))) || PetscMemcpy(t->V,
V,(r*r)*sizeof(*(t->V))));;CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),488,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

489

ierr = PetscArraycpy(t->S,S,r )((sizeof(*(t->S)) != sizeof(*(S))) || PetscMemcpy(t->S,
S,(r)*sizeof(*(t->S))));;CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),489,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

490

ierr = PetscArraycpy(t->F,F,r )((sizeof(*(t->F)) != sizeof(*(F))) || PetscMemcpy(t->F,
F,(r)*sizeof(*(t->F))));;CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),490,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

491

if (c) {

492

ierr = PetscArraycpy(t->c,c,s)((sizeof(*(t->c)) != sizeof(*(c))) || PetscMemcpy(t->c,
c,(s)*sizeof(*(t->c))));;CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),492,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

493

} else {

494

for (i=0; i<s; i++) for (j=0,t->c[i]=0; j<s; j++) t->c[i] += A[i*s+j];

495

}

496

ierr = PetscMalloc1(r,&t->Fembed)PetscMallocA(1,PETSC_FALSE,496,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,(size_t)(r)*sizeof(**(&t->Fembed)),(&t->Fembed
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),496,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

497

ierr = PetscMalloc1(r,&t->Ferror)PetscMallocA(1,PETSC_FALSE,497,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,(size_t)(r)*sizeof(**(&t->Ferror)),(&t->Ferror
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),497,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

498

ierr = PetscMalloc1(r,&t->Serror)PetscMallocA(1,PETSC_FALSE,498,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,(size_t)(r)*sizeof(**(&t->Serror)),(&t->Serror
));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),498,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

499

ierr = PetscArraycpy(t->Fembed,Fembed,r)((sizeof(*(t->Fembed)) != sizeof(*(Fembed))) || PetscMemcpy
(t->Fembed,Fembed,(r)*sizeof(*(t->Fembed))));;CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),499,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

500

ierr = PetscArraycpy(t->Ferror,Ferror,r)((sizeof(*(t->Ferror)) != sizeof(*(Ferror))) || PetscMemcpy
(t->Ferror,Ferror,(r)*sizeof(*(t->Ferror))));;CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),500,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

501

ierr = PetscArraycpy(t->Serror,Serror,r)((sizeof(*(t->Serror)) != sizeof(*(Serror))) || PetscMemcpy
(t->Serror,Serror,(r)*sizeof(*(t->Serror))));;CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),501,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

502

t->pinterp = pinterp;

503

ierr = PetscMalloc1(s*pinterp,&t->binterp)PetscMallocA(1,PETSC_FALSE,503,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,(size_t)(s*pinterp)*sizeof(**(&t->binterp)),(&t->
binterp));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),503,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

504

ierr = PetscArraycpy(t->binterp,binterp,s*pinterp)((sizeof(*(t->binterp)) != sizeof(*(binterp))) || PetscMemcpy
(t->binterp,binterp,(s*pinterp)*sizeof(*(t->binterp))))
;;CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),504,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

505

506

link->next = GLEETableauList;

507

GLEETableauList = link;

508

509

}

510

511

static PetscErrorCode TSEvaluateStep_GLEE(TS ts,PetscInt order,Vec X,PetscBool *done)

512

{

513

TS_GLEE *glee = (TS_GLEE*) ts->data;

514

GLEETableau tab = glee->tableau;

515

PetscReal h, *B = tab->B, *V = tab->V,

516

*F = tab->F,

517

*Fembed = tab->Fembed;

518

PetscInt s = tab->s, r = tab->r, i, j;

519

Vec *Y = glee->Y, *YdotStage = glee->YdotStage;

520

PetscScalar *ws = glee->swork, *wr = glee->rwork;

521

PetscErrorCode ierr;

522

523

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 523; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

524

525

switch (glee->status) {

526

case TS_STEP_INCOMPLETE:

527

case TS_STEP_PENDING:

528

h = ts->time_step; break;

529

case TS_STEP_COMPLETE:

530

h = ts->ptime - ts->ptime_prev; break;

531

default: SETERRQ(PetscObjectComm((PetscObject)ts),PETSC_ERR_PLIB,"Invalid TSStepStatus")return PetscError(PetscObjectComm((PetscObject)ts),531,__func__
,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c",77,
PETSC_ERROR_INITIAL,"Invalid TSStepStatus");

532

}

533

534

if (order == tab->order) {

535

536

/* Note: Irrespective of whether status is TS_STEP_INCOMPLETE

537

or TS_STEP_COMPLETE, glee->X has the solution at the

538

beginning of the time step. So no need to roll-back.

539

540

if (glee->status == TS_STEP_INCOMPLETE) {

541

for (i=0; i<r; i++) {

542

ierr = VecZeroEntries(Y[i]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),542,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

543

for (j=0; j<r; j++) wr[j] = V[i*r+j];

544

ierr = VecMAXPY(Y[i],r,wr,glee->X);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),544,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

545

for (j=0; j<s; j++) ws[j] = h*B[i*s+j];

546

ierr = VecMAXPY(Y[i],s,ws,YdotStage);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),546,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

547

}

548

ierr = VecZeroEntries(X);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),548,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

549

for (j=0; j<r; j++) wr[j] = F[j];

550

ierr = VecMAXPY(X,r,wr,Y);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),550,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

551

} else {ierr = VecCopy(ts->vec_sol,X);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),551,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);}

552

553

554

} else if (order == tab->order-1) {

555

556

/* Complete with the embedded method (Fembed) */

557

for (i=0; i<r; i++) {

558

ierr = VecZeroEntries(Y[i]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),558,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

559

for (j=0; j<r; j++) wr[j] = V[i*r+j];

560

ierr = VecMAXPY(Y[i],r,wr,glee->X);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),560,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

561

for (j=0; j<s; j++) ws[j] = h*B[i*s+j];

562

ierr = VecMAXPY(Y[i],s,ws,YdotStage);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),562,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

563

}

564

ierr = VecZeroEntries(X);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),564,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

565

for (j=0; j<r; j++) wr[j] = Fembed[j];

566

ierr = VecMAXPY(X,r,wr,Y);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),566,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

567

568

if (done) *done = PETSC_TRUE;

569

570

}

571

if (done) *done = PETSC_FALSE;

572

else SETERRQ3(PetscObjectComm((PetscObject)ts),PETSC_ERR_SUP,"GLEE '%s' of order %D cannot evaluate step at order %D",tab->name,tab->order,order)return PetscError(PetscObjectComm((PetscObject)ts),572,__func__
,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c",56,
PETSC_ERROR_INITIAL,"GLEE '%s' of order %D cannot evaluate step at order %D"
,tab->name,tab->order,order);

573

574

}

575

576

static PetscErrorCode TSStep_GLEE(TS ts)

577

{

578

TS_GLEE *glee = (TS_GLEE*)ts->data;

579

GLEETableau tab = glee->tableau;

580

const PetscInt s = tab->s, r = tab->r;

581

PetscReal *A = tab->A, *U = tab->U,

582

*F = tab->F,

583

*c = tab->c;

584

Vec *Y = glee->Y, *X = glee->X,

585

*YStage = glee->YStage,

586

*YdotStage = glee->YdotStage,

587

W = glee->W;

588

SNES snes;

589

PetscScalar *ws = glee->swork, *wr = glee->rwork;

590

TSAdapt adapt;

591

PetscInt i,j,reject,next_scheme,its,lits;

592

PetscReal next_time_step;

593

PetscReal t;

594

PetscBool accept;

595

PetscErrorCode ierr;

596

597

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 597; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

598

ierr = PetscCitationsRegister(citation,&cited);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),598,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

599

600

for (i=0; i<r; i++) { ierr = VecCopy(Y[i],X[i]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),600,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0); }

601

602

ierr = TSGetSNES(ts,&snes);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),602,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

603

next_time_step = ts->time_step;

604

t = ts->ptime;

605

accept = PETSC_TRUE;

606

glee->status = TS_STEP_INCOMPLETE;

607

608

for (reject=0; reject<ts->max_reject && !ts->reason; reject++,ts->reject++) {

609

610

PetscReal h = ts->time_step;

611

ierr = TSPreStep(ts);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),611,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

612

613

for (i=0; i<s; i++) {

614

615

glee->stage_time = t + h*c[i];

616

ierr = TSPreStage(ts,glee->stage_time);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),616,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

617

618

if (A[i*s+i] == 0) { /* Explicit stage */

619

ierr = VecZeroEntries(YStage[i]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),619,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

620

for (j=0; j<r; j++) wr[j] = U[i*r+j];

621

ierr = VecMAXPY(YStage[i],r,wr,X);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),621,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

622

for (j=0; j<i; j++) ws[j] = h*A[i*s+j];

623

ierr = VecMAXPY(YStage[i],i,ws,YdotStage);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),623,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

624

} else { /* Implicit stage */

625

glee->scoeff = 1.0/A[i*s+i];

626

/* compute right-hand-side */

627

ierr = VecZeroEntries(W);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),627,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

628

for (j=0; j<r; j++) wr[j] = U[i*r+j];

629

ierr = VecMAXPY(W,r,wr,X);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),629,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

630

for (j=0; j<i; j++) ws[j] = h*A[i*s+j];

631

ierr = VecMAXPY(W,i,ws,YdotStage);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),631,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

632

ierr = VecScale(W,glee->scoeff/h);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),632,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

633

/* set initial guess */

634

ierr = VecCopy(i>0 ? YStage[i-1] : ts->vec_sol,YStage[i]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),634,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

635

/* solve for this stage */

636

ierr = SNESSolve(snes,W,YStage[i]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),636,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

637

ierr = SNESGetIterationNumber(snes,&its);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),637,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

638

ierr = SNESGetLinearSolveIterations(snes,&lits);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),638,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

639

ts->snes_its += its; ts->ksp_its += lits;

640

}

641

ierr = TSGetAdapt(ts,&adapt);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),641,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

642

ierr = TSAdaptCheckStage(adapt,ts,glee->stage_time,YStage[i],&accept);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),642,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

643

if (!accept) goto reject_step;

644

ierr = TSPostStage(ts,glee->stage_time,i,YStage);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),644,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

645

ierr = TSComputeRHSFunction(ts,t+h*c[i],YStage[i],YdotStage[i]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),645,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

646

}

647

ierr = TSEvaluateStep(ts,tab->order,ts->vec_sol,NULL((void*)0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),647,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

648

glee->status = TS_STEP_PENDING;

649

650

/* Register only the current method as a candidate because we're not supporting multiple candidates yet. */

651

ierr = TSGetAdapt(ts,&adapt);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),651,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

652

ierr = TSAdaptCandidatesClear(adapt);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),652,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

653

ierr = TSAdaptCandidateAdd(adapt,tab->name,tab->order,1,tab->ccfl,(PetscReal)tab->s,PETSC_TRUE);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),653,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

654

ierr = TSAdaptChoose(adapt,ts,ts->time_step,&next_scheme,&next_time_step,&accept);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),654,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

655

if (accept) {

656

/* ignore next_scheme for now */

657

ts->ptime += ts->time_step;

658

ts->time_step = next_time_step;

659

glee->status = TS_STEP_COMPLETE;

660

/* compute and store the global error */

661

/* Note: this is not needed if TSAdaptGLEE is not used */

662

ierr = TSGetTimeError(ts,0,&(glee->yGErr));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),662,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

663

ierr = PetscObjectComposedDataSetReal((PetscObject)ts->vec_sol,explicit_stage_time_id,ts->ptime)(((((PetscObject)ts->vec_sol)->real_idmax < PetscObjectComposedDataMax
) && PetscObjectComposedDataIncreaseReal((PetscObject
)ts->vec_sol)) || (((PetscObject)ts->vec_sol)->realcomposeddata
[explicit_stage_time_id] = ts->ptime,((PetscObject)ts->
vec_sol)->realcomposedstate[explicit_stage_time_id] = ((PetscObject
)ts->vec_sol)->state, 0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),663,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

664

break;

665

} else { /* Roll back the current step */

666

for (j=0; j<r; j++) wr[j] = F[j];

667

ierr = VecMAXPY(ts->vec_sol,r,wr,X);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),667,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

668

ts->time_step = next_time_step;

669

glee->status = TS_STEP_INCOMPLETE;

670

}

671

reject_step: continue;

672

}

673

if (glee->status != TS_STEP_COMPLETE && !ts->reason) ts->reason = TS_DIVERGED_STEP_REJECTED;

674

675

}

676

677

static PetscErrorCode TSInterpolate_GLEE(TS ts,PetscReal itime,Vec X)

678

{

679

TS_GLEE *glee = (TS_GLEE*)ts->data;

680

PetscInt s=glee->tableau->s, pinterp=glee->tableau->pinterp,i,j;

681

PetscReal h,tt,t;

682

PetscScalar *b;

683

const PetscReal *B = glee->tableau->binterp;

684

PetscErrorCode ierr;

685

686

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 686; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

687

if (!B) SETERRQ1(PetscObjectComm((PetscObject)ts),PETSC_ERR_SUP,"TSGLEE %s does not have an interpolation formula",glee->tableau->name)return PetscError(PetscObjectComm((PetscObject)ts),687,__func__
,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c",56,
PETSC_ERROR_INITIAL,"TSGLEE %s does not have an interpolation formula"
,glee->tableau->name);

688

switch (glee->status) {

689

case TS_STEP_INCOMPLETE:

690

case TS_STEP_PENDING:

691

h = ts->time_step;

692

t = (itime - ts->ptime)/h;

693

break;

694

case TS_STEP_COMPLETE:

695

h = ts->ptime - ts->ptime_prev;

696

t = (itime - ts->ptime)/h + 1; /* In the interval [0,1] */

697

break;

698

default: SETERRQ(PetscObjectComm((PetscObject)ts),PETSC_ERR_PLIB,"Invalid TSStepStatus")return PetscError(PetscObjectComm((PetscObject)ts),698,__func__
,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c",77,
PETSC_ERROR_INITIAL,"Invalid TSStepStatus");

699

}

700

ierr = PetscMalloc1(s,&b)PetscMallocA(1,PETSC_FALSE,700,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,(size_t)(s)*sizeof(**(&b)),(&b));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),700,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

701

for (i=0; i<s; i++) b[i] = 0;

702

for (j=0,tt=t; j<pinterp; j++,tt*=t) {

703

for (i=0; i<s; i++) {

704

b[i] += h * B[i*pinterp+j] * tt;

705

}

706

}

707

ierr = VecCopy(glee->YStage[0],X);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),707,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

708

ierr = VecMAXPY(X,s,b,glee->YdotStage);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),708,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

709

ierr = PetscFree(b)((*PetscTrFree)((void*)(b),709,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
) || ((b) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),709,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

710

711

}

712

713

/*------------------------------------------------------------*/

714

static PetscErrorCode TSReset_GLEE(TS ts)

715

{

716

TS_GLEE *glee = (TS_GLEE*)ts->data;

717

PetscInt s, r;

718

PetscErrorCode ierr;

719

720

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 720; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

721

if (!glee->tableau) PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);

722

s = glee->tableau->s;

723

r = glee->tableau->r;

724

ierr = VecDestroyVecs(r,&glee->Y);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),724,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

725

ierr = VecDestroyVecs(r,&glee->X);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),725,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

726

ierr = VecDestroyVecs(s,&glee->YStage);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),726,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

727

ierr = VecDestroyVecs(s,&glee->YdotStage);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),727,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

728

ierr = VecDestroy(&glee->Ydot);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),728,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

729

ierr = VecDestroy(&glee->yGErr);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),729,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

730

ierr = VecDestroy(&glee->W);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),730,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

731

ierr = PetscFree2(glee->swork,glee->rwork)PetscFreeA(2,731,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,&(glee->swork),&(glee->rwork));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),731,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

732

733

}

734

735

static PetscErrorCode TSGLEEGetVecs(TS ts,DM dm,Vec *Ydot)

736

{

737

TS_GLEE *glee = (TS_GLEE*)ts->data;

738

PetscErrorCode ierr;

739

740

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 740; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

741

if (Ydot) {

742

if (dm && dm != ts->dm) {

743

ierr = DMGetNamedGlobalVector(dm,"TSGLEE_Ydot",Ydot);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),743,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

744

} else *Ydot = glee->Ydot;

745

}

746

747

}

748

749

750

static PetscErrorCode TSGLEERestoreVecs(TS ts,DM dm,Vec *Ydot)

751

{

752

PetscErrorCode ierr;

753

754

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 754; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

755

if (Ydot) {

756

if (dm && dm != ts->dm) {

757

ierr = DMRestoreNamedGlobalVector(dm,"TSGLEE_Ydot",Ydot);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),757,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

758

}

759

}

760

761

}

762

763

764

This defines the nonlinear equation that is to be solved with SNES

765

766

static PetscErrorCode SNESTSFormFunction_GLEE(SNES snes,Vec X,Vec F,TS ts)

767

{

768

TS_GLEE *glee = (TS_GLEE*)ts->data;

769

DM dm,dmsave;

770

Vec Ydot;

771

PetscReal shift = glee->scoeff / ts->time_step;

772

PetscErrorCode ierr;

773

774

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 774; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

775

ierr = SNESGetDM(snes,&dm);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),775,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

776

ierr = TSGLEEGetVecs(ts,dm,&Ydot);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),776,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

777

/* Set Ydot = shift*X */

778

ierr = VecCopy(X,Ydot);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),778,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

779

ierr = VecScale(Ydot,shift);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),779,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

780

dmsave = ts->dm;

781

ts->dm = dm;

782

783

ierr = TSComputeIFunction(ts,glee->stage_time,X,Ydot,F,PETSC_FALSE);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),783,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

784

785

ts->dm = dmsave;

786

ierr = TSGLEERestoreVecs(ts,dm,&Ydot);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),786,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

787

788

}

789

790

static PetscErrorCode SNESTSFormJacobian_GLEE(SNES snes,Vec X,Mat A,Mat B,TS ts)

791

{

792

TS_GLEE *glee = (TS_GLEE*)ts->data;

793

DM dm,dmsave;

794

Vec Ydot;

795

PetscReal shift = glee->scoeff / ts->time_step;

796

PetscErrorCode ierr;

797

798

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 798; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

799

ierr = SNESGetDM(snes,&dm);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),799,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

800

ierr = TSGLEEGetVecs(ts,dm,&Ydot);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),800,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

801

/* glee->Ydot has already been computed in SNESTSFormFunction_GLEE (SNES guarantees this) */

802

dmsave = ts->dm;

803

ts->dm = dm;

804

805

ierr = TSComputeIJacobian(ts,glee->stage_time,X,Ydot,shift,A,B,PETSC_FALSE);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),805,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

806

807

ts->dm = dmsave;

808

ierr = TSGLEERestoreVecs(ts,dm,&Ydot);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),808,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

809

810

}

811

812

static PetscErrorCode DMCoarsenHook_TSGLEE(DM fine,DM coarse,void *ctx)

813

{

814

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 814; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

815

816

}

817

818

static PetscErrorCode DMRestrictHook_TSGLEE(DM fine,Mat restrct,Vec rscale,Mat inject,DM coarse,void *ctx)

819

{

820

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 820; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

821

822

}

823

824

825

static PetscErrorCode DMSubDomainHook_TSGLEE(DM dm,DM subdm,void *ctx)

826

{

827

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 827; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

828

829

}

830

831

static PetscErrorCode DMSubDomainRestrictHook_TSGLEE(DM dm,VecScatter gscat,VecScatter lscat,DM subdm,void *ctx)

832

{

833

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 833; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

834

835

}

836

837

static PetscErrorCode TSSetUp_GLEE(TS ts)

838

{

839

TS_GLEE *glee = (TS_GLEE*)ts->data;

840

GLEETableau tab;

841

PetscInt s,r;

842

PetscErrorCode ierr;

843

DM dm;

844

845

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 845; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

846

if (!glee->tableau) {

Assuming pointer value is null

→

←

Taking true branch

→

847

ierr = TSGLEESetType(ts,TSGLEEDefaultType);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),847,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

848

}

849

tab = glee->tableau;

←

Null pointer value stored to 'tab'

→

850

s = tab->s;

←

Access to field 's' results in a dereference of a null pointer (loaded from variable 'tab')

851

r = tab->r;

852

ierr = VecDuplicateVecs(ts->vec_sol,r,&glee->Y);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),852,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

853

ierr = VecDuplicateVecs(ts->vec_sol,r,&glee->X);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),853,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

854

ierr = VecDuplicateVecs(ts->vec_sol,s,&glee->YStage);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),854,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

855

ierr = VecDuplicateVecs(ts->vec_sol,s,&glee->YdotStage);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),855,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

856

ierr = VecDuplicate(ts->vec_sol,&glee->Ydot);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),856,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

857

ierr = VecDuplicate(ts->vec_sol,&glee->yGErr);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),857,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

858

ierr = VecZeroEntries(glee->yGErr);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),858,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

859

ierr = VecDuplicate(ts->vec_sol,&glee->W);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),859,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

860

ierr = PetscMalloc2(s,&glee->swork,r,&glee->rwork)PetscMallocA(2,PETSC_FALSE,860,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,(size_t)(s)*sizeof(**(&glee->swork)),(&glee->swork
),(size_t)(r)*sizeof(**(&glee->rwork)),(&glee->
rwork));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),860,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

861

ierr = TSGetDM(ts,&dm);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),861,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

862

ierr = DMCoarsenHookAdd(dm,DMCoarsenHook_TSGLEE,DMRestrictHook_TSGLEE,ts);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),862,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

863

ierr = DMSubDomainHookAdd(dm,DMSubDomainHook_TSGLEE,DMSubDomainRestrictHook_TSGLEE,ts);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),863,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

864

865

}

866

867

PetscErrorCode TSStartingMethod_GLEE(TS ts)

868

{

869

TS_GLEE *glee = (TS_GLEE*)ts->data;

870

GLEETableau tab = glee->tableau;

871

PetscInt r=tab->r,i;

872

PetscReal *S=tab->S;

873

PetscErrorCode ierr;

874

875

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 875; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

876

for (i=0; i<r; i++) {

877

ierr = VecZeroEntries(glee->Y[i]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),877,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

878

ierr = VecAXPY(glee->Y[i],S[i],ts->vec_sol);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),878,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

879

}

880

881

882

}

883

884

885

/*------------------------------------------------------------*/

886

887

static PetscErrorCode TSSetFromOptions_GLEE(PetscOptionItems *PetscOptionsObject,TS ts)

888

{

889

PetscErrorCode ierr;

890

char gleetype[256];

891

892

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 892; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

893

ierr = PetscOptionsHead(PetscOptionsObject,"GLEE ODE solver options");CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),893,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

894

{

895

GLEETableauLink link;

896

PetscInt count,choice;

897

PetscBool flg;

898

const char **namelist;

899

900

ierr = PetscStrncpy(gleetype,TSGLEEDefaultType,sizeof(gleetype));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),900,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

901

for (link=GLEETableauList,count=0; link; link=link->next,count++) ;

902

ierr = PetscMalloc1(count,(char***)&namelist)PetscMallocA(1,PETSC_FALSE,902,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,(size_t)(count)*sizeof(**((char***)&namelist)),((char***
)&namelist));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),902,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

903

for (link=GLEETableauList,count=0; link; link=link->next,count++) namelist[count] = link->tab.name;

904

ierr = PetscOptionsEList("-ts_glee_type","Family of GLEE method","TSGLEESetType",(const char*const*)namelist,count,gleetype,&choice,&flg)PetscOptionsEList_Private(PetscOptionsObject,"-ts_glee_type",
"Family of GLEE method","TSGLEESetType",(const char*const*)namelist
,count,gleetype,&choice,&flg);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),904,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

905

ierr = TSGLEESetType(ts,flg ? namelist[choice] : gleetype);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),905,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

906

ierr = PetscFree(namelist)((*PetscTrFree)((void*)(namelist),906,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
) || ((namelist) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),906,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

907

}

908

ierr = PetscOptionsTail()0; do {if (PetscOptionsObject->count != 1) do { do { ; if (
petscstack && petscstack->currentsize > 0) { petscstack
->currentsize--; petscstack->function[petscstack->currentsize
] = 0; petscstack->file[petscstack->currentsize] = 0; petscstack
->line[petscstack->currentsize] = 0; petscstack->petscroutine
[petscstack->currentsize] = PETSC_FALSE; } if (petscstack)
{ petscstack->hotdepth = (((petscstack->hotdepth-1)<
(0)) ? (0) : (petscstack->hotdepth-1)); } ; } while (0); return
(0);} while (0);} while(0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),908,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

909

910

}

911

912

static PetscErrorCode TSView_GLEE(TS ts,PetscViewer viewer)

913

{

914

TS_GLEE *glee = (TS_GLEE*)ts->data;

915

GLEETableau tab = glee->tableau;

916

PetscBool iascii;

917

PetscErrorCode ierr;

918

919

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 919; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

920

ierr = PetscObjectTypeCompare((PetscObject)viewer,PETSCVIEWERASCII"ascii",&iascii);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),920,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

921

if (iascii) {

922

TSGLEEType gleetype;

923

char buf[512];

924

ierr = TSGLEEGetType(ts,&gleetype);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),924,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

925

ierr = PetscViewerASCIIPrintf(viewer," GLEE type %s\n",gleetype);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),925,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

926

ierr = PetscFormatRealArray(buf,sizeof(buf),"% 8.6f",tab->s,tab->c);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),926,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

927

ierr = PetscViewerASCIIPrintf(viewer," Abscissa c = %s\n",buf);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),927,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

928

/* Note: print out r as well */

929

}

930

931

}

932

933

static PetscErrorCode TSLoad_GLEE(TS ts,PetscViewer viewer)

934

{

935

PetscErrorCode ierr;

936

SNES snes;

937

TSAdapt tsadapt;

938

939

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 939; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

940

ierr = TSGetAdapt(ts,&tsadapt);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),940,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

941

ierr = TSAdaptLoad(tsadapt,viewer);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),941,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

942

ierr = TSGetSNES(ts,&snes);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),942,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

943

ierr = SNESLoad(snes,viewer);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),943,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

944

/* function and Jacobian context for SNES when used with TS is always ts object */

945

ierr = SNESSetFunction(snes,NULL((void*)0),NULL((void*)0),ts);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),945,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

946

ierr = SNESSetJacobian(snes,NULL((void*)0),NULL((void*)0),NULL((void*)0),ts);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),946,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

947

948

}

949

950

/*@C

951

TSGLEESetType - Set the type of GLEE scheme

952

953

Logically collective

954

955

Input Parameter:

956

+ ts - timestepping context

957

- gleetype - type of GLEE-scheme

958

959

Level: intermediate

960

961

.seealso: TSGLEEGetType(), TSGLEE

962

@*/

963

PetscErrorCode TSGLEESetType(TS ts,TSGLEEType gleetype)

964

{

965

PetscErrorCode ierr;

966

967

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 967; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

968

PetscValidHeaderSpecific(ts,TS_CLASSID,1)do { if (!ts) return PetscError(((MPI_Comm)0x44000001),968,__func__
,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c",85,
PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",1); if (!PetscCheckPointer
(ts,PETSC_OBJECT)) return PetscError(((MPI_Comm)0x44000001),968
,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,1); if (((PetscObject)(ts))->classid != TS_CLASSID) { if (
((PetscObject)(ts))->classid == -1) return PetscError(((MPI_Comm
)0x44000001),968,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,1); else return PetscError(((MPI_Comm)0x44000001),968,__func__
,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c",62,
PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d",1)
; } } while (0);

969

PetscValidCharPointer(gleetype,2)do { if (!gleetype) return PetscError(((MPI_Comm)0x44000001),
969,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,85,PETSC_ERROR_INITIAL,"Null Pointer: Parameter # %d",2); if
(!PetscCheckPointer(gleetype,PETSC_CHAR)) return PetscError(
((MPI_Comm)0x44000001),969,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,68,PETSC_ERROR_INITIAL,"Invalid Pointer to char: Parameter # %d"
,2); } while (0);

970

ierr = PetscTryMethod(ts,"TSGLEESetType_C",(TS,TSGLEEType),(ts,gleetype))0; do { PetscErrorCode (*f)(TS,TSGLEEType), __ierr; __ierr = PetscObjectQueryFunction_Private
(((PetscObject)ts),("TSGLEESetType_C"),(PetscVoidFunction*)(&
f));do {if (__builtin_expect(!!(__ierr),0)) return PetscError
(((MPI_Comm)0x44000001),970,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0); if (f) {__ierr =
(*f)(ts,gleetype);do {if (__builtin_expect(!!(__ierr),0)) return
PetscError(((MPI_Comm)0x44000001),970,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0);} } while(0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),970,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

971

972

}

973

974

/*@C

975

TSGLEEGetType - Get the type of GLEE scheme

976

977

Logically collective

978

979

Input Parameter:

980

. ts - timestepping context

981

982

Output Parameter:

983

. gleetype - type of GLEE-scheme

984

985

Level: intermediate

986

987

.seealso: TSGLEESetType()

988

@*/

989

PetscErrorCode TSGLEEGetType(TS ts,TSGLEEType *gleetype)

990

{

991

PetscErrorCode ierr;

992

993

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 993; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

994

PetscValidHeaderSpecific(ts,TS_CLASSID,1)do { if (!ts) return PetscError(((MPI_Comm)0x44000001),994,__func__
,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c",85,
PETSC_ERROR_INITIAL,"Null Object: Parameter # %d",1); if (!PetscCheckPointer
(ts,PETSC_OBJECT)) return PetscError(((MPI_Comm)0x44000001),994
,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,64,PETSC_ERROR_INITIAL,"Invalid Pointer to Object: Parameter # %d"
,1); if (((PetscObject)(ts))->classid != TS_CLASSID) { if (
((PetscObject)(ts))->classid == -1) return PetscError(((MPI_Comm
)0x44000001),994,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,64,PETSC_ERROR_INITIAL,"Object already free: Parameter # %d"
,1); else return PetscError(((MPI_Comm)0x44000001),994,__func__
,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c",62,
PETSC_ERROR_INITIAL,"Wrong type of object: Parameter # %d",1)
; } } while (0);

995

ierr = PetscUseMethod(ts,"TSGLEEGetType_C",(TS,TSGLEEType*),(ts,gleetype))0; do { PetscErrorCode (*f)(TS,TSGLEEType*), __ierr; __ierr =
PetscObjectQueryFunction_Private(((PetscObject)ts),("TSGLEEGetType_C"
),(PetscVoidFunction*)(&f));do {if (__builtin_expect(!!(__ierr
),0)) return PetscError(((MPI_Comm)0x44000001),995,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0); if (f) {__ierr =
(*f)(ts,gleetype);do {if (__builtin_expect(!!(__ierr),0)) return
PetscError(((MPI_Comm)0x44000001),995,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0);} else return PetscError
(PetscObjectComm((PetscObject)ts),995,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,56,PETSC_ERROR_INITIAL,"Cannot locate function %s in object"
,"TSGLEEGetType_C"); } while(0);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),995,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

996

997

}

998

999

PetscErrorCode TSGLEEGetType_GLEE(TS ts,TSGLEEType *gleetype)

1000

{

1001

TS_GLEE *glee = (TS_GLEE*)ts->data;

1002

PetscErrorCode ierr;

1003

1004

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 1004; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1005

if (!glee->tableau) {

1006

ierr = TSGLEESetType(ts,TSGLEEDefaultType);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1006,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1007

}

1008

*gleetype = glee->tableau->name;

1009

1010

}

1011

PetscErrorCode TSGLEESetType_GLEE(TS ts,TSGLEEType gleetype)

1012

{

1013

TS_GLEE *glee = (TS_GLEE*)ts->data;

1014

PetscErrorCode ierr;

1015

PetscBool match;

1016

GLEETableauLink link;

1017

1018

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 1018; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1019

if (glee->tableau) {

1020

ierr = PetscStrcmp(glee->tableau->name,gleetype,&match);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1020,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1021

if (match) PetscFunctionReturn(0)do { do { ; if (petscstack && petscstack->currentsize
> 0) { petscstack->currentsize--; petscstack->function
[petscstack->currentsize] = 0; petscstack->file[petscstack
->currentsize] = 0; petscstack->line[petscstack->currentsize
] = 0; petscstack->petscroutine[petscstack->currentsize
] = PETSC_FALSE; } if (petscstack) { petscstack->hotdepth =
(((petscstack->hotdepth-1)<(0)) ? (0) : (petscstack->
hotdepth-1)); } ; } while (0); return(0);} while (0);

1022

}

1023

for (link = GLEETableauList; link; link=link->next) {

1024

ierr = PetscStrcmp(link->tab.name,gleetype,&match);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1024,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1025

if (match) {

1026

ierr = TSReset_GLEE(ts);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1026,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1027

glee->tableau = &link->tab;

1028

1029

}

1030

}

1031

SETERRQ1(PetscObjectComm((PetscObject)ts),PETSC_ERR_ARG_UNKNOWN_TYPE,"Could not find '%s'",gleetype)return PetscError(PetscObjectComm((PetscObject)ts),1031,__func__
,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c",86,
PETSC_ERROR_INITIAL,"Could not find '%s'",gleetype);

1032

1033

}

1034

1035

static PetscErrorCode TSGetStages_GLEE(TS ts,PetscInt *ns,Vec **Y)

1036

{

1037

TS_GLEE *glee = (TS_GLEE*)ts->data;

1038

1039

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 1039; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1040

if (ns) *ns = glee->tableau->s;

1041

if (Y) *Y = glee->YStage;

1042

1043

}

1044

1045

PetscErrorCode TSGetSolutionComponents_GLEE(TS ts,PetscInt *n,Vec *Y)

1046

{

1047

TS_GLEE *glee = (TS_GLEE*)ts->data;

1048

GLEETableau tab = glee->tableau;

1049

PetscErrorCode ierr;

1050

1051

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 1051; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1052

if (!Y) *n = tab->r;

1053

else {

1054

if ((*n >= 0) && (*n < tab->r)) {

1055

ierr = VecCopy(glee->Y[*n],*Y);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1055,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1056

} else SETERRQ3(PetscObjectComm((PetscObject)ts),PETSC_ERR_ARG_OUTOFRANGE,"Second argument (%d) out of range[%d,%d].",*n,0,tab->r-1)return PetscError(PetscObjectComm((PetscObject)ts),1056,__func__
,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c",63,
PETSC_ERROR_INITIAL,"Second argument (%d) out of range[%d,%d]."
,*n,0,tab->r-1);

1057

}

1058

1059

}

1060

1061

PetscErrorCode TSGetAuxSolution_GLEE(TS ts,Vec *X)

1062

{

1063

TS_GLEE *glee = (TS_GLEE*)ts->data;

1064

GLEETableau tab = glee->tableau;

1065

PetscReal *F = tab->Fembed;

1066

PetscInt r = tab->r;

1067

Vec *Y = glee->Y;

1068

PetscScalar *wr = glee->rwork;

1069

PetscInt i;

1070

PetscErrorCode ierr;

1071

1072

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 1072; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1073

ierr = VecZeroEntries(*X);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1073,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1074

for (i=0; i<r; i++) wr[i] = F[i];

1075

ierr = VecMAXPY((*X),r,wr,Y);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1075,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1076

1077

}

1078

1079

PetscErrorCode TSGetTimeError_GLEE(TS ts,PetscInt n,Vec *X)

1080

{

1081

TS_GLEE *glee = (TS_GLEE*)ts->data;

1082

GLEETableau tab = glee->tableau;

1083

PetscReal *F = tab->Ferror;

1084

PetscInt r = tab->r;

1085

Vec *Y = glee->Y;

1086

PetscScalar *wr = glee->rwork;

1087

PetscInt i;

1088

PetscErrorCode ierr;

1089

1090

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 1090; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1091

ierr = VecZeroEntries(*X);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1091,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1092

if(n==0){

1093

for (i=0; i<r; i++) wr[i] = F[i];

1094

ierr = VecMAXPY((*X),r,wr,Y);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1094,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1095

} else if(n==-1) {

1096

*X=glee->yGErr;

1097

}

1098

1099

}

1100

1101

PetscErrorCode TSSetTimeError_GLEE(TS ts,Vec X)

1102

{

1103

TS_GLEE *glee = (TS_GLEE*)ts->data;

1104

GLEETableau tab = glee->tableau;

1105

PetscReal *S = tab->Serror;

1106

PetscInt r = tab->r,i;

1107

Vec *Y = glee->Y;

1108

PetscErrorCode ierr;

1109

1110

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 1110; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1111

if (r != 2) SETERRQ2(PetscObjectComm((PetscObject)ts),PETSC_ERR_SUP,"TSSetTimeError_GLEE not supported for '%s' with r=%D.",tab->name,tab->r)return PetscError(PetscObjectComm((PetscObject)ts),1111,__func__
,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c",56,
PETSC_ERROR_INITIAL,"TSSetTimeError_GLEE not supported for '%s' with r=%D."
,tab->name,tab->r);

1112

for (i=1; i<r; i++) {

1113

ierr = VecCopy(ts->vec_sol,Y[i]);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1113,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1114

ierr = VecAXPBY(Y[i],S[0],S[1],X);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1114,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1115

ierr = VecCopy(X,glee->yGErr);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1115,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1116

}

1117

1118

}

1119

1120

static PetscErrorCode TSDestroy_GLEE(TS ts)

1121

{

1122

PetscErrorCode ierr;

1123

1124

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 1124; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1125

ierr = TSReset_GLEE(ts);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1125,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1126

if (ts->dm) {

1127

ierr = DMCoarsenHookRemove(ts->dm,DMCoarsenHook_TSGLEE,DMRestrictHook_TSGLEE,ts);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1127,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1128

ierr = DMSubDomainHookRemove(ts->dm,DMSubDomainHook_TSGLEE,DMSubDomainRestrictHook_TSGLEE,ts);CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1128,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1129

}

1130

ierr = PetscFree(ts->data)((*PetscTrFree)((void*)(ts->data),1130,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
) || ((ts->data) = 0,0));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1130,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1131

ierr = PetscObjectComposeFunction((PetscObject)ts,"TSGLEEGetType_C",NULL)PetscObjectComposeFunction_Private((PetscObject)ts,"TSGLEEGetType_C"
,(PetscVoidFunction)(((void*)0)));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1131,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1132

ierr = PetscObjectComposeFunction((PetscObject)ts,"TSGLEESetType_C",NULL)PetscObjectComposeFunction_Private((PetscObject)ts,"TSGLEESetType_C"
,(PetscVoidFunction)(((void*)0)));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1132,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1133

1134

}

1135

1136

/* ------------------------------------------------------------ */

1137

/*MC

1138

TSGLEE - ODE and DAE solver using General Linear with Error Estimation schemes

1139

1140

The user should provide the right hand side of the equation

1141

using TSSetRHSFunction().

1142

1143

Notes:

1144

The default is TSGLEE35, it can be changed with TSGLEESetType() or -ts_glee_type

1145

1146

Level: beginner

1147

1148

.seealso: TSCreate(), TS, TSSetType(), TSGLEESetType(), TSGLEEGetType(),

1149

TSGLEE23, TTSGLEE24, TSGLEE35, TSGLEE25I, TSGLEEEXRK2A,

1150

TSGLEERK32G1, TSGLEERK285EX, TSGLEEType, TSGLEERegister()

1151

1152

M*/

1153

PETSC_EXTERNextern __attribute__((visibility ("default"))) PetscErrorCode TSCreate_GLEE(TS ts)

1154

{

1155

TS_GLEE *th;

1156

PetscErrorCode ierr;

1157

1158

PetscFunctionBegindo { do { ; if (petscstack && (petscstack->currentsize
< 64)) { petscstack->function[petscstack->currentsize
] = __func__; petscstack->file[petscstack->currentsize]
= "/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"; petscstack
->line[petscstack->currentsize] = 1158; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1159

#if !defined(PETSC_USE_DYNAMIC_LIBRARIES)

1160

ierr = TSGLEEInitializePackage();CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1160,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1161

#endif

1162

1163

ts->ops->reset = TSReset_GLEE;

1164

ts->ops->destroy = TSDestroy_GLEE;

1165

ts->ops->view = TSView_GLEE;

1166

ts->ops->load = TSLoad_GLEE;

1167

ts->ops->setup = TSSetUp_GLEE;

1168

ts->ops->step = TSStep_GLEE;

1169

ts->ops->interpolate = TSInterpolate_GLEE;

1170

ts->ops->evaluatestep = TSEvaluateStep_GLEE;

1171

ts->ops->setfromoptions = TSSetFromOptions_GLEE;

1172

ts->ops->getstages = TSGetStages_GLEE;

1173

ts->ops->snesfunction = SNESTSFormFunction_GLEE;

1174

ts->ops->snesjacobian = SNESTSFormJacobian_GLEE;

1175

ts->ops->getsolutioncomponents = TSGetSolutionComponents_GLEE;

1176

ts->ops->getauxsolution = TSGetAuxSolution_GLEE;

1177

ts->ops->gettimeerror = TSGetTimeError_GLEE;

1178

ts->ops->settimeerror = TSSetTimeError_GLEE;

1179

ts->ops->startingmethod = TSStartingMethod_GLEE;

1180

ts->default_adapt_type = TSADAPTGLEE"glee";

1181

1182

ts->usessnes = PETSC_TRUE;

1183

1184

ierr = PetscNewLog(ts,&th)(PetscMallocA(1,PETSC_TRUE,1184,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,(size_t)(1)*sizeof(**(((&th)))),(((&th)))) || PetscLogObjectMemory
((PetscObject)ts,sizeof(**(&th))));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1184,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1185

ts->data = (void*)th;

1186

1187

ierr = PetscObjectComposeFunction((PetscObject)ts,"TSGLEEGetType_C",TSGLEEGetType_GLEE)PetscObjectComposeFunction_Private((PetscObject)ts,"TSGLEEGetType_C"
,(PetscVoidFunction)(TSGLEEGetType_GLEE));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1187,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1188

ierr = PetscObjectComposeFunction((PetscObject)ts,"TSGLEESetType_C",TSGLEESetType_GLEE)PetscObjectComposeFunction_Private((PetscObject)ts,"TSGLEESetType_C"
,(PetscVoidFunction)(TSGLEESetType_GLEE));CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1188,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1189

1190

}