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

Timestamp:

01/15/15 12:22:32 (10 years ago)

Author:

Erica Kaminski

Comment:

—

u/erica/CylindricalGravity

@@ -113 +113 @@
 Dx for the patch is easy, that is just,
 
-dx = levels(Info%level)% dx
+[[latex($dx = levels(Info\%level)\% dx$)]]
 
 To get all of the x_i's of the cells in the patch we can allocate an array 'x_i' that has length = number of cells along x axis of patch,
 
-Allocate x(Info%mGlobal(1,1): Info%mGlobal(1,2))
+[[latex($Allocate ~x(Info\%mGlobal(1,1): Info\%mGlobal(1,2))$)]]
 
 (For a refresher page describing patch specific quantities, refer here)
@@ -123 +123 @@
 Now to populate this array, we loop through:
 
-Do i = Info%mGlobal(1,1), Info%mGlobal(1,2)[[br]]
- x(i) = GxBounds(1,1) + (i - half)*dx[[br]]
-End Do
+[[latex($Do ~i = Info\%mGlobal(1,1), Info\%mGlobal(1,2)$)]]
+ [[latex($~x(i) = GxBounds(1,1) + (i - half)*dx$)]]
+[[latex($End ~Do$)]]
 
 This starts at the  ''global'' lower left hand corner of the entire box (gXbounds(1,1)) and marches over i*dx_patch spaces (subtracting off the .5dx to get to the cell center). 
@@ -134 +134 @@
 Once we have these 2 things calculated, we can then add in the meat and potatoes: 
 
-If (iCylindrical > 0) Then [[br]]
- tempStencilValues(1) = 1-dx/2x(i) [[br]]
- tempStencilValues(2) = 1+dx/2x(i) [[br]]
-End If [[br]]
+[[latex($If ~(iCylindrical > 0)~Then$)]]
+ [[latex($~tempStencilValues(1) = 1-dx/2x(i)$)]] [[br]]
+ [[latex($tempStencilValues(2) = 1+dx/2x(i)$)]]
+[[latex($End~ If$)]][[br]]
 
 Then set the matrixvalues for the stencil of the ith cell, by doing
 
-matrixvalues(m:m+4) = tempStencilValues(1:5)
+[[latex($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)
+[[latex($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.