Gentlemen, Thanks for all your help so far. I have been using fipy Trunk successfully with the development version of numpy (1.3.0dev5945). I get the occasional warning about deprecated names, but everything seems to work fine. I've also gotten PyTrilinos working and the Aztec00 solver with MultiLevelPreconditioner is very fast. I did a quick comparison using a simple Poisson problem between fipy with Aztec00 and a stock Poisson matrix builder solved with pyamg. The fipy code was only about 5% slower, and given all the flexibility in fipy versus roll-your-own solutions, that's really impressive.
The problem I'm having now is that I want to save my simulation results using the h5 file format with pytables. I can extract all the problem related arrays without much difficulty, but the boundary conditions are more vexing since in the forms I can readily access, the locations where each boundary condition is applied live as unresolved union and intersection of boolean arrays buried inside the FixedValue or FixedFlux instances. Is there an approved method for extracting and resolving this location information? If not, is there an easy hack I can use? Best Regards, Angus Hendrick
