QuoteOriginally posted by: outrunQuoteOriginally posted by: CuchulainnImo this is a clear candidate for a reusable library.Indeed! I see two separate (orthogonal) concepts / libraries* generic numerical derivatives* mesh / grid coordinates mappings / transforms. This will be some mapping between S,t "space" and mesh/grids space i,j (and in doing that also allow for shadow points, 0 or 1 based addressing)The last one has two sub elements: grid, and coordinate transform. For FDM you'll probably have cube-shaped grids, for trees (binomial, trinonimal) you'll have different (triangular) shapes. Both can have the same coordinate transforms. So we have 3 sub libraries:1. numerical derivatives ala dfdx(x1,x2,x3,f1,f2,f3)2. grids, meshes, trees. Basically integer coordinate regions in N dimensions e.g. things like 0<=i<N 3. coordinate tranforms x,y,z <-> f(i,j,k) what do you think?Bang on. All we need now is someone to do it BTW 1,2 and 3 are closely related to Alan/Financial Alex post-processing error analysis requirements, as they will use 1, 2 and 3.
Last edited by Cuchulainn
on November 9th, 2011, 11:00 pm, edited 1 time in total.