doc/spline3.lyx

   1 #LyX 1.6.8 created this file. For more info see http://www.lyx.org/
   2 \lyxformat 345
   3 \begin_document
   4 \begin_header
   5 \textclass amsart
   6 \use_default_options false
   7 \begin_modules
   8 theorems-ams
   9 eqs-within-sections
  10 figs-within-sections
  11 \end_modules
  12 \language english
  13 \inputencoding auto
  14 \font_roman ae
  15 \font_sans default
  16 \font_typewriter default
  17 \font_default_family default
  18 \font_sc false
  19 \font_osf false
  20 \font_sf_scale 100
  21 \font_tt_scale 100
  22
  23 \graphics default
  24 \paperfontsize default
  25 \spacing single
  26 \use_hyperref false
  27 \papersize letterpaper
  28 \use_geometry false
  29 \use_amsmath 1
  30 \use_esint 1
  31 \cite_engine basic
  32 \use_bibtopic false
  33 \paperorientation portrait
  34 \secnumdepth 3
  35 \tocdepth 3
  36 \paragraph_separation skip
  37 \defskip bigskip
  38 \quotes_language english
  39 \papercolumns 1
  40 \papersides 1
  41 \paperpagestyle default
  42 \tracking_changes false
  43 \output_changes false
  44 \author ""
  45 \author ""
  46 \end_header
  47
  48 \begin_body
  49
  50 \begin_layout Title
  51 Spline3
  52 \end_layout
  53
  54 \begin_layout Author
  55 Michael Orlitzky
  56 \end_layout
  57
  58 \begin_layout Example*
  59 First, will will define our grid
  60 \begin_inset Formula $G$
  61 \end_inset
  62
  63  to be over a
  64 \begin_inset Formula $3\times3\times3$
  65 \end_inset
  66
  67  cube having grid size
  68 \begin_inset Formula $h=1$
  69 \end_inset
  70
  71 .
  72 \end_layout
  73
  74 \begin_layout Example*
  75 To test the reproduction of trilinears, we will take the function (per Sorokina
  76  and Zeilfelder, p.
  77  88),
  78 \end_layout
  79
  80 \begin_layout Example*
  81 \begin_inset Formula \[
  82 f\left(x,y,z\right)=1+x+xy+xyz\]
  83
  84 \end_inset
  85
  86
  87 \end_layout
  88
  89 \begin_layout Example*
  90 and compute its value at the 27 points,
  91 \begin_inset Formula $\left(0,0,0\right),\left(0,0,1\right),\dots\left(2,2,2\right)$
  92 \end_inset
  93
  94 .
  95 \end_layout
  96
  97 \begin_layout Example*
  98 \begin_inset listings
  99 inline false
 100 status open
 101
 102 \begin_layout Plain Layout
 103
 104 sage: x,y,z = var('x,y,z')
 105 \end_layout
 106
 107 \begin_layout Plain Layout
 108
 109 sage: f = 1 + x + x*y + x*y*z
 110 \end_layout
 111
 112 \begin_layout Plain Layout
 113
 114 sage: f(x=0, y=0, z=0)
 115 \end_layout
 116
 117 \begin_layout Plain Layout
 118
 119 1
 120 \end_layout
 121
 122 \begin_layout Plain Layout
 123
 124 sage: f(x=1, y=0, z=0)
 125 \end_layout
 126
 127 \begin_layout Plain Layout
 128
 129 2
 130 \end_layout
 131
 132 \begin_layout Plain Layout
 133
 134 sage: f(x=2, y=0, z=0)
 135 \end_layout
 136
 137 \begin_layout Plain Layout
 138
 139 3
 140 \end_layout
 141
 142 \begin_layout Plain Layout
 143
 144 sage: f(x=0, y=1, z=0)
 145 \end_layout
 146
 147 \begin_layout Plain Layout
 148
 149 1
 150 \end_layout
 151
 152 \begin_layout Plain Layout
 153
 154 sage: f(x=1, y=1, z=0)
 155 \end_layout
 156
 157 \begin_layout Plain Layout
 158
 159 3
 160 \end_layout
 161
 162 \begin_layout Plain Layout
 163
 164 sage: f(x=2, y=1, z=0)
 165 \end_layout
 166
 167 \begin_layout Plain Layout
 168
 169 5
 170 \end_layout
 171
 172 \begin_layout Plain Layout
 173
 174 sage: f(x=0, y=2, z=0)
 175 \end_layout
 176
 177 \begin_layout Plain Layout
 178
 179 1
 180 \end_layout
 181
 182 \begin_layout Plain Layout
 183
 184 sage: f(x=1, y=2, z=0)
 185 \end_layout
 186
 187 \begin_layout Plain Layout
 188
 189 4
 190 \end_layout
 191
 192 \begin_layout Plain Layout
 193
 194 sage: f(x=2, y=2, z=0)
 195 \end_layout
 196
 197 \begin_layout Plain Layout
 198
 199 7
 200 \end_layout
 201
 202 \begin_layout Plain Layout
 203
 204 sage: f(x=0, y=0, z=1)
 205 \end_layout
 206
 207 \begin_layout Plain Layout
 208
 209 1
 210 \end_layout
 211
 212 \begin_layout Plain Layout
 213
 214 sage: f(x=1, y=0, z=1)
 215 \end_layout
 216
 217 \begin_layout Plain Layout
 218
 219 2
 220 \end_layout
 221
 222 \begin_layout Plain Layout
 223
 224 sage: f(x=2, y=0, z=1)
 225 \end_layout
 226
 227 \begin_layout Plain Layout
 228
 229 3
 230 \end_layout
 231
 232 \begin_layout Plain Layout
 233
 234 sage: f(x=0, y=1, z=1)
 235 \end_layout
 236
 237 \begin_layout Plain Layout
 238
 239 1
 240 \end_layout
 241
 242 \begin_layout Plain Layout
 243
 244 sage: f(x=1, y=1, z=1)
 245 \end_layout
 246
 247 \begin_layout Plain Layout
 248
 249 4
 250 \end_layout
 251
 252 \begin_layout Plain Layout
 253
 254 sage: f(x=2, y=1, z=1)
 255 \end_layout
 256
 257 \begin_layout Plain Layout
 258
 259 7
 260 \end_layout
 261
 262 \begin_layout Plain Layout
 263
 264 sage: f(x=0, y=2, z=1)
 265 \end_layout
 266
 267 \begin_layout Plain Layout
 268
 269 1
 270 \end_layout
 271
 272 \begin_layout Plain Layout
 273
 274 sage: f(x=1, y=2, z=1)
 275 \end_layout
 276
 277 \begin_layout Plain Layout
 278
 279 6
 280 \end_layout
 281
 282 \begin_layout Plain Layout
 283
 284 sage: f(x=2, y=2, z=1)
 285 \end_layout
 286
 287 \begin_layout Plain Layout
 288
 289 11
 290 \end_layout
 291
 292 \begin_layout Plain Layout
 293
 294 sage: f(x=0, y=0, z=2)
 295 \end_layout
 296
 297 \begin_layout Plain Layout
 298
 299 1
 300 \end_layout
 301
 302 \begin_layout Plain Layout
 303
 304 sage: f(x=1, y=0, z=2)
 305 \end_layout
 306
 307 \begin_layout Plain Layout
 308
 309 2
 310 \end_layout
 311
 312 \begin_layout Plain Layout
 313
 314 sage: f(x=2, y=0, z=2)
 315 \end_layout
 316
 317 \begin_layout Plain Layout
 318
 319 3
 320 \end_layout
 321
 322 \begin_layout Plain Layout
 323
 324 sage: f(x=0, y=1, z=2)
 325 \end_layout
 326
 327 \begin_layout Plain Layout
 328
 329 1
 330 \end_layout
 331
 332 \begin_layout Plain Layout
 333
 334 sage: f(x=1, y=1, z=2)
 335 \end_layout
 336
 337 \begin_layout Plain Layout
 338
 339 5
 340 \end_layout
 341
 342 \begin_layout Plain Layout
 343
 344 sage: f(x=2, y=1, z=2)
 345 \end_layout
 346
 347 \begin_layout Plain Layout
 348
 349 9
 350 \end_layout
 351
 352 \begin_layout Plain Layout
 353
 354 sage: f(x=0, y=2, z=2)
 355 \end_layout
 356
 357 \begin_layout Plain Layout
 358
 359 1
 360 \end_layout
 361
 362 \begin_layout Plain Layout
 363
 364 sage: f(x=1, y=2, z=2)
 365 \end_layout
 366
 367 \begin_layout Plain Layout
 368
 369 8
 370 \end_layout
 371
 372 \begin_layout Plain Layout
 373
 374 sage: f(x=2, y=2, z=2)
 375 \end_layout
 376
 377 \begin_layout Plain Layout
 378
 379 15
 380 \end_layout
 381
 382 \end_inset
 383
 384
 385 \end_layout
 386
 387 \begin_layout Example*
 388 We are most interested in the
 389 \begin_inset Quotes eld
 390 \end_inset
 391
 392 interior
 393 \begin_inset Quotes erd
 394 \end_inset
 395
 396  cube centered at
 397 \begin_inset Formula $\left(1,1,1\right)$
 398 \end_inset
 399
 400  with datum
 401 \begin_inset Formula $4$
 402 \end_inset
 403
 404 .
 405  We can enter the data above into a list,
 406 \end_layout
 407
 408 \begin_layout Example*
 409 \begin_inset listings
 410 inline false
 411 status open
 412
 413 \begin_layout Plain Layout
 414
 415 sage: g = [[[ f(x=a, y=b, z=c) for a in range(0,3) ] for b in range(0,3)
 416  ] for c in range(0,3) ]
 417 \end_layout
 418
 419 \begin_layout Plain Layout
 420
 421 sage: g
 422 \end_layout
 423
 424 \begin_layout Plain Layout
 425
 426 [[[1, 2, 3], [1, 3, 5], [1, 4, 7]],
 427 \end_layout
 428
 429 \begin_layout Plain Layout
 430
 431  [[1, 2, 3], [1, 4, 7], [1, 6, 11]],
 432 \end_layout
 433
 434 \begin_layout Plain Layout
 435
 436  [[1, 2, 3], [1, 5, 9], [1, 8, 15]]]
 437 \end_layout
 438
 439 \end_inset
 440
 441
 442 \end_layout
 443
 444 \begin_layout Example*
 445 although the list will be indexed by
 446 \begin_inset Formula $\left(z,y,x\right)$
 447 \end_inset
 448
 449  so we define a function to access it by
 450 \begin_inset Formula $\left(x,y,z\right)$
 451 \end_inset
 452
 453 ,
 454 \end_layout
 455
 456 \begin_layout Example*
 457 \begin_inset listings
 458 inline false
 459 status open
 460
 461 \begin_layout Plain Layout
 462
 463 sage: def grid(x,y,z):
 464 \end_layout
 465
 466 \begin_layout Plain Layout
 467
 468 ....:     return g[z][y][x]
 469 \end_layout
 470
 471 \end_inset
 472
 473
 474 \end_layout
 475
 476 \begin_layout Example*
 477 and define directional functions according to Sorokina and Zeilfelder, p.
 478  81.
 479 \end_layout
 480
 481 \begin_layout Example*
 482 \begin_inset listings
 483 inline false
 484 status open
 485
 486 \begin_layout Plain Layout
 487
 488 sage: def I(x,y,z):
 489 \end_layout
 490
 491 \begin_layout Plain Layout
 492
 493 ....:     return grid(x,y,z)
 494 \end_layout
 495
 496 \begin_layout Plain Layout
 497
 498 \end_layout
 499
 500 \begin_layout Plain Layout
 501
 502 sage: def F(x,y,z):
 503 \end_layout
 504
 505 \begin_layout Plain Layout
 506
 507 ....:     return grid(x-1, y, z)
 508 \end_layout
 509
 510 \begin_layout Plain Layout
 511
 512 \end_layout
 513
 514 \begin_layout Plain Layout
 515
 516 sage: def B(x,y,z):
 517 \end_layout
 518
 519 \begin_layout Plain Layout
 520
 521 ....:     return grid(x+1, y, z)
 522 \end_layout
 523
 524 \begin_layout Plain Layout
 525
 526 \end_layout
 527
 528 \begin_layout Plain Layout
 529
 530 sage: def L(x,y,z):
 531 \end_layout
 532
 533 \begin_layout Plain Layout
 534
 535 ....:     return grid(x, y-1, z)
 536 \end_layout
 537
 538 \begin_layout Plain Layout
 539
 540 \end_layout
 541
 542 \begin_layout Plain Layout
 543
 544 sage: def R(x,y,z):
 545 \end_layout
 546
 547 \begin_layout Plain Layout
 548
 549 ....:     return grid(x, y+1, z)
 550 \end_layout
 551
 552 \begin_layout Plain Layout
 553
 554 \end_layout
 555
 556 \begin_layout Plain Layout
 557
 558 sage: def T(x,y,z):
 559 \end_layout
 560
 561 \begin_layout Plain Layout
 562
 563 ....:     return grid(x, y, z+1)
 564 \end_layout
 565
 566 \begin_layout Plain Layout
 567
 568 \end_layout
 569
 570 \begin_layout Plain Layout
 571
 572 sage: def D(x,y,z):
 573 \end_layout
 574
 575 \begin_layout Plain Layout
 576
 577 ....:     return grid(x, y, z-1)
 578 \end_layout
 579
 580 \begin_layout Plain Layout
 581
 582 \end_layout
 583
 584 \begin_layout Plain Layout
 585
 586 sage: def FL(x,y,z):
 587 \end_layout
 588
 589 \begin_layout Plain Layout
 590
 591 ....:     return grid(x-1, y-1, z)
 592 \end_layout
 593
 594 \begin_layout Plain Layout
 595
 596 \end_layout
 597
 598 \begin_layout Plain Layout
 599
 600 sage: def FR(x,y,z):
 601 \end_layout
 602
 603 \begin_layout Plain Layout
 604
 605 ....:     return grid(x-1, y+1, z)
 606 \end_layout
 607
 608 \begin_layout Plain Layout
 609
 610 \end_layout
 611
 612 \begin_layout Plain Layout
 613
 614 sage: def FD(x,y,z):
 615 \end_layout
 616
 617 \begin_layout Plain Layout
 618
 619 ....:     return grid(x-1, y, z-1)
 620 \end_layout
 621
 622 \begin_layout Plain Layout
 623
 624 \end_layout
 625
 626 \begin_layout Plain Layout
 627
 628 sage: def FT(x,y,z):
 629 \end_layout
 630
 631 \begin_layout Plain Layout
 632
 633 ....:     return grid(x-1, y, z+1)
 634 \end_layout
 635
 636 \begin_layout Plain Layout
 637
 638 \end_layout
 639
 640 \begin_layout Plain Layout
 641
 642 sage: def BL(x,y,z):
 643 \end_layout
 644
 645 \begin_layout Plain Layout
 646
 647 ....:     return grid(x+1, y-1, z)
 648 \end_layout
 649
 650 \begin_layout Plain Layout
 651
 652 \end_layout
 653
 654 \begin_layout Plain Layout
 655
 656 sage: def BR(x,y,z):
 657 \end_layout
 658
 659 \begin_layout Plain Layout
 660
 661 ....:     return grid(x+1, y+1, z)
 662 \end_layout
 663
 664 \begin_layout Plain Layout
 665
 666 \end_layout
 667
 668 \begin_layout Plain Layout
 669
 670 sage: def BD(x,y,z):
 671 \end_layout
 672
 673 \begin_layout Plain Layout
 674
 675 ....:     return grid(x+1, y, z-1)
 676 \end_layout
 677
 678 \begin_layout Plain Layout
 679
 680 \end_layout
 681
 682 \begin_layout Plain Layout
 683
 684 sage: def BT(x,y,z):
 685 \end_layout
 686
 687 \begin_layout Plain Layout
 688
 689 ....:     return grid(x+1, y, z+1)
 690 \end_layout
 691
 692 \begin_layout Plain Layout
 693
 694 \end_layout
 695
 696 \begin_layout Plain Layout
 697
 698 sage: def LD(x,y,z):
 699 \end_layout
 700
 701 \begin_layout Plain Layout
 702
 703 ....:     return grid(x, y-1, z-1)
 704 \end_layout
 705
 706 \begin_layout Plain Layout
 707
 708 \end_layout
 709
 710 \begin_layout Plain Layout
 711
 712 sage: def LT(x,y,z):
 713 \end_layout
 714
 715 \begin_layout Plain Layout
 716
 717 ....:     return grid(x, y-1, z+1)
 718 \end_layout
 719
 720 \begin_layout Plain Layout
 721
 722 \end_layout
 723
 724 \begin_layout Plain Layout
 725
 726 sage: def RD(x,y,z):
 727 \end_layout
 728
 729 \begin_layout Plain Layout
 730
 731 ....:     return grid(x, y+1, z-1)
 732 \end_layout
 733
 734 \begin_layout Plain Layout
 735
 736 \end_layout
 737
 738 \begin_layout Plain Layout
 739
 740 sage: def RT(x,y,z):
 741 \end_layout
 742
 743 \begin_layout Plain Layout
 744
 745 ....:     return grid(x, y+1, z+1)
 746 \end_layout
 747
 748 \begin_layout Plain Layout
 749
 750 \end_layout
 751
 752 \begin_layout Plain Layout
 753
 754 sage: def FLD(x,y,z):
 755 \end_layout
 756
 757 \begin_layout Plain Layout
 758
 759 ....:     return grid(x-1, y-1, z-1)
 760 \end_layout
 761
 762 \begin_layout Plain Layout
 763
 764 \end_layout
 765
 766 \begin_layout Plain Layout
 767
 768 sage: def FLT(x,y,z):
 769 \end_layout
 770
 771 \begin_layout Plain Layout
 772
 773 ....:     return grid(x-1, y-1, z+1)
 774 \end_layout
 775
 776 \begin_layout Plain Layout
 777
 778 \end_layout
 779
 780 \begin_layout Plain Layout
 781
 782 sage: def FRD(x,y,z):
 783 \end_layout
 784
 785 \begin_layout Plain Layout
 786
 787 ....:     return grid(x-1, y+1, z-1)
 788 \end_layout
 789
 790 \begin_layout Plain Layout
 791
 792 \end_layout
 793
 794 \begin_layout Plain Layout
 795
 796 sage: def FRT(x,y,z):
 797 \end_layout
 798
 799 \begin_layout Plain Layout
 800
 801 ....:     return grid(x-1, y+1, z+1)
 802 \end_layout
 803
 804 \begin_layout Plain Layout
 805
 806 \end_layout
 807
 808 \begin_layout Plain Layout
 809
 810 sage: def BLD(x,y,z):
 811 \end_layout
 812
 813 \begin_layout Plain Layout
 814
 815 ....:     return grid(x+1, y-1, z-1)
 816 \end_layout
 817
 818 \begin_layout Plain Layout
 819
 820 \end_layout
 821
 822 \begin_layout Plain Layout
 823
 824 sage: def BLT(x,y,z):
 825 \end_layout
 826
 827 \begin_layout Plain Layout
 828
 829 ....:     return grid(x+1, y-1, z+1)
 830 \end_layout
 831
 832 \begin_layout Plain Layout
 833
 834 \end_layout
 835
 836 \begin_layout Plain Layout
 837
 838 sage: def BRD(x,y,z):
 839 \end_layout
 840
 841 \begin_layout Plain Layout
 842
 843 ....:     return grid(x+1, y+1, z-1)
 844 \end_layout
 845
 846 \begin_layout Plain Layout
 847
 848 \end_layout
 849
 850 \begin_layout Plain Layout
 851
 852 sage: def BRT(x,y,z):
 853 \end_layout
 854
 855 \begin_layout Plain Layout
 856
 857 ....:     return grid(x+1, y+1, z+1)
 858 \end_layout
 859
 860 \end_inset
 861
 862
 863 \end_layout
 864
 865 \begin_layout Example*
 866 Next, we define the coefficients for the cube centered on
 867 \begin_inset Formula $\left(1,1,1\right)$
 868 \end_inset
 869
 870  based on these directional functions.
 871 \end_layout
 872
 873 \begin_layout Example*
 874 \begin_inset listings
 875 inline false
 876 status open
 877
 878 \begin_layout Plain Layout
 879
 880 sage: c0030 = (1/8)*( I(1,1,1) + F(1,1,1) + L(1,1,1) + T(1,1,1) +
 881 \end_layout
 882
 883 \begin_layout Plain Layout
 884
 885                       LT(1,1,1) + FL(1,1,1) + FT(1,1,1) + FLT(1,1,1) )
 886 \end_layout
 887
 888 \begin_layout Plain Layout
 889
 890 \end_layout
 891
 892 \begin_layout Plain Layout
 893
 894 sage: c0003 = (1/8)*( I(1,1,1) + F(1,1,1) + R(1,1,1) + T(1,1,1) +
 895 \end_layout
 896
 897 \begin_layout Plain Layout
 898
 899                 RT(1,1,1) + FR(1,1,1) + FT(1,1,1) + FRT(1,1,1) )
 900 \end_layout
 901
 902 \begin_layout Plain Layout
 903
 904 \end_layout
 905
 906 \begin_layout Plain Layout
 907
 908 sage: c0021 = (5/24)*(I(1,1,1) + F(1,1,1) + T(1,1,1) + FT(1,1,1)) +
 909 \end_layout
 910
 911 \begin_layout Plain Layout
 912
 913               (1/24)*(L(1,1,1) + FL(1,1,1) + LT(1,1,1) + FLT(1,1,1))
 914 \end_layout
 915
 916 \begin_layout Plain Layout
 917
 918 \end_layout
 919
 920 \begin_layout Plain Layout
 921
 922 sage: c0012 = (5/24)*(I(1,1,1) + F(1,1,1) + T(1,1,1) + FT(1,1,1)) +
 923 \end_layout
 924
 925 \begin_layout Plain Layout
 926
 927               (1/24)*(R(1,1,1) + FR(1,1,1) + RT(1,1,1) + FRT(1,1,1))
 928 \end_layout
 929
 930 \begin_layout Plain Layout
 931
 932 \end_layout
 933
 934 \begin_layout Plain Layout
 935
 936 \end_layout
 937
 938 \begin_layout Plain Layout
 939
 940 \end_layout
 941
 942 \end_inset
 943
 944
 945 \end_layout
 946
 947 \begin_layout Example*
 948 We can see what the constant values are now:
 949 \end_layout
 950
 951 \begin_layout Example*
 952 \begin_inset listings
 953 inline false
 954 status open
 955
 956 \begin_layout Plain Layout
 957
 958 sage: c0030
 959 \end_layout
 960
 961 \begin_layout Plain Layout
 962
 963 17/8
 964 \end_layout
 965
 966 \begin_layout Plain Layout
 967
 968 \end_layout
 969
 970 \begin_layout Plain Layout
 971
 972 sage: c0003
 973 \end_layout
 974
 975 \begin_layout Plain Layout
 976
 977 27/8
 978 \end_layout
 979
 980 \begin_layout Plain Layout
 981
 982 \end_layout
 983
 984 \begin_layout Plain Layout
 985
 986 sage: c0021
 987 \end_layout
 988
 989 \begin_layout Plain Layout
 990
 991 61/24
 992 \end_layout
 993
 994 \begin_layout Plain Layout
 995
 996 \end_layout
 997
 998 \begin_layout Plain Layout
 999
1000 sage: c0012
1001 \end_layout
1002
1003 \begin_layout Plain Layout
1004
1005 71/24
1006 \end_layout
1007
1008 \end_inset
1009
1010
1011 \end_layout
1012
1013 \end_body
1014 \end_document