[Gfs-devel] GfsOutputSpectra

[Daniel F. <df...@gm...>]

Hello Stephane

I have pushed again my patch for GfsOutputSpectra in my gerris-bubbles
repository. For now it does not work in parallel. Before doing that I
would like to have some small problems fixed:

1) I do not know how to destroy the object (I know it is simple but
what I do gives me segmentation fault...., for now I just do not
destroy the object)

2) I am not sure if the hierarchy I chose is the correct one. I have
derived GfsOutputSpectra2D and GfsOutputSpectra3D from a general class
called GfsOutputSpectra. Users should not be allowed to use directly
GfsOutputSpectra, where the dimension is not set. I do not know very
well how to avoid users to use GfsOutputSpectra.

If you feel it is fine, I will send you the parallel version (there
are also some minor issues we should discuss)

In any case, the current patch works already in serial. You can test
it with the following example:

GModule fft

1 1 GfsSimulation GfsBox GfsGEdge { } {

   Time { iend = 1 dtmax = 1.e-4 }

   Refine 7

   GfsInit {} { U = sin(4.*M_PI*x)*cos(2.*M_PI*y)
                V = -cos(4.*M_PI*x)*sin(2.*M_PI*y)
       }

   GfsOutputSpectra2D { istep = 1 } file-%ld.dat U { x=0 y=0 z=0
Lx=0.99 Ly=0.99 Lz=0 } 4
   GfsOutputTime { istep =1 } stderr
}
GfsBox {
}
1 1 right
1 1 top
~


-- 
www.danielfuster.com