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

Bug Summary

File:	ts/impls/glee/glee.c
Warning:	line 858, 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

.keywords: TS, TSGLEE, register, all

145

146

.seealso: TSGLEERegisterDestroy()

147

@*/

148

PetscErrorCode TSGLEERegisterAll(void)

149

{

150

PetscErrorCode ierr;

151

152

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] = 152; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

153

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);

154

TSGLEERegisterAllCalled = PETSC_TRUE;

155

156

{

157

#define GAMMA0.25 0.5

158

/* y-eps form */

159

const PetscInt

160

p = 1,

161

s = 3,

162

r = 2;

163

const PetscReal

164

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

165

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

166

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

167

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

168

S[2] = {1,0},

169

F[2] = {1,0},

170

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

171

Ferror[2] = {0,1},

172

Serror[2] = {1,0};

173

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),173,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

174

}

175

{

176

#undef GAMMA0.25

177

#define GAMMA0.25 0.0

178

/* y-eps form */

179

const PetscInt

180

p = 2,

181

s = 3,

182

r = 2;

183

const PetscReal

184

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

185

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}},

186

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

187

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

188

S[2] = {1,0},

189

F[2] = {1,0},

190

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

191

Ferror[2] = {0,1},

192

Serror[2] = {1,0};

193

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),193,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

194

}

195

{

196

#undef GAMMA0.25

197

#define GAMMA0.25 0.0

198

/* y-y~ form */

199

const PetscInt

200

p = 2,

201

s = 4,

202

r = 2;

203

const PetscReal

204

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}},

205

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}},

206

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

207

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

208

S[2] = {1,1},

209

F[2] = {1,0},

210

Fembed[2] = {0,1},

211

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

212

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

213

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),213,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

214

}

215

{

216

#undef GAMMA0.25

217

#define GAMMA0.25 0.0

218

/* y-y~ form */

219

const PetscInt

220

p = 2,

221

s = 5,

222

r = 2;

223

const PetscReal

224

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

225

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

226

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

227

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

228

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

229

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

230

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

231

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

232

{0.53638465733199574340,0.46361534266800425660},

233

{0.39901579167169582526,0.60098420832830417474},

234

{0.87689005530618575480,0.12310994469381424520},

235

{0.99056100455550913009,0.0094389954444908699092}},

236

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

237

S[2] = {1,1},

238

F[2] = {1,0},

239

Fembed[2] = {0,1},

240

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

241

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

242

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),242,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

243

}

244

{

245

#undef GAMMA0.25

246

#define GAMMA0.25 0.0

247

/* y-y~ form */

248

const PetscInt

249

p = 3,

250

s = 5,

251

r = 2;

252

const PetscReal

253

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

254

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

255

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

256

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

257

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

258

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

259

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

260

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

261

{161694774978034105510.0 / 106187653640211060371.0,-55507121337823045139.0 / 106187653640211060371.0},

262

{78486094644566264568.0 / 88171030896733822981.0,9684936252167558413.0 / 88171030896733822981.0},

263

{65394922146334854435.0 / 84570853840405479554.0,19175931694070625119.0 / 84570853840405479554.0},

264

{8607282770183754108.0 / 108658046436496925911.0,100050763666313171803.0 / 108658046436496925911.0}},

265

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

266

S[2] = {1,1},

267

F[2] = {1,0},

268

Fembed[2] = {0,1},

269

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

270

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

271

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),271,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

272

}

273

{

274

#undef GAMMA0.25

275

#define GAMMA0.25 0.25

276

/* y-eps form */

277

const PetscInt

278

p = 2,

279

s = 6,

280

r = 2;

281

const PetscReal

282

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

283

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

284

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

285

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

286

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

287

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

288

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

289

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

290

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

291

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

292

S[2] = {1,0},

293

F[2] = {1,0},

294

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

295

Ferror[2] = {0,1},

296

Serror[2] = {1,0};

297

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),297,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

298

}

299

{

300

#undef GAMMA0.25

301

#define GAMMA0.25 0.0

302

/* y-eps form */

303

const PetscInt

304

p = 3,

305

s = 8,

306

r = 2;

307

const PetscReal

308

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

309

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

310

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

311

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

312

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

313

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

314

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

315

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

316

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

317

{-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}},

318

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

319

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

320

S[2] = {1,0},

321

F[2] = {1,0},

322

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

323

Ferror[2] = {0,1},

324

Serror[2] = {1,0};

325

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),325,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

326

}

327

{

328

#undef GAMMA0.25

329

#define GAMMA0.25 0.25

330

/* y-eps form */

331

const PetscInt

332

p = 2,

333

s = 9,

334

r = 2;

335

const PetscReal

336

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

337

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

338

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

339

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

340

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

341

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

342

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

343

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

344

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

345

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

346

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

347

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}},

348

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

349

S[2] = {1,0},

350

F[2] = {1,0},

351

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

352

Ferror[2] = {0,1},

353

Serror[2] = {1,0};

354

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),354,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

355

}

356

357

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);

358

}

359

360

/*@C

361

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

362

363

Not Collective

364

365

Level: advanced

366

367

.keywords: TSGLEE, register, destroy

368

.seealso: TSGLEERegister(), TSGLEERegisterAll()

369

@*/

370

PetscErrorCode TSGLEERegisterDestroy(void)

371

{

372

PetscErrorCode ierr;

373

GLEETableauLink link;

374

375

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] = 375; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

376

while ((link = GLEETableauList)) {

377

GLEETableau t = &link->tab;

378

GLEETableauList = link->next;

379

ierr = PetscFree5(t->A,t->B,t->U,t->V,t->c)PetscFreeA(5,379,__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),379,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

380

ierr = PetscFree2(t->S,t->F)PetscFreeA(2,380,__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),380,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

381

ierr = PetscFree (t->Fembed)((*PetscTrFree)((void*)(t->Fembed),381,__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),381,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

382

ierr = PetscFree (t->Ferror)((*PetscTrFree)((void*)(t->Ferror),382,__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),382,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

383

ierr = PetscFree (t->Serror)((*PetscTrFree)((void*)(t->Serror),383,__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),383,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

384

ierr = PetscFree (t->binterp)((*PetscTrFree)((void*)(t->binterp),384,__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),384,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

385

ierr = PetscFree (t->name)((*PetscTrFree)((void*)(t->name),385,__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),385,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

386

ierr = PetscFree (link)((*PetscTrFree)((void*)(link),386,__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),386,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

387

}

388

TSGLEERegisterAllCalled = PETSC_FALSE;

389

390

}

391

392

/*@C

393

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

394

from TSInitializePackage().

395

396

Level: developer

397

398

.keywords: TS, TSGLEE, initialize, package

399

.seealso: PetscInitialize()

400

@*/

401

PetscErrorCode TSGLEEInitializePackage(void)

402

{

403

PetscErrorCode ierr;

404

405

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] = 405; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

406

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);

407

TSGLEEPackageInitialized = PETSC_TRUE;

408

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

409

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

410

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

411

412

}

413

414

/*@C

415

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

416

called from PetscFinalize().

417

418

Level: developer

419

420

.keywords: Petsc, destroy, package

421

.seealso: PetscFinalize()

422

@*/

423

PetscErrorCode TSGLEEFinalizePackage(void)

424

{

425

PetscErrorCode ierr;

426

427

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] = 427; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

428

TSGLEEPackageInitialized = PETSC_FALSE;

429

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

430

431

}

432

433

/*@C

434

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

435

436

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

437

438

Input Parameters:

439

+ name - identifier for method

440

. order - order of method

441

. s - number of stages

442

. r - number of steps

443

. gamma - LTE ratio

444

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

445

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

446

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

447

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

448

. S - starting coefficients

449

. F - finishing coefficients

450

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

451

. Fembed - step completion coefficients for embedded method

452

. Ferror - error computation coefficients

453

. Serror - error initialization coefficients

454

. pinterp - order of interpolation (0 if unavailable)

455

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

456

457

Notes:

458

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

459

460

Level: advanced

461

462

.keywords: TS, register

463

464

.seealso: TSGLEE

465

@*/

466

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

467

PetscReal gamma,

468

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

469

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

470

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

471

const PetscReal c[],

472

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

473

const PetscReal Serror[],

474

PetscInt pinterp, const PetscReal binterp[])

475

{

476

PetscErrorCode ierr;

477

GLEETableauLink link;

478

GLEETableau t;

479

PetscInt i,j;

480

481

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] = 481; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

482

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

483

ierr = PetscMalloc(sizeof(*link),&link)((*PetscTrMalloc)((sizeof(*link)),483,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,(void**)(&link)));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 = PetscMemzero(link,sizeof(*link));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

t = &link->tab;

486

ierr = PetscStrallocpy(name,&t->name);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

t->order = order;

488

t->s = s;

489

t->r = r;

490

t->gamma = gamma;

491

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,491,__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),491,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

492

ierr = PetscMalloc2(r,&t->S,r,&t->F)PetscMallocA(2,PETSC_FALSE,492,__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),492,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

493

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

494

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

495

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

496

ierr = PetscMemcpy(t->V,V,r*r*sizeof(V[0]));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 = PetscMemcpy(t->S,S,r *sizeof(S[0]));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 = PetscMemcpy(t->F,F,r *sizeof(F[0]));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

if (c) {

500

ierr = PetscMemcpy(t->c,c,s*sizeof(c[0])); 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

} else {

502

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

503

}

504

ierr = PetscMalloc1(r,&t->Fembed)PetscMallocA(1,PETSC_FALSE,504,__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),504,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

505

ierr = PetscMalloc1(r,&t->Ferror)PetscMallocA(1,PETSC_FALSE,505,__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),505,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

506

ierr = PetscMalloc1(r,&t->Serror)PetscMallocA(1,PETSC_FALSE,506,__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),506,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

507

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

508

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

509

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

510

t->pinterp = pinterp;

511

ierr = PetscMalloc1(s*pinterp,&t->binterp)PetscMallocA(1,PETSC_FALSE,511,__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),511,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

512

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

513

514

link->next = GLEETableauList;

515

GLEETableauList = link;

516

517

}

518

519

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

520

{

521

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

522

GLEETableau tab = glee->tableau;

523

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

524

*F = tab->F,

525

*Fembed = tab->Fembed;

526

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

527

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

528

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

529

PetscErrorCode ierr;

530

531

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] = 531; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

532

533

switch (glee->status) {

534

case TS_STEP_INCOMPLETE:

535

case TS_STEP_PENDING:

536

h = ts->time_step; break;

537

case TS_STEP_COMPLETE:

538

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

539

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

540

}

541

542

if (order == tab->order) {

543

544

/* Note: Irrespective of whether status is TS_STEP_INCOMPLETE

545

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

546

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

547

548

if (glee->status == TS_STEP_INCOMPLETE) {

549

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

550

ierr = VecZeroEntries(Y[i]); 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

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

552

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

553

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

554

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

555

}

556

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

557

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

558

ierr = VecMAXPY(X,r,wr,Y);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

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

560

561

562

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

563

564

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

565

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

566

ierr = VecZeroEntries(Y[i]); 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

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

568

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

569

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

570

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

571

}

572

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

573

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

574

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

575

576

if (done) *done = PETSC_TRUE;

577

578

}

579

if (done) *done = PETSC_FALSE;

580

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),580,__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);

581

582

}

583

584

static PetscErrorCode TSStep_GLEE(TS ts)

585

{

586

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

587

GLEETableau tab = glee->tableau;

588

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

589

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

590

*F = tab->F,

591

*c = tab->c;

592

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

593

*YStage = glee->YStage,

594

*YdotStage = glee->YdotStage,

595

W = glee->W;

596

SNES snes;

597

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

598

TSAdapt adapt;

599

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

600

PetscReal next_time_step;

601

PetscReal t;

602

PetscBool accept;

603

PetscErrorCode ierr;

604

605

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] = 605; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

606

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

607

608

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),608,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0); }

609

610

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

611

next_time_step = ts->time_step;

612

t = ts->ptime;

613

accept = PETSC_TRUE;

614

glee->status = TS_STEP_INCOMPLETE;

615

616

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

617

618

PetscReal h = ts->time_step;

619

ierr = TSPreStep(ts);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

621

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

622

623

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

624

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

625

626

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

627

ierr = VecZeroEntries(YStage[i]);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(YStage[i],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(YStage[i],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

} else { /* Implicit stage */

633

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

634

/* compute right-hand-side */

635

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

636

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

637

ierr = VecMAXPY(W,r,wr,X);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

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

639

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

640

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

641

/* set initial guess */

642

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),642,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

643

/* solve for this stage */

644

ierr = SNESSolve(snes,W,YStage[i]);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 = SNESGetIterationNumber(snes,&its);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

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

647

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

648

}

649

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

650

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

651

if (!accept) goto reject_step;

652

ierr = TSPostStage(ts,glee->stage_time,i,YStage); 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 = TSComputeRHSFunction(ts,t+h*c[i],YStage[i],YdotStage[i]);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

}

655

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

656

glee->status = TS_STEP_PENDING;

657

658

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

659

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

660

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

661

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),661,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

662

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),662,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

663

if (accept) {

664

/* ignore next_scheme for now */

665

ts->ptime += ts->time_step;

666

ts->time_step = next_time_step;

667

glee->status = TS_STEP_COMPLETE;

668

/* compute and store the global error */

669

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

670

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

671

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),671,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

672

break;

673

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

674

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

675

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

676

ts->time_step = next_time_step;

677

glee->status = TS_STEP_INCOMPLETE;

678

}

679

reject_step: continue;

680

}

681

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

682

683

}

684

685

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

686

{

687

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

688

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

689

PetscReal h,tt,t;

690

PetscScalar *b;

691

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

692

PetscErrorCode ierr;

693

694

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] = 694; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

695

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),695,__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);

696

switch (glee->status) {

697

case TS_STEP_INCOMPLETE:

698

case TS_STEP_PENDING:

699

h = ts->time_step;

700

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

701

break;

702

case TS_STEP_COMPLETE:

703

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

704

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

705

break;

706

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

707

}

708

ierr = PetscMalloc1(s,&b)PetscMallocA(1,PETSC_FALSE,708,__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),708,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

709

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

710

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

711

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

712

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

713

}

714

}

715

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

716

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

717

ierr = PetscFree(b)((*PetscTrFree)((void*)(b),717,__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),717,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

718

719

}

720

721

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

722

static PetscErrorCode TSReset_GLEE(TS ts)

723

{

724

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

725

PetscInt s, r;

726

PetscErrorCode ierr;

727

728

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] = 728; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

729

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);

730

s = glee->tableau->s;

731

r = glee->tableau->r;

732

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

733

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

734

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

735

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

736

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

737

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

738

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

739

ierr = PetscFree2(glee->swork,glee->rwork)PetscFreeA(2,739,__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),739,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

740

741

}

742

743

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

744

{

745

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

746

PetscErrorCode ierr;

747

748

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] = 748; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

749

if (Ydot) {

750

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

751

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

752

} else *Ydot = glee->Ydot;

753

}

754

755

}

756

757

758

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

759

{

760

PetscErrorCode ierr;

761

762

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] = 762; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

763

if (Ydot) {

764

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

765

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

766

}

767

}

768

769

}

770

771

772

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

773

774

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

775

{

776

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

777

DM dm,dmsave;

778

Vec Ydot;

779

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

780

PetscErrorCode ierr;

781

782

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] = 782; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

783

ierr = SNESGetDM(snes,&dm);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

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

785

/* Set Ydot = shift*X */

786

ierr = VecCopy(X,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

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

788

dmsave = ts->dm;

789

ts->dm = dm;

790

791

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

792

793

ts->dm = dmsave;

794

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

795

796

}

797

798

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

799

{

800

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

801

DM dm,dmsave;

802

Vec Ydot;

803

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

804

PetscErrorCode ierr;

805

806

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] = 806; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

807

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

808

ierr = TSGLEEGetVecs(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

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

810

dmsave = ts->dm;

811

ts->dm = dm;

812

813

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),813,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

814

815

ts->dm = dmsave;

816

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

817

818

}

819

820

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

821

{

822

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] = 822; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

823

824

}

825

826

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

827

{

828

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] = 828; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

829

830

}

831

832

833

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

834

{

835

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] = 835; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

836

837

}

838

839

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

840

{

841

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] = 841; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

842

843

}

844

845

static PetscErrorCode TSSetUp_GLEE(TS ts)

846

{

847

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

848

GLEETableau tab;

849

PetscInt s,r;

850

PetscErrorCode ierr;

851

DM dm;

852

853

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] = 853; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

854

if (!glee->tableau) {

Assuming pointer value is null

→

←

Taking true branch

→

855

ierr = TSGLEESetType(ts,TSGLEEDefaultType);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

}

857

tab = glee->tableau;

←

Null pointer value stored to 'tab'

→

858

s = tab->s;

←

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

859

r = tab->r;

860

ierr = VecDuplicateVecs(ts->vec_sol,r,&glee->Y);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 = VecDuplicateVecs(ts->vec_sol,r,&glee->X);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 = VecDuplicateVecs(ts->vec_sol,s,&glee->YStage);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 = VecDuplicateVecs(ts->vec_sol,s,&glee->YdotStage);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

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

865

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

866

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

867

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

868

ierr = PetscMalloc2(s,&glee->swork,r,&glee->rwork)PetscMallocA(2,PETSC_FALSE,868,__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),868,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

869

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

870

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

871

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

872

873

}

874

875

PetscErrorCode TSStartingMethod_GLEE(TS ts)

876

{

877

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

878

GLEETableau tab = glee->tableau;

879

PetscInt r=tab->r,i;

880

PetscReal *S=tab->S;

881

PetscErrorCode ierr;

882

883

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] = 883; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

884

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

885

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

886

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

887

}

888

889

890

}

891

892

893

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

894

895

static PetscErrorCode TSSetFromOptions_GLEE(PetscOptionItems *PetscOptionsObject,TS ts)

896

{

897

PetscErrorCode ierr;

898

char gleetype[256];

899

900

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] = 900; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

901

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

902

{

903

GLEETableauLink link;

904

PetscInt count,choice;

905

PetscBool flg;

906

const char **namelist;

907

908

ierr = PetscStrncpy(gleetype,TSGLEEDefaultType,sizeof(gleetype));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

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

910

ierr = PetscMalloc1(count,(char***)&namelist)PetscMallocA(1,PETSC_FALSE,910,__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),910,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

911

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

912

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),912,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

913

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

914

ierr = PetscFree(namelist)((*PetscTrFree)((void*)(namelist),914,__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),914,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

915

}

916

ierr = PetscOptionsTail()0; {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);};CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),916,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

917

918

}

919

920

static PetscErrorCode TSView_GLEE(TS ts,PetscViewer viewer)

921

{

922

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

923

GLEETableau tab = glee->tableau;

924

PetscBool iascii;

925

PetscErrorCode ierr;

926

927

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] = 927; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

928

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

929

if (iascii) {

930

TSGLEEType gleetype;

931

char buf[512];

932

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

933

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

934

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

935

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

936

/* Note: print out r as well */

937

}

938

939

}

940

941

static PetscErrorCode TSLoad_GLEE(TS ts,PetscViewer viewer)

942

{

943

PetscErrorCode ierr;

944

SNES snes;

945

TSAdapt tsadapt;

946

947

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] = 947; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

948

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

949

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

950

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

951

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

952

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

953

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

954

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),954,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

955

956

}

957

958

/*@C

959

TSGLEESetType - Set the type of GLEE scheme

960

961

Logically collective

962

963

Input Parameter:

964

+ ts - timestepping context

965

- gleetype - type of GLEE-scheme

966

967

Level: intermediate

968

969

.seealso: TSGLEEGetType(), TSGLEE

970

@*/

971

PetscErrorCode TSGLEESetType(TS ts,TSGLEEType gleetype)

972

{

973

PetscErrorCode ierr;

974

975

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] = 975; petscstack->petscroutine
[petscstack->currentsize] = PETSC_TRUE; petscstack->currentsize
++; } if (petscstack) { petscstack->hotdepth += (PETSC_FALSE
|| petscstack->hotdepth); } ; } while (0); ; } while (0);

976

PetscValidHeaderSpecific(ts,TS_CLASSID,1)do { if (!ts) return PetscError(((MPI_Comm)0x44000001),976,__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),976
,__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),976,__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),976,__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);

977

PetscValidCharPointer(gleetype,2)do { if (!gleetype) return PetscError(((MPI_Comm)0x44000001),
977,__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),977,__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);

978

ierr = PetscTryMethod(ts,"TSGLEESetType_C",(TS,TSGLEEType),(ts,gleetype))0;{ 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),978,__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),978,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0);} };CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),978,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

979

980

}

981

982

/*@C

983

TSGLEEGetType - Get the type of GLEE scheme

984

985

Logically collective

986

987

Input Parameter:

988

. ts - timestepping context

989

990

Output Parameter:

991

. gleetype - type of GLEE-scheme

992

993

Level: intermediate

994

995

.seealso: TSGLEESetType()

996

@*/

997

PetscErrorCode TSGLEEGetType(TS ts,TSGLEEType *gleetype)

998

{

999

PetscErrorCode ierr;

1000

1001

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] = 1001; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1002

PetscValidHeaderSpecific(ts,TS_CLASSID,1)do { if (!ts) return PetscError(((MPI_Comm)0x44000001),1002,__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),1002
,__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),1002,__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),1002,__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);

1003

ierr = PetscUseMethod(ts,"TSGLEEGetType_C",(TS,TSGLEEType*),(ts,gleetype))0;{ 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),1003,__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),1003,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,__ierr,PETSC_ERROR_REPEAT," ");} while (0);} else return PetscError
(PetscObjectComm((PetscObject)ts),1003,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,56,PETSC_ERROR_INITIAL,"Cannot locate function %s in object"
,"TSGLEEGetType_C"); };CHKERRQ(ierr)do {if (__builtin_expect(!!(ierr),0)) return PetscError(((MPI_Comm
)0x44000001),1003,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1004

1005

}

1006

1007

PetscErrorCode TSGLEEGetType_GLEE(TS ts,TSGLEEType *gleetype)

1008

{

1009

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

1010

PetscErrorCode ierr;

1011

1012

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] = 1012; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1013

if (!glee->tableau) {

1014

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

1015

}

1016

*gleetype = glee->tableau->name;

1017

1018

}

1019

PetscErrorCode TSGLEESetType_GLEE(TS ts,TSGLEEType gleetype)

1020

{

1021

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

1022

PetscErrorCode ierr;

1023

PetscBool match;

1024

GLEETableauLink link;

1025

1026

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] = 1026; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1027

if (glee->tableau) {

1028

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

1029

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);

1030

}

1031

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

1032

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

1033

if (match) {

1034

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

1035

glee->tableau = &link->tab;

1036

1037

}

1038

}

1039

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

1040

1041

}

1042

1043

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

1044

{

1045

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

1046

1047

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] = 1047; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1048

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

1049

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

1050

1051

}

1052

1053

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

1054

{

1055

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

1056

GLEETableau tab = glee->tableau;

1057

PetscErrorCode ierr;

1058

1059

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] = 1059; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1060

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

1061

else {

1062

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

1063

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

1064

} 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),1064,__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);

1065

}

1066

1067

}

1068

1069

PetscErrorCode TSGetAuxSolution_GLEE(TS ts,Vec *X)

1070

{

1071

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

1072

GLEETableau tab = glee->tableau;

1073

PetscReal *F = tab->Fembed;

1074

PetscInt r = tab->r;

1075

Vec *Y = glee->Y;

1076

PetscScalar *wr = glee->rwork;

1077

PetscInt i;

1078

PetscErrorCode ierr;

1079

1080

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] = 1080; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1081

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

1082

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

1083

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

1084

1085

}

1086

1087

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

1088

{

1089

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

1090

GLEETableau tab = glee->tableau;

1091

PetscReal *F = tab->Ferror;

1092

PetscInt r = tab->r;

1093

Vec *Y = glee->Y;

1094

PetscScalar *wr = glee->rwork;

1095

PetscInt i;

1096

PetscErrorCode ierr;

1097

1098

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] = 1098; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1099

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

1100

if(n==0){

1101

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

1102

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

1103

} else if(n==-1) {

1104

*X=glee->yGErr;

1105

}

1106

1107

}

1108

1109

PetscErrorCode TSSetTimeError_GLEE(TS ts,Vec X)

1110

{

1111

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

1112

GLEETableau tab = glee->tableau;

1113

PetscReal *S = tab->Serror;

1114

PetscInt r = tab->r,i;

1115

Vec *Y = glee->Y;

1116

PetscErrorCode ierr;

1117

1118

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] = 1118; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1119

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),1119,__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);

1120

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

1121

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

1122

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

1123

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

1124

}

1125

1126

}

1127

1128

static PetscErrorCode TSDestroy_GLEE(TS ts)

1129

{

1130

PetscErrorCode ierr;

1131

1132

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] = 1132; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1133

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

1134

if (ts->dm) {

1135

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

1136

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

1137

}

1138

ierr = PetscFree(ts->data)((*PetscTrFree)((void*)(ts->data),1138,__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),1138,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1139

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),1139,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1140

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),1140,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1141

1142

}

1143

1144

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

1145

/*MC

1146

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

1147

1148

The user should provide the right hand side of the equation

1149

using TSSetRHSFunction().

1150

1151

Notes:

1152

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

1153

1154

Level: beginner

1155

1156

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

1157

TSGLEE23, TTSGLEE24, TSGLEE35, TSGLEE25I, TSGLEEEXRK2A,

1158

TSGLEERK32G1, TSGLEERK285EX, TSGLEEType, TSGLEERegister()

1159

1160

M*/

1161

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

1162

{

1163

TS_GLEE *th;

1164

PetscErrorCode ierr;

1165

1166

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] = 1166; petscstack->
petscroutine[petscstack->currentsize] = PETSC_TRUE; petscstack
->currentsize++; } if (petscstack) { petscstack->hotdepth
+= (PETSC_FALSE || petscstack->hotdepth); } ; } while (0)
; ; } while (0);

1167

#if !defined(PETSC_USE_DYNAMIC_LIBRARIES)

1168

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

1169

#endif

1170

1171

ts->ops->reset = TSReset_GLEE;

1172

ts->ops->destroy = TSDestroy_GLEE;

1173

ts->ops->view = TSView_GLEE;

1174

ts->ops->load = TSLoad_GLEE;

1175

ts->ops->setup = TSSetUp_GLEE;

1176

ts->ops->step = TSStep_GLEE;

1177

ts->ops->interpolate = TSInterpolate_GLEE;

1178

ts->ops->evaluatestep = TSEvaluateStep_GLEE;

1179

ts->ops->setfromoptions = TSSetFromOptions_GLEE;

1180

ts->ops->getstages = TSGetStages_GLEE;

1181

ts->ops->snesfunction = SNESTSFormFunction_GLEE;

1182

ts->ops->snesjacobian = SNESTSFormJacobian_GLEE;

1183

ts->ops->getsolutioncomponents = TSGetSolutionComponents_GLEE;

1184

ts->ops->getauxsolution = TSGetAuxSolution_GLEE;

1185

ts->ops->gettimeerror = TSGetTimeError_GLEE;

1186

ts->ops->settimeerror = TSSetTimeError_GLEE;

1187

ts->ops->startingmethod = TSStartingMethod_GLEE;

1188

ts->default_adapt_type = TSADAPTGLEE"glee";

1189

1190

ts->usessnes = PETSC_TRUE;

1191

1192

ierr = PetscNewLog(ts,&th)(PetscMallocA(1,PETSC_TRUE,1192,__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),1192,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1193

ts->data = (void*)th;

1194

1195

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),1195,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1196

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),1196,__func__,"/sandbox/petsc/petsc.next-tmp/src/ts/impls/glee/glee.c"
,ierr,PETSC_ERROR_REPEAT," ");} while (0);

1197

1198

}