Parent Directory
|
Revision Log
Links to HEAD: | (view) (annotate) |
Sticky Revision: |
added missing ifdefs for library presence
added ripley rectangle reading from compressed files, added unit tests for binaryGrid reads and writes (only single ranks for compressed thus far) (#104)
added exceptions for conflicting assembler types, added tests for lame assemblers
vastly decreased ripley domain setup time under MPI (approximately from O(n*m) to O(n+m), n being shared nodes, m being the rank DOF)
adding skeleton of fast Lame assemblers adjusted automatic ripley domain subdivision ( should solve issue #94 ) more tab->space conversions
Set the randomfill back to generating randoms and updated doco to tell people how to use it.
This version has been tested in some detail. It does not actually do randoms though (since that wrecks the testing). In will do a commit soon which puts the randomness back
Correcting python libname for sage. Fixing bug where the number of values in the shape was not considered for buffer and message size.
Remove randomFill python method from ripley domains. All random data objects (for all domain types) should be generated using esys.escript.RandomData() The only filtered random we have is gaussian on ripley but it is triggered by passing the tuple as the last arg of RandomData(). While the interface is a bit more complicated (in that you always need to pass in shape and functionspace) it does mean we have a common interface for all domains. Removed randomFill from DataExpanded. The reasoning behind this is to force domains to call the util function themselves and enforce whatever consistancy requirements they have. Added version of blocktools to deal with 2D case in Ripley. Use blocktools for the 2D transfers [This was cleaner than modifying the previous implementation to deal with variable shaped points]. Note that under MPI, ripley can not generate random data (even unfiltered) if any of its per rank dimensions is <4 elements on any side. Unit tests for these calls are in but some extra checks still needed.
expanded tab chars for consistentcy and removed trailing whitespace
fixed dirac point->node mapping issues in ripley when using mpi
I changed some files. Updated copyright notices, added GeoComp.
Fixed a spelling error and missing virtual (not sure that one matters). Added _untested_ 3D gaussian smoothed random data
Fixes necessary to get a clang/macports build working on 10.9
Fixed failure for lazy data when interpolating.
added dirac point interpolation to ripley (and fixed dirac points in Brick) added fast wave assembler for Brick and minor correction to fast wave assembler for Rectangle
Eliminated all const_cast<Data*> hacks in ripley and finley now that Data.getSampleDataRO returns a const pointer.
Added dirac support to ripley, added interface for custom assemblers for ripleydomains (also added the custom assembler for 2D VTI waves), changed synthetic_VTI example to use the new, faster custom assembler
Implemented reverse reading of grid data from netcdf files. Fixes #44.
Step 1 for #31: wrap parameters in an object.
merged 3.4.1 release changes into trunk
unused variable fixed
2D guassian smoothing of random data done. Added more memory non-leaks to the supp file
uint -> unsigned int
more buildbot commit testing
deliberate error to test buildbot
Added libpython version of suppressions to the aggressive file. Rectangle now passes correct values around. Random still to come.
Some more work on random Added a more agressive (And smaller) supressions file for python
Some more work on this
Some work on rectangle random MPI
Ripley can now read grids in non-native byte order. This is in preparation for Gocad Voxet import.
Single rank rectangle randoms.
Remove bool_t Part of random.
Added support for more input data types in ER Mapper files. Tests will follow soon.
More leak fixes and suppression file additions.
escript: use scoped_array instead of scoped_ptr ripley: netcdf returns a copy of an array so use scoped_array downunder: s/Sperical/Spherical/ scripts: added a few more suppressions
ripley work: Almost done - one more round in Brick left....
More ripley code cleaning/optimisation.
WIP: ripley hand-optimisation & further rules to generator -> drastically reduced number of constants -> compile time of Brick.cpp less than half of what it was on savanna -> additional runtime savings -> to be continued...
put back std::
templated writeBinaryGrid and a bit of cleanup.
For some reason std::isnan<float> is not available on NCI so changed this to ::isnan. Added options file for NCI vayu.
Bringing the changes from doubleplusgood branch. Can't merge directly because svn doesn't transfer changes to renamed files (mutter grumble).
Preparations for Voxet (GOCAD) binary grid output + code clean up. Changed node id's order to remove code duplication in weipa. More work to do...
Assorted spelling fixes.
Implemented claimed support for mixing data with different resolutions.
ripley's readBinaryGrid and readNcGrid now filter out NaN leaving the relevant destination data points untouched.
Round 1 of copyright fixes
Implemented netCDF grid reader analogous to binary grid reader in ripley.
Fixed a potential div by zero in ripley's automatic domain subdivider....
Fixed a memory leak (freeing paso coupler) in ripley.
First pass of updating copyright notices
Implemented binary grid reader in ripley which is now used in ER Mapper source. Interface/details likely to change but need to get this running now.
Added heuristics for number of domain subdivisions to ripley so ripley domains can now be created without specifying the values manually (but a warning will be printed).
Fixed misuse of getSampleDataRO with lazy data in ripley.
ripley::getSize() also returns circumcircle now.
Fixed integration with constant data on ripley.
Had to remove 'nowait' directive from PDE boundary conditions since all 4 (2D) / 8 (3D) values are updated the way it is implemented. This caused race conditions when updating RHS.
new implementation of FCT solver with some modifications to the python interface
Merged ripley rectangular domain into trunk.
Getting ready to merge with trunk.
Fixed early deallocation of matrix pattern.
Since setX() is not supported the constructors now take location and extent. Python wrapper allows construction the old way with l0,l1 or with tuples (x0,x1),(y0,y1).
3D tests pass, yeah :-)
Finished 2D - all 2D PDE (non-lumped) tests pass now
Adjusted interface for boundary conditions.
BC added
Preparations for PDE boundary element assembly. Removed some code that slipped in before.
weipa processes ripley 3D boundary elements now.
Checkpoint. Getting weipa to like ripley again.
Fixed integration now that elements overlap + code cleanup.
PDE assemblage with reduced order done.
Can solve full order systems of PDEs in 2D and 3D now.
-implemented addPDEToRHS() and setToSize() -added a few missing calls to requireWrite() -added assemblePDESystem() to Rectangle but haven't tested yet
Implemented ownSample for face elements.
Implemented ownSample() for DOF and Elements.
Fixed interpolation from DOF to nodes and moved common code to the base class.
Made element IDs globally unique so ownSample can work (to be implemented). Fixed couple block in 3D so PDEs seem to work now with MPI :-)
checkpoint. 3D single PDEs almost working
Took over Node->DOF interpolation to Brick and parallelised some loops.
can solve PDEs in 2D with MPI now, boundary conditions still missing.
Checkpoint - reinstated DOF FS and tried to fix couple blocks...
PDE assembly in serial and 2D seems to be doing what it's supposed to when boundary elements are not involved.
In Ripley Solution==ContinuousFunction and ReducedSolution==ReducedCF now. Removed a test in escript that relied on failure when trying to tag Data on Solution/ReducedSolution.
Implemented (Face)Elements->Reduced(Face)Elements interpolation and started separating DOF and nodes. Also, implemented operator==() so that a==b for different domain objects a and b which are in the same state.
Set initial node/element/face tags to 0 and updated tests.
escript on ripley tests pass. However, returning 0 for number of point elements and DOF id's are bogus.
Edits to generated code in Rectangle (renamed constants & removed duplicates, more nowait optimizations)
All "util on ripley" tests pass now :-) -added support for node/element/face tags -implemented setToNormal() -added a bunch of omp nowait
Integrals in 2D and 3D for rectangle & brick elements & face elements in regular and reduced orders.
Interpolation & Gradient for reduced elements & faces.
Fixes to the face element code.
Interpolation on faces.
Gradient for Rectangle elements.
Prepared interpolation methods.
getPattern() for Rectangle. Couple node numbering is likely to change once we get to the PDE assembly.
[x] Rectangle... [x] Brick node id's and weipa compatibility
[x] Rectangle node id's and weipa compatibility [ ] Brick...
Next iteration
This form allows you to request diffs between any two revisions of this file. For each of the two "sides" of the diff, enter a numeric revision.
ViewVC Help | |
Powered by ViewVC 1.1.26 |