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

Timestamp:

01/13/15 17:44:15 (10 years ago)

Author:

Erica Kaminski

Comment:

—

u/erica/CylindricalGravity

@@ -109 +109 @@
 ==== 2.5D Mod ====
 
-What we then need to add is a modification in these loops for cylindrical. This modification is going to need both the dx of the patch, and the x position of the cell (also patch dependent). Instead of calculating for each iteration of the loop, we can try and save some computation time and put the calculation of x and dx before the loop. 
+What we then need to add is a modification in these loops for cylindrical. As stated above, this modification is going to need both the dx of the patch, and the x position of the cell (also patch dependent - because x(dx)). Instead of calculating for each iteration of the loop, we can try and save some computational power and put the calculation of x and dx before the loop. 
 
 Dx for the patch is easy, that is just,
 
-''' levels(Info%level)% dx'''
+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))'''
+Allocate x(Info%mGlobal(1,1): Info%mGlobal(1,2))
 
-For a page describing patch specific quantities, refer here.
+(For a refresher page describing patch specific quantities, refer here)
 
 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'''
+Do i = Info%mGlobal(1,1), Info%mGlobal(1,2)[[br]]
+ x(i) = GxBounds(1,1) + (i - half)*dx[[br]]
+End Do
 
-This starts at the lower left hand corner of the mesh (gXbounds(1,1)) and marches over i*dx_patch times (subtracting off the .5dx to get to the cell center). 
+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). 
 
 Another way to do this would have been to use the lower x position of the patch itself, say:
 
+
+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]]
+
+Then set the matrixvalues for the stencil of the ith cell = tempStencilValues