Sparse Engine

[Jean Marshall]

Hi guys:

John and I are starting to write a sparse engine -- for sparse storage of 
our material-dependent fields.  We think we have come up with an 
optimization for our Eulerian code that should really make it scream.  We 
have been studying the IndirectionEngine example, which is very similar to 
what we need.

Unfortunately, the IndirectionEngine example problem, indirect_test1.cpp, 
only demonstrates how to build the engine, not an array or a field.  Could 
someone please show us how to move forward with this example to building a 
complete array and field version?

Unlike the IndirectionEngine example, we only need local communications 
(not all-to-all) along with the same type of guard cell update found in a 
normal field.

Any help will be greatly appreciated!

Jean, John, Don

Jean and Don Marshall

84250 Indio Springs Dr #291
Indio, CA 92203-3413

760-775-1576  home
760-574-0182  Jean's cell
760-574-0192  Don's cell