Stephen Mardson McQuay
456ecf759f
added one more 2D gmsh test case
2010-11-01 15:52:45 -06:00
Stephen Mardson McQuay
37f056e49b
added 3D gmsh case
2010-11-01 15:52:03 -06:00
Stephen Mardson McQuay
17b22c173d
added a 3d gmsh test script
2010-11-01 15:51:40 -06:00
Stephen Mardson McQuay
81f1c6557a
added interactive prompt bootstraping code
2010-11-01 15:51:16 -06:00
Stephen Mardson McQuay
3aefd7b75c
made one test approximate
2010-11-01 15:50:54 -06:00
Stephen Mardson McQuay
656394e747
updated for directory structure changes
2010-10-29 21:34:15 -06:00
Stephen Mardson McQuay
22a8855a99
removed the old baker error routines
2010-10-29 21:05:21 -06:00
Stephen Mardson McQuay
bcf1622f35
added a generic get simplex and points function to the grid
2010-10-29 12:43:18 -06:00
Stephen Mardson McQuay
57685a1c04
smblog -> log ... again
2010-10-29 12:42:05 -06:00
Stephen Mardson McQuay
e9fb1eab03
the gmsh grid object now instantiates its own KDTree
2010-10-29 12:41:39 -06:00
Stephen Mardson McQuay
967170d765
changed points to verts in baker
2010-10-29 12:40:32 -06:00
Stephen Mardson McQuay
699d0c1796
started testing the parametric varying of the solution scheme
2010-10-29 12:39:45 -06:00
fe1564881a
more refactoring
2010-10-23 17:06:57 -06:00
7d32f4391f
updated imports, and changed smberror to Exception
2010-10-23 17:06:11 -06:00
02b311f72d
removed the smberror, and corrected a typo in the exact function
2010-10-23 17:05:33 -06:00
47a840c4b4
updated the interp.grid object to be a better base class
2010-10-23 17:01:10 -06:00
2eaec8268c
created a distince gmsh 2D plugin
2010-10-23 16:57:44 -06:00
71761dd5a1
removed uneccesarry import from blender plotting file
2010-10-23 16:56:52 -06:00
885c5c34c1
beefed up the setup file so that people can find me
2010-10-23 16:56:30 -06:00
eadfa46e57
middle of refactor
2010-10-23 13:49:15 -06:00
431f3a4e44
forgot the __init__.py file
2010-10-22 16:21:41 -06:00
2db4169bfa
major refactoring
...
--HG--
rename : lib/baker/__init__.py => interp/baker/__init__.py
rename : lib/grid/DD.py => interp/grid/DD.py
rename : lib/grid/DDD.py => interp/grid/DDD.py
rename : lib/grid/__init__.py => interp/grid/__init__.py
rename : lib/grid/qhull.py => interp/grid/qhull.py
rename : lib/grid/simplex.py => interp/grid/simplex.py
rename : lib/grid/smcqdelaunay.py => interp/grid/smcqdelaunay.py
rename : lib/baker/tools.py => interp/tools.py
2010-10-22 16:10:58 -06:00
98d7247497
minor: doc change
2010-10-22 09:21:04 -06:00
91c8a661f4
working 2D gmsh parser
2010-10-22 09:20:59 -06:00
e4988fd150
implementing a 2D mesh object for gmsh files
2010-10-06 12:15:09 -06:00
b3680848d5
testing push
2010-09-04 16:52:16 -06:00
73322f7554
testing push
2010-09-04 16:51:05 -06:00
df35cb174b
made a get_error function that is generic and functional. I'm still getting unsavory numbers for the cubic, 2-D case on my single test case. I will have to try out more refined grids, albeit the quadratic is consistently improving results for my case.
2010-05-05 23:03:12 -06:00
3a1c13bcac
implemented a function that will be used to provide a generic n-th order, nth-dimension error approximation function
2010-05-04 23:03:07 -06:00
Stephen Mardson McQuay
c5adab295b
minor: but i'm on the verge of implementing the cubic error term for 3D
2010-04-30 16:23:29 -06:00
Stephen Mardson McQuay
a1c88a186e
something about the changes I made (abs -> np.abs) made my stuff succeed more often
2010-04-30 11:56:01 -06:00
2cbd92e15b
wrapping up a night. there isn't enough consistent improvement to merit using this method. i must have a bug somewhere.
2010-04-29 23:29:35 -06:00
Stephen Mardson McQuay
3dcc10ea0e
working on the containing simplex problem. the get containing simplex function works, if the point is in the domain. I am trying to force the domain to contain the point (placing points along perimiter of the mesh). I am not doing that right, but have to run
...
--HG--
rename : bin/qhull-029.txt => data/qhull-029.txt
rename : bin/qhull-029.txt.gv => data/qhull-029.txt.gv
2010-04-24 18:26:44 -06:00
Stephen Mardson McQuay
cfcd1e15a2
working on the guarantee simplex routine. I have a test case to try it out against, and need to look at visiting neighbors of faces adjacent to points.
2010-04-23 16:29:31 -06:00
Stephen Mardson McQuay
0a388ff1b5
trying to guarantee a containing simplex
2010-04-23 09:58:30 -06:00
Stephen Mardson McQuay
79880acad4
fixed discrepancy with Baker's eq 17.
2010-04-09 15:03:16 -06:00
Stephen Mardson McQuay
3efe335563
moved the quadratic error calculation out of the baker method, and added a cubic version too
2010-04-02 11:06:17 -06:00
Stephen Mardson McQuay
700ccc8c25
created a test for tesing the cubic interpolation
2010-03-20 20:09:46 -06:00
Stephen Mardson McQuay
1fc988428d
putting together a simple test case so that I can test my quad/cubic interpolator
2010-03-20 17:10:02 -06:00
Stephen Mardson McQuay
993c901cec
properly catching the LinAlgError. also, prepping to change the run_baker method (more generic)
2010-03-20 11:34:24 -06:00
Stephen Mardson McQuay
bd8f24d64f
working on implementing the cubic interpolation. need to split that routine out of the run_baker method
2010-03-19 21:27:51 -06:00
Stephen Mardson McQuay
f089848a1b
i think part of my problem (why i'm not always getting win), is that i do not guarantee that the simplex contains the point (especially with the random clounds; i think that the connectivity-based lookup is doing that fine
2010-03-18 23:50:40 -06:00
Stephen Mardson McQuay
9836fdcbfe
added files. still not consistently getting better numbers
2010-03-18 23:32:24 -06:00
Stephen Mardson McQuay
b33159f8a9
working on getting baker to work in 3D. the code runs, but the numbers are odd. I suspect the rectangular grid, and am going to try a random cloud of points.
2010-03-18 23:18:59 -06:00
Stephen Mardson McQuay
a2d7b3f063
updated library files so that i don't have to edit the __init__.py files
...
--HG--
rename : lib/baker/__init__.py => lib/baker/baker.py
rename : lib/grid/__init__.py => lib/grid/grid.py
2010-03-08 14:30:22 -07:00
Stephen Mardson McQuay
b9ea6a3ac2
modified code to use Distribute (python's distribution framework)
...
one should simply run `python setup.py install` to install these librarires
2010-03-08 13:31:51 -07:00
smcquay@cfdviz2
c5d8fff883
moved some files around to conform to python best practicies
...
--HG--
rename : lib/baker.py => lib/baker/__init__.py
rename : lib/tools.py => lib/baker/tools.py
rename : lib/grid.py => lib/grid/__init__.py
rename : lib/smcqdelaunay.py => lib/grid/smcqdelaunay.py
2010-03-08 13:05:42 -07:00
Stephen Mardson McQuay
9a1b8d14b2
finally added proper (in my opinon) exception handling to the grid's run_baker method. it should properly autocorrect, and switch the extra point lookup scheme to connectivity-based if solving the system of equations yields a singular solution
2010-03-05 08:58:07 -07:00
smcquay@cfdviz2
a25e6d03d1
partially implemented using new data structure in point lookup
...
I spoke with Dr. Gorrell about first trying to call the nearest-neighbor routine, releasing the trapped exception where the matrix ends up being singular, and trapping it in the grid object's run_baker, where I could then call the connectivity-based lookup routine. Those are my next steps, as well as starting to think about pushing the 3D front ahead this week.
2010-02-26 14:44:29 -07:00
Stephen Mardson McQuay
2792329eaa
cleaned up the method (removed print statements) i'm ready to now add this method as an option to the mesh object
2010-02-21 20:52:51 -07:00
