Changes between Version 25 and Version 26 of u/erica/CylindricalGravity

Timestamp:

01/15/15 12:17:26 (10 years ago)

Author:

Erica Kaminski

Comment:

—

u/erica/CylindricalGravity

@@ -139 +139 @@
 End If [[br]]
 
-Then set the matrixvalues for the stencil of the ith cell = tempStencilValues
+Then set the matrixvalues for the stencil of the ith cell, by doing
+
+matrixvalues(m:m+4) = tempStencilValues(1:5)
+
+where m increases by 5 for each cell in the mesh. For example, 
+
+matrixvalues(:)=(-4,1,1,1,1, -4, 1, 1, 1, 1, -4, 1, 1, 1, 1, ...., -4, 1, 1, 1, 1)
+
+for regular 2D Cartesian geometry. 
+
+This means that matrixvalues(1) is the center cell stencil coefficient of the 1st cell in the mesh, matrix(6) is the center cell stencil coefficient of the 2nd cell in the mesh, and so on.