Richardson Extrapolation-Based Verification Method of Scientific Calculation Program without the Oracles: A Case Study

<table class="fixed-width table-group" id="tab1"><tr><td><table class="table"><colgroup><col style="width:6.74em"/><col style="width:5.25em"/><col style="width:5.21em"/><col style="width:5.67em"/><col style="width:4.93em"/><col style="width:5.96em"/></colgroup><tr><td class="thead-hr" colspan="6"><hr/></td></tr><tr class="thead"><td class="align_left">Energy group</td><td class="align_center"><span style="width: 31.0643ptpx;"><svg height="11.439pt" id="M99" style="vertical-align:-2.15067pt" version="1.1" viewbox="-0.0498162 -9.28833 31.0643 11.439" width="31.0643pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M743 375C743 570 607 650 407 650H136L130 622C214 612 220 606 204 526L125 122C110 43 105 36 24 27L17 0H250C382 0 482 23 571 76C680 141 743 248 743 375ZM644 379C644 188 520 36 305 36C277 36 249 38 229 46C205 55 202 71 212 126L293 553C300 588 306 600 317 606C330 613 356 617 390 617C547 617 644 537 644 379Z"></path></g><g transform="matrix(.013,0,0,-0.013,9.958,0)"><path d="M368 703H309L44 -163H104L368 703Z"></path></g><g transform="matrix(.013,0,0,-0.013,15.319,0)"><path d="M383 397C383 424 344 448 299 448C244 448 172 409 132 375C66 319 23 227 23 146C23 42 74 -12 146 -12C208 -12 298 30 359 103L343 124C315 95 248 48 192 48C145 48 111 85 111 163C111 228 129 294 151 330C171 363 201 401 241 401C275 401 302 384 325 356C332 347 339 344 348 348C373 360 383 381 383 397Z"></path></g><g transform="matrix(.013,0,0,-0.013,20.602,0)"><path d="M766 88L752 113C719 83 690 64 681 64C674 64 672 74 679 103L724 292C758 436 724 448 701 448C680 448 666 442 639 429C594 407 514 350 441 252H439L447 289C476 423 450 448 419 448C398 448 379 441 355 427C307 400 234 344 162 249H160L180 324C203 409 197 448 170 448C144 448 82 413 23 349L35 321C57 343 96 374 108 374C115 374 117 371 111 341C87 227 57 112 24 -6L32 -12C53 -4 81 4 108 6C119 68 134 128 149 171C177 229 309 383 364 383C387 383 388 355 373 282C354 190 330 92 303 -6L309 -12C332 -4 356 3 386 6C396 63 411 122 424 171C458 236 590 383 642 383C658 383 664 369 652 315L603 91C587 20 593 -12 619 -12C642 -12 708 23 766 88Z"></path></g></svg></span></td><td class="align_center"> <svg height="14.8825pt" id="M100" style="vertical-align:-3.2911pt" version="1.1" viewbox="-0.0498162 -11.5914 43.0671 14.8825" width="43.0671pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M570 168L541 175C510 88 481 71 415 71H135L378 358L167 614H348C422 614 447 608 463 588C478 571 488 546 500 500L529 504C524 559 519 633 518 650H47V628L295 324L33 23V0H531L570 168Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,7.774,3.132)"><path d="M329 433H203L239 587L230 596L147 534L123 433H57L30 395L34 388H115L61 129C37 16 59 -12 85 -12C147 -12 222 58 260 98L241 125C212 95 160 62 144 62C132 62 127 71 138 126L192 386L305 394L329 433Z"></path></g><g transform="matrix(.013,0,0,-0.013,11.442,0)"><path d="M368 703H309L44 -163H104L368 703Z"></path></g><g transform="matrix(.013,0,0,-0.013,16.803,0)"><path d="M383 397C383 424 344 448 299 448C244 448 172 409 132 375C66 319 23 227 23 146C23 42 74 -12 146 -12C208 -12 298 30 359 103L343 124C315 95 248 48 192 48C145 48 111 85 111 163C111 228 129 294 151 330C171 363 201 401 241 401C275 401 302 384 325 356C332 347 339 344 348 348C373 360 383 381 383 397Z"></path></g><g transform="matrix(.013,0,0,-0.013,22.085,0)"><path d="M766 88L752 113C719 83 690 64 681 64C674 64 672 74 679 103L724 292C758 436 724 448 701 448C680 448 666 442 639 429C594 407 514 350 441 252H439L447 289C476 423 450 448 419 448C398 448 379 441 355 427C307 400 234 344 162 249H160L180 324C203 409 197 448 170 448C144 448 82 413 23 349L35 321C57 343 96 374 108 374C115 374 117 371 111 341C87 227 57 112 24 -6L32 -12C53 -4 81 4 108 6C119 68 134 128 149 171C177 229 309 383 364 383C387 383 388 355 373 282C354 190 330 92 303 -6L309 -12C332 -4 356 3 386 6C396 63 411 122 424 171C458 236 590 383 642 383C658 383 664 369 652 315L603 91C587 20 593 -12 619 -12C642 -12 708 23 766 88Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,32.342,-5.741)"><path d="M556 236V289H56V236H556Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,37.902,-5.741)"><path d="M389 0V32C297 38 291 46 291 118V635C234 613 175 595 109 583V556L161 554C203 552 207 547 207 497V118C207 46 201 38 110 32V0H389Z"></path></g></svg></td><td class="align_center"><span style="width: 51.6649ptpx;"><svg height="17.1204pt" id="M101" style="vertical-align:-5.528999pt" version="1.1" viewbox="-0.0498162 -11.5914 51.6649 17.1204" width="51.6649pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M219 86C216 168 211 250 206 337C201 410 189 448 163 448C131 448 79 396 43 344L60 322C91 359 110 375 118 375S132 358 136 298C141 238 152 81 155 -12H182C242 62 331 177 390 258C435 321 451 360 451 391C450 424 432 448 408 448C390 448 372 435 366 419C362 410 362 401 366 394C373 383 376 367 376 350C376 283 262 138 221 86H219Z"></path></g><g transform="matrix(.013,0,0,-0.013,5.863,0)"><path d="M570 168L541 175C510 88 481 71 415 71H135L378 358L167 614H348C422 614 447 608 463 588C478 571 488 546 500 500L529 504C524 559 519 633 518 650H47V628L295 324L33 23V0H531L570 168Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,13.637,3.132)"><path d="M620 667C620 684 593 710 552 710C515 710 459 687 409 634C360 582 334 506 320 433H247L215 400L220 388H310L259 80C224 -129 199 -161 185 -174C175 -183 160 -192 141 -192C123 -192 90 -181 75 -165C69 -158 64 -156 54 -164C39 -178 24 -195 24 -208C24 -231 61 -257 94 -257C124 -257 162 -243 209 -194C267 -134 303 -66 343 107C369 220 378 279 393 385L503 397L523 433H404C432 619 464 660 500 660C522 660 542 646 566 621C573 614 583 615 593 621C603 628 620 648 620 667Z"></path></g><g transform="matrix(.013,0,0,-0.013,20.018,0)"><path d="M368 703H309L44 -163H104L368 703Z"></path></g><g transform="matrix(.013,0,0,-0.013,25.379,0)"><path d="M383 397C383 424 344 448 299 448C244 448 172 409 132 375C66 319 23 227 23 146C23 42 74 -12 146 -12C208 -12 298 30 359 103L343 124C315 95 248 48 192 48C145 48 111 85 111 163C111 228 129 294 151 330C171 363 201 401 241 401C275 401 302 384 325 356C332 347 339 344 348 348C373 360 383 381 383 397Z"></path></g><g transform="matrix(.013,0,0,-0.013,30.661,0)"><path d="M766 88L752 113C719 83 690 64 681 64C674 64 672 74 679 103L724 292C758 436 724 448 701 448C680 448 666 442 639 429C594 407 514 350 441 252H439L447 289C476 423 450 448 419 448C398 448 379 441 355 427C307 400 234 344 162 249H160L180 324C203 409 197 448 170 448C144 448 82 413 23 349L35 321C57 343 96 374 108 374C115 374 117 371 111 341C87 227 57 112 24 -6L32 -12C53 -4 81 4 108 6C119 68 134 128 149 171C177 229 309 383 364 383C387 383 388 355 373 282C354 190 330 92 303 -6L309 -12C332 -4 356 3 386 6C396 63 411 122 424 171C458 236 590 383 642 383C658 383 664 369 652 315L603 91C587 20 593 -12 619 -12C642 -12 708 23 766 88Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,40.918,-5.741)"><path d="M556 236V289H56V236H556Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,46.478,-5.741)"><path d="M389 0V32C297 38 291 46 291 118V635C234 613 175 595 109 583V556L161 554C203 552 207 547 207 497V118C207 46 201 38 110 32V0H389Z"></path></g></svg></span></td><td class="align_center"><span style="width: 31.091ptpx;"><svg height="14.7625pt" id="M102" style="vertical-align:-5.47417pt" version="1.1" viewbox="-0.0498162 -9.28833 31.091 14.7625" width="31.091pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M565 434L556 448L482 423L322 189L287 373C280 410 261 448 236 448C207 448 161 415 110 338L131 317C160 355 181 374 193 374C204 374 208 366 217 330C230 279 247 185 259 116L50 -175C18 -223 29 -260 41 -278H51C97 -198 207 -33 270 57L313 -165C328 -235 349 -261 376 -261S445 -244 510 -150L487 -127C444 -178 423 -186 411 -186C402 -186 398 -177 387 -130C371 -62 349 39 332 135C390 217 524 392 565 434Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,6.981,3.132)"><path d="M573 302C573 402 527 451 414 451C386 451 343 446 313 437L330 513L320 522C295 508 261 484 243 463L230 415C194 400 159 383 126 359L131 330C159 344 187 357 222 368L109 -147C96 -204 80 -214 18 -223L13 -255L256 -244L259 -212L236 -210C184 -205 180 -195 191 -141L219 -1C240 -10 268 -12 284 -12C352 4 431 48 484 104C543 166 573 240 573 302ZM481 290C481 165 381 37 305 37C280 37 249 56 235 71L302 395C328 399 353 402 372 402C427 402 481 378 481 290Z"></path></g><g transform="matrix(.013,0,0,-0.013,12.932,0)"><path d="M368 703H309L44 -163H104L368 703Z"></path></g><g transform="matrix(.013,0,0,-0.013,18.293,0)"><path d="M565 434L556 448L482 423L322 189L287 373C280 410 261 448 236 448C207 448 161 415 110 338L131 317C160 355 181 374 193 374C204 374 208 366 217 330C230 279 247 185 259 116L50 -175C18 -223 29 -260 41 -278H51C97 -198 207 -33 270 57L313 -165C328 -235 349 -261 376 -261S445 -244 510 -150L487 -127C444 -178 423 -186 411 -186C402 -186 398 -177 387 -130C371 -62 349 39 332 135C390 217 524 392 565 434Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,25.274,3.132)"><path d="M538 684C542 703 539 710 529 710C512 710 455 682 363 674L361 644H396C444 644 448 642 438 594L405 438C373 450 347 451 336 451C288 451 200 416 147 374C69 312 24 207 24 113C24 29 60 -12 94 -12C126 -12 156 4 196 30C231 53 296 100 346 166H348L328 81C311 10 325 -12 348 -12C380 -12 450 24 512 98L495 125C466 97 431 70 417 70C408 70 408 80 411 98L538 684ZM390 373L361 240C336 196 210 56 145 56C132 56 113 73 113 131C113 216 156 338 220 379C243 394 266 402 303 402C331 402 375 388 390 373Z"></path></g></svg></span></td><td class="align_center"><span style="width: 62.0478ptpx;"><svg height="14.7729pt" id="M103" style="vertical-align:-3.181499pt" version="1.1" viewbox="-0.0498162 -11.5914 62.0478 14.7729" width="62.0478pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M570 168L541 175C510 88 481 71 415 71H135L378 358L167 614H348C422 614 447 608 463 588C478 571 488 546 500 500L529 504C524 559 519 633 518 650H47V628L295 324L33 23V0H531L570 168Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,7.774,3.132)"><path d="M389 0V32C297 38 291 46 291 118V635C234 613 175 595 109 583V556L161 554C203 552 207 547 207 497V118C207 46 201 38 110 32V0H389Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,12.206,3.132)"><path d="M588 236V289H72V236H588Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,16.394,3.132)"><path d="M935 263C895 299 821 372 745 456L719 439L824 289H72V236H823L719 87L745 70C821 154 895 227 935 263Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,25.459,3.132)"><path d="M414 144C384 79 371 75 317 75H135L276 221C367 316 408 376 408 465C408 570 327 635 237 635C179 635 131 609 100 575L42 494L67 471C94 510 138 565 205 565C277 565 321 517 321 435C321 348 258 270 195 195C146 137 88 81 33 26V0H411C423 44 433 88 446 135L414 144Z"></path></g><g transform="matrix(.013,0,0,-0.013,30.422,0)"><path d="M368 703H309L44 -163H104L368 703Z"></path></g><g transform="matrix(.013,0,0,-0.013,35.783,0)"><path d="M383 397C383 424 344 448 299 448C244 448 172 409 132 375C66 319 23 227 23 146C23 42 74 -12 146 -12C208 -12 298 30 359 103L343 124C315 95 248 48 192 48C145 48 111 85 111 163C111 228 129 294 151 330C171 363 201 401 241 401C275 401 302 384 325 356C332 347 339 344 348 348C373 360 383 381 383 397Z"></path></g><g transform="matrix(.013,0,0,-0.013,41.066,0)"><path d="M766 88L752 113C719 83 690 64 681 64C674 64 672 74 679 103L724 292C758 436 724 448 701 448C680 448 666 442 639 429C594 407 514 350 441 252H439L447 289C476 423 450 448 419 448C398 448 379 441 355 427C307 400 234 344 162 249H160L180 324C203 409 197 448 170 448C144 448 82 413 23 349L35 321C57 343 96 374 108 374C115 374 117 371 111 341C87 227 57 112 24 -6L32 -12C53 -4 81 4 108 6C119 68 134 128 149 171C177 229 309 383 364 383C387 383 388 355 373 282C354 190 330 92 303 -6L309 -12C332 -4 356 3 386 6C396 63 411 122 424 171C458 236 590 383 642 383C658 383 664 369 652 315L603 91C587 20 593 -12 619 -12C642 -12 708 23 766 88Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,51.323,-5.741)"><path d="M556 236V289H56V236H556Z"></path></g><g transform="matrix(.0091,0,0,-0.0091,56.883,-5.741)"><path d="M389 0V32C297 38 291 46 291 118V635C234 613 175 595 109 583V556L161 554C203 552 207 547 207 497V118C207 46 201 38 110 32V0H389Z"></path></g></svg></span></td></tr><tr><td class="thead-hr" colspan="6"><hr/></td></tr><tr><td class="align_left">1</td><td class="align_center">1.264E-00</td><td class="align_center">8.154E-03</td><td class="align_center">0.000E-00</td><td class="align_center">1.00E-00</td><td class="align_center">7.368E-03</td></tr><tr><td class="align_left">2</td><td class="align_center">9.328E-01</td><td class="align_center">4.100E-03</td><td class="align_center">4.562E-03</td><td class="align_center">0.00E-00</td><td class="align_center">7.368E-03</td></tr><tr class="table-tr"><td colspan="6"><hr class="tbody-hr"/></td></tr></table></td></tr></table>

<div>Parameters of SLAB_1D_2G one-dimensional slab benchmark problem.</div>

Mathematical Problems in Engineering

tab1

Table 1

Table 1: Richardson Extrapolation-Based Verification Method of Scientific Calculation Program without the Oracles: A Case Study 