84:
74:
53:
1337:
22:
785:
1332:{\displaystyle b={\begin{bmatrix}-dx^{2}g_{22}+u_{12}+u_{21}\\-dx^{2}g_{32}+u_{31}~~~~~~~~\\-dx^{2}g_{42}+u_{52}+u_{41}\\-dx^{2}g_{23}+u_{13}~~~~~~~~\\-dx^{2}g_{33}~~~~~~~~~~~~~~~~\\-dx^{2}g_{43}+u_{53}~~~~~~~~\\-dx^{2}g_{24}+u_{14}+u_{25}\\-dx^{2}g_{34}+u_{35}~~~~~~~~\\-dx^{2}g_{44}+u_{54}+u_{45}\\\end{bmatrix}}}
221:
I might not be the right person to address some of the things mentioned above which is beyond what I have seen with this subject. For instance, in terms of the eignevalues of this system, I am not aware if there is an expression that easily gives them. I don't see anything mentioned in my numerical
710:
Recalculating the approach given on the wiki page, I realized that there is a sign flaw regarding the right-hand-side of the equation system. Since the block tridiagonal notion of A uses signs flipped with respect to the discretetized 2 dimensional
Poisson equation (given in the first formula on the
1709:
I really don't want to get into the middle of this, but I would like to re-emphasize one point Sławomir made: for scientific and mathematical articles in
Knowledge, the general rule is general reference in cases like this. One or a small number of references are given that cover a large portion of
626:
This is a good example of what I mean. I have gone though my numerical methods books as well as papers I Xeroxed out of some journals, and I don't see this expression. It seems to me having all the eigenvalues should be very helpful in solving the
Poisson equation, so I am wondering if there is a
1710:
the article, with only a few inline citations on individual claims (where material is not from those general sources). This is unlike political or culture articles where there is not a general reference for the article and so many/most sentences are individually cited.
1611:
706:
I fortunately found this page, when I was trying to implement a solver for a boundary value problem in image processing discipline. I noted, however, that the solutions produced by the solver showed an oposite sign as expected.
778:
140:
603:
235:
There is no expression for the eigenvalues of the discrete poisson equation for arbitrary domains, but over a rectangular or square grid, with uniform spacing, it is pretty simple. Let
1399:
412:
157:
It would be nice to expand this page a bit so that it has real information, not just as a page that gets people "started in the right direction". Things that would be nice to add:
1688:
a system of equations. It doesn't mean that any article with the phrase "how to solve" should be decimated. Rather it means that we don't provide instruction manuals. That is
490:
1677:
original research. Anyone who argues otherwise either doesn't know enough to make that assessment, or is not acting in good faith: this is completely standard material.
1647:
I have reverted the blanking of the article a number of times. First off, not every sentence needs a citation. In fact, for standard material like that in the section
739:
198:
550:
530:
510:
452:
432:
253:
1784:
130:
1659:
does not demand that each and every statement be cited individually, and I think that is clearly the case here. We don't even delete material that is
1779:
1756:
The matrix A for poisson equation is written wrongly. The diagonal (i,j+1(where i=j)) should be -1 for whole matrix and same with (i,j-1). Dekay315
1732:
I agree with the comment above and also wanted to emphasize another point. The section on Method of solution, which keeps getting removed per
106:
1623:
Without the negative in front, what is on the page after BR's correction is correct. I suspect that there was a typo that lead to the problem.
222:
methods books. I have seen discussion of FFT as a solution method, but I want to apply it before I am comfortable elaborating more on it here.
552:
are small, this is close to the spectrum of the continuous laplacian. You can derive this if you assume the eigenfunctions are in the form
1629:
164:
for the case of a square grid and disc, mention or derive the fast poisson methods (those involving the FFT; see, eg. Arieh
Iserles' book).
674:
This article is pretty poor IMHO. How about explaining what it is without resorting to algebra, and also explaining its applications... --
1722:
97:
58:
1673:, in hopes that someone will provide a citation. We don't have a scorched-earth policy towards uncited material, and this is
1699:
685:
33:
648:
In this page the A is written wrongly. The diagonal ( i,j+1 (when i=j) ) should be -1 for all. And same with (i,j-1).
688:, which is the only place I have used this discretization. I don't think I can do anything about the algebra though.
747:
1684:
doesn't mean that we shouldn't include a discussion of methods for solution of a system in an encyclopedia article
1606:{\displaystyle -({\nabla }^{2}u)_{ij}=-{\frac {1}{dx^{2}}}(u_{i+1,j}+u_{i-1,j}+u_{i,j+1}+u_{i,j-1}-4u_{ij})=g_{ij}}
1648:
1693:
555:
255:
be the number of interior grid points, and let the domain be the unit grid, all the eigenvalues are in the form
1667:
1633:
21:
1736:
is not a how to. It neutrally lists multiple different algorithms that might be used to solve this system.
260:
1718:
1761:
206:
perhaps a mention to the form of the discrete
Laplacian, e.g. block TST (Toeplitz Symmetric Tridiagonal)
39:
83:
1378:
1351:
654:
1757:
1741:
649:
457:
105:
on
Knowledge. If you would like to participate, please visit the project page, where you can join
1347:
689:
628:
223:
89:
627:
method out there that takes advantage of the above. Obviously, I have to research this further.
73:
52:
1713:
1656:
714:
173:
mention the discrete laplace equation and how it is simpler (and how the 5 pt stencil gets
176:
1737:
170:
various stencils and their accuracy, e.g. standard 5 pt, and the 9 pt and modified 9 pt
1765:
1745:
1727:
1703:
1637:
1382:
1355:
692:
678:
658:
631:
609:
535:
515:
495:
437:
417:
238:
226:
214:
1773:
1733:
1681:
1663:
verifiable, that lacks in citations. If necessary, we tag it with a template like
675:
1652:
1374:
606:
211:
102:
79:
711:
page) - the constant vector b has to use oposite signs on the derivatives (
161:
for the case of a square grid, derive the condition number of the matrix
203:
make the matrices a bit more general, using kronecker product notation
1651:, a general reference should be enough, as long as the material is
684:
I added an "Applications" section for where it is encountered in
167:
for the case of a square grid, the eignevalues and eigenfunctions
1649:
Discrete
Poisson equation#On a two-dimensional rectangular grid
1368:,Brooks/Cole Publishing Company, Pacific Grove, 1985, page 448
15:
1373:
and you are right. How did no one notice this until now?
1394:
Generally, the
Discrete Poisson Equation takes the form
800:
756:
1402:
788:
750:
717:
558:
538:
518:
498:
460:
440:
420:
263:
241:
179:
101:, a collaborative effort to improve the coverage of
1605:
1331:
772:
733:
597:
544:
524:
504:
484:
446:
426:
406:
247:
192:
773:{\displaystyle {\begin{bmatrix}b\end{bmatrix}}}
744:So, the correct notion for the constant vector
200:accuracy, similar to the modified 9 pt stencil)
1341:I changed the signs on the page accordingly.
8:
1366:Numerical Mathematics and Computing 2nd Ed.
598:{\displaystyle \sin(\alpha x)\sin(\beta y)}
19:
652:
47:
1594:
1575:
1547:
1522:
1497:
1472:
1456:
1443:
1428:
1415:
1410:
1401:
1315:
1302:
1289:
1279:
1243:
1230:
1220:
1200:
1187:
1174:
1164:
1128:
1115:
1105:
1053:
1043:
1007:
994:
984:
964:
951:
938:
928:
892:
879:
869:
849:
836:
823:
813:
795:
787:
751:
749:
722:
716:
670:less algebra, and more explanation please
557:
537:
517:
497:
459:
439:
419:
364:
352:
310:
298:
268:
262:
240:
184:
178:
49:
407:{\displaystyle \lambda _{a,b}=-4\left}
7:
95:This article is within the scope of
38:It is of interest to the following
1411:
14:
1785:Low-priority mathematics articles
115:Knowledge:WikiProject Mathematics
1780:Start-Class mathematics articles
118:Template:WikiProject Mathematics
82:
72:
51:
20:
1364:Cheny, Ward and David Kincaid,
135:This article has been rated as
1584:
1465:
1425:
1406:
592:
583:
574:
565:
396:
390:
378:
361:
342:
336:
324:
307:
215:02:30, 18 September 2006 (UTC)
1:
1746:07:02, 22 November 2011 (UTC)
1728:00:19, 21 November 2011 (UTC)
1704:22:37, 20 November 2011 (UTC)
1383:17:37, 13 February 2008 (UTC)
1356:13:32, 13 February 2008 (UTC)
693:20:46, 19 November 2006 (UTC)
679:10:21, 19 November 2006 (UTC)
632:23:18, 20 November 2006 (UTC)
610:18:02, 20 November 2006 (UTC)
485:{\displaystyle 1,2,\cdots ,m}
227:20:46, 19 November 2006 (UTC)
109:and see a list of open tasks.
686:Computational fluid dynamics
1801:
1766:12:04, 15 June 2018 (UTC)
1361:I checked the reference:
659:12:01, 15 June 2018 (UTC)
134:
67:
46:
699:Correct Algebraic Signs?
141:project's priority scale
1638:19:21, 6 May 2009 (UTC)
98:WikiProject Mathematics
1607:
1333:
774:
735:
734:{\displaystyle g_{ij}}
599:
546:
526:
506:
486:
448:
428:
408:
249:
194:
28:This article is rated
1608:
1334:
775:
736:
600:
547:
527:
507:
487:
449:
429:
409:
250:
195:
193:{\displaystyle h^{4}}
1400:
786:
748:
715:
556:
536:
516:
496:
458:
438:
418:
261:
239:
177:
121:mathematics articles
1603:
1329:
1323:
770:
764:
731:
595:
542:
522:
502:
482:
444:
424:
404:
245:
190:
90:Mathematics portal
34:content assessment
1726:
1463:
1264:
1262:
1260:
1258:
1256:
1254:
1252:
1250:
1149:
1147:
1145:
1143:
1141:
1139:
1137:
1135:
1090:
1088:
1086:
1084:
1082:
1080:
1078:
1076:
1074:
1072:
1070:
1068:
1066:
1064:
1062:
1060:
1028:
1026:
1024:
1022:
1020:
1018:
1016:
1014:
913:
911:
909:
907:
905:
903:
901:
899:
661:
545:{\displaystyle b}
525:{\displaystyle a}
505:{\displaystyle m}
447:{\displaystyle b}
427:{\displaystyle a}
394:
340:
248:{\displaystyle m}
155:
154:
151:
150:
147:
146:
1792:
1716:
1696:
1672:
1666:
1612:
1610:
1609:
1604:
1602:
1601:
1583:
1582:
1564:
1563:
1539:
1538:
1514:
1513:
1489:
1488:
1464:
1462:
1461:
1460:
1444:
1436:
1435:
1420:
1419:
1414:
1338:
1336:
1335:
1330:
1328:
1327:
1320:
1319:
1307:
1306:
1294:
1293:
1284:
1283:
1263:
1261:
1259:
1257:
1255:
1253:
1251:
1249:
1248:
1247:
1235:
1234:
1225:
1224:
1205:
1204:
1192:
1191:
1179:
1178:
1169:
1168:
1148:
1146:
1144:
1142:
1140:
1138:
1136:
1134:
1133:
1132:
1120:
1119:
1110:
1109:
1089:
1087:
1085:
1083:
1081:
1079:
1077:
1075:
1073:
1071:
1069:
1067:
1065:
1063:
1061:
1059:
1058:
1057:
1048:
1047:
1027:
1025:
1023:
1021:
1019:
1017:
1015:
1013:
1012:
1011:
999:
998:
989:
988:
969:
968:
956:
955:
943:
942:
933:
932:
912:
910:
908:
906:
904:
902:
900:
898:
897:
896:
884:
883:
874:
873:
854:
853:
841:
840:
828:
827:
818:
817:
779:
777:
776:
771:
769:
768:
740:
738:
737:
732:
730:
729:
604:
602:
601:
596:
551:
549:
548:
543:
531:
529:
528:
523:
511:
509:
508:
503:
491:
489:
488:
483:
453:
451:
450:
445:
433:
431:
430:
425:
413:
411:
410:
405:
403:
399:
395:
393:
373:
365:
357:
356:
341:
339:
319:
311:
303:
302:
279:
278:
254:
252:
251:
246:
199:
197:
196:
191:
189:
188:
123:
122:
119:
116:
113:
92:
87:
86:
76:
69:
68:
63:
55:
48:
31:
25:
24:
16:
1800:
1799:
1795:
1794:
1793:
1791:
1790:
1789:
1770:
1769:
1754:
1694:
1670:
1668:citation needed
1664:
1645:
1590:
1571:
1543:
1518:
1493:
1468:
1452:
1448:
1424:
1409:
1398:
1397:
1322:
1321:
1311:
1298:
1285:
1275:
1266:
1265:
1239:
1226:
1216:
1207:
1206:
1196:
1183:
1170:
1160:
1151:
1150:
1124:
1111:
1101:
1092:
1091:
1049:
1039:
1030:
1029:
1003:
990:
980:
971:
970:
960:
947:
934:
924:
915:
914:
888:
875:
865:
856:
855:
845:
832:
819:
809:
796:
784:
783:
763:
762:
752:
746:
745:
718:
713:
712:
701:
672:
554:
553:
534:
533:
514:
513:
494:
493:
456:
455:
436:
435:
416:
415:
374:
366:
348:
320:
312:
294:
293:
289:
264:
259:
258:
237:
236:
180:
175:
174:
120:
117:
114:
111:
110:
88:
81:
61:
32:on Knowledge's
29:
12:
11:
5:
1798:
1796:
1788:
1787:
1782:
1772:
1771:
1753:
1750:
1749:
1748:
1730:
1711:
1695:Sławomir Biały
1692:what this is.
1644:
1641:
1630:98.201.170.218
1627:
1626:
1625:
1624:
1618:
1617:
1616:
1615:
1614:
1613:
1600:
1597:
1593:
1589:
1586:
1581:
1578:
1574:
1570:
1567:
1562:
1559:
1556:
1553:
1550:
1546:
1542:
1537:
1534:
1531:
1528:
1525:
1521:
1517:
1512:
1509:
1506:
1503:
1500:
1496:
1492:
1487:
1484:
1481:
1478:
1475:
1471:
1467:
1459:
1455:
1451:
1447:
1442:
1439:
1434:
1431:
1427:
1423:
1418:
1413:
1408:
1405:
1388:
1386:
1385:
1371:
1370:
1369:
1326:
1318:
1314:
1310:
1305:
1301:
1297:
1292:
1288:
1282:
1278:
1274:
1271:
1268:
1267:
1246:
1242:
1238:
1233:
1229:
1223:
1219:
1215:
1212:
1209:
1208:
1203:
1199:
1195:
1190:
1186:
1182:
1177:
1173:
1167:
1163:
1159:
1156:
1153:
1152:
1131:
1127:
1123:
1118:
1114:
1108:
1104:
1100:
1097:
1094:
1093:
1056:
1052:
1046:
1042:
1038:
1035:
1032:
1031:
1010:
1006:
1002:
997:
993:
987:
983:
979:
976:
973:
972:
967:
963:
959:
954:
950:
946:
941:
937:
931:
927:
923:
920:
917:
916:
895:
891:
887:
882:
878:
872:
868:
864:
861:
858:
857:
852:
848:
844:
839:
835:
831:
826:
822:
816:
812:
808:
805:
802:
801:
799:
794:
791:
767:
761:
758:
757:
755:
728:
725:
721:
700:
697:
696:
695:
671:
668:
667:
666:
665:
664:
663:
662:
657:comment added
641:
640:
639:
638:
637:
636:
635:
634:
617:
616:
615:
614:
613:
612:
594:
591:
588:
585:
582:
579:
576:
573:
570:
567:
564:
561:
541:
521:
501:
481:
478:
475:
472:
469:
466:
463:
443:
423:
402:
398:
392:
389:
386:
383:
380:
377:
372:
369:
363:
360:
355:
351:
347:
344:
338:
335:
332:
329:
326:
323:
318:
315:
309:
306:
301:
297:
292:
288:
285:
282:
277:
274:
271:
267:
244:
230:
229:
218:
217:
208:
207:
204:
201:
187:
183:
171:
168:
165:
162:
153:
152:
149:
148:
145:
144:
133:
127:
126:
124:
107:the discussion
94:
93:
77:
65:
64:
56:
44:
43:
37:
26:
13:
10:
9:
6:
4:
3:
2:
1797:
1786:
1783:
1781:
1778:
1777:
1775:
1768:
1767:
1763:
1759:
1751:
1747:
1743:
1739:
1735:
1731:
1729:
1724:
1720:
1715:
1712:
1708:
1707:
1706:
1705:
1701:
1697:
1691:
1687:
1683:
1678:
1676:
1669:
1662:
1658:
1654:
1650:
1642:
1640:
1639:
1635:
1631:
1622:
1621:
1620:
1619:
1598:
1595:
1591:
1587:
1579:
1576:
1572:
1568:
1565:
1560:
1557:
1554:
1551:
1548:
1544:
1540:
1535:
1532:
1529:
1526:
1523:
1519:
1515:
1510:
1507:
1504:
1501:
1498:
1494:
1490:
1485:
1482:
1479:
1476:
1473:
1469:
1457:
1453:
1449:
1445:
1440:
1437:
1432:
1429:
1421:
1416:
1403:
1396:
1395:
1393:
1392:
1391:
1390:
1389:
1384:
1380:
1376:
1372:
1367:
1363:
1362:
1360:
1359:
1358:
1357:
1353:
1349:
1345:
1342:
1339:
1324:
1316:
1312:
1308:
1303:
1299:
1295:
1290:
1286:
1280:
1276:
1272:
1269:
1244:
1240:
1236:
1231:
1227:
1221:
1217:
1213:
1210:
1201:
1197:
1193:
1188:
1184:
1180:
1175:
1171:
1165:
1161:
1157:
1154:
1129:
1125:
1121:
1116:
1112:
1106:
1102:
1098:
1095:
1054:
1050:
1044:
1040:
1036:
1033:
1008:
1004:
1000:
995:
991:
985:
981:
977:
974:
965:
961:
957:
952:
948:
944:
939:
935:
929:
925:
921:
918:
893:
889:
885:
880:
876:
870:
866:
862:
859:
850:
846:
842:
837:
833:
829:
824:
820:
814:
810:
806:
803:
797:
792:
789:
781:
765:
759:
753:
742:
726:
723:
719:
708:
704:
698:
694:
691:
687:
683:
682:
681:
680:
677:
669:
660:
656:
651:
647:
646:
645:
644:
643:
642:
633:
630:
625:
624:
623:
622:
621:
620:
619:
618:
611:
608:
589:
586:
580:
577:
571:
568:
562:
559:
539:
519:
499:
492:. For large
479:
476:
473:
470:
467:
464:
461:
454:ranging from
441:
421:
400:
387:
384:
381:
375:
370:
367:
358:
353:
349:
345:
333:
330:
327:
321:
316:
313:
304:
299:
295:
290:
286:
283:
280:
275:
272:
269:
265:
257:
256:
242:
234:
233:
232:
231:
228:
225:
220:
219:
216:
213:
210:
209:
205:
202:
185:
181:
172:
169:
166:
163:
160:
159:
158:
142:
138:
132:
129:
128:
125:
108:
104:
100:
99:
91:
85:
80:
78:
75:
71:
70:
66:
60:
57:
54:
50:
45:
41:
35:
27:
23:
18:
17:
1755:
1734:WP:NOT#HOWTO
1714:CRGreathouse
1689:
1685:
1682:WP:NOT#HOWTO
1679:
1674:
1661:in principle
1660:
1646:
1628:
1387:
1365:
1346:
1343:
1340:
782:
743:
709:
705:
702:
673:
156:
137:Low-priority
136:
96:
62:Low‑priority
40:WikiProjects
741:) as well.
653:—Preceding
112:Mathematics
103:mathematics
59:Mathematics
30:Start-class
1774:Categories
1752:Correction
1680:Secondly,
1657:WP:SCICITE
1653:verifiable
780:would be:
1738:Thenub314
512:and when
1758:Dekay315
1643:Blanking
650:Dekay315
676:Rebroad
655:undated
139:on the
1375:CFDFEM
1348:Smader
690:Slffea
629:Slffea
607:Lavaka
224:Slffea
212:Lavaka
36:scale.
1686:about
1762:talk
1742:talk
1700:talk
1634:talk
1379:talk
1352:talk
703:Hi,
532:and
434:and
414:for
1690:not
1675:not
1655:.
1344:BR
578:sin
560:sin
350:sin
296:sin
131:Low
1776::
1764:)
1744:)
1721:|
1702:)
1671:}}
1665:{{
1636:)
1566:−
1558:−
1502:−
1441:−
1412:∇
1404:−
1381:)
1354:)
1317:45
1304:54
1291:44
1270:−
1245:35
1232:34
1211:−
1202:25
1189:14
1176:24
1155:−
1130:53
1117:43
1096:−
1055:33
1034:−
1009:13
996:23
975:−
966:41
953:52
940:42
919:−
894:31
881:32
860:−
851:21
838:12
825:22
804:−
605:.
587:β
581:
569:α
563:
474:⋯
371:π
359:
317:π
305:
284:−
266:λ
1760:(
1740:(
1725:)
1723:c
1719:t
1717:(
1698:(
1632:(
1599:j
1596:i
1592:g
1588:=
1585:)
1580:j
1577:i
1573:u
1569:4
1561:1
1555:j
1552:,
1549:i
1545:u
1541:+
1536:1
1533:+
1530:j
1527:,
1524:i
1520:u
1516:+
1511:j
1508:,
1505:1
1499:i
1495:u
1491:+
1486:j
1483:,
1480:1
1477:+
1474:i
1470:u
1466:(
1458:2
1454:x
1450:d
1446:1
1438:=
1433:j
1430:i
1426:)
1422:u
1417:2
1407:(
1377:(
1350:(
1325:]
1313:u
1309:+
1300:u
1296:+
1287:g
1281:2
1277:x
1273:d
1241:u
1237:+
1228:g
1222:2
1218:x
1214:d
1198:u
1194:+
1185:u
1181:+
1172:g
1166:2
1162:x
1158:d
1126:u
1122:+
1113:g
1107:2
1103:x
1099:d
1051:g
1045:2
1041:x
1037:d
1005:u
1001:+
992:g
986:2
982:x
978:d
962:u
958:+
949:u
945:+
936:g
930:2
926:x
922:d
890:u
886:+
877:g
871:2
867:x
863:d
847:u
843:+
834:u
830:+
821:g
815:2
811:x
807:d
798:[
793:=
790:b
766:]
760:b
754:[
727:j
724:i
720:g
593:)
590:y
584:(
575:)
572:x
566:(
540:b
520:a
500:m
480:m
477:,
471:,
468:2
465:,
462:1
442:b
422:a
401:]
397:)
391:)
388:1
385:+
382:m
379:(
376:2
368:b
362:(
354:2
346:+
343:)
337:)
334:1
331:+
328:m
325:(
322:2
314:a
308:(
300:2
291:[
287:4
281:=
276:b
273:,
270:a
243:m
186:4
182:h
143:.
42::
Text is available under the Creative Commons Attribution-ShareAlike License. Additional terms may apply.