Adding CORRECT initial condition, with energy source position and
wavelet defined as a time array.
Well seams that I am loosing the touch... or I've never known that, not working :P . Or python or pylab somehow no plotting working, no convergence.
Propagation just in the first z level and not going down... Movie not working!
Well seams that I am loosing the touch... or I've never known that, not working :P . Or python or pylab somehow no plotting working, no convergence.
Propagation just in the first z level and not going down... Movie not working!
uhuuu its going. Theres something with matplotlib for not being able to plot animation with pyshow, gotta have to fix that, but solution seams to be working.
TODO: Try to find inverse of mUt matrix for other means like LU that will increase performance a lot.
somewhere we are going, now with LU decomposition must fast approach.
Reorganizing the tree, python and testing one side. C code and examples other side
Updating file structure, also 1d finite explicit solution 4order time Lax-Wendroff 2order space working.
4 order in time (Lax-Wendroff), 4 order space
Both taylor finite differences explicit,
TODO 1) add calculation of Dt based on velocity field and convergence criteria
2) better first aproximation of time derivatives at the source location
Besides all that: Perfect working!!
Added functionality to save the solution at every t steps.
Some renames and gif, avi support animation. Working better yet.
1) Need to figure out why velocities are twice bigger than expected
2) Need to figure out how to put sponge boundaries.
implicit implementation is ok, the matrix was being assembled wrongly. figured out witha
a test of grid 10x10.
Merge
Adding pyshape module used as per-before.
To convert cost data from nooa.
Everything solved, pyshape and python code before.
Now everything togueter.
Some modifications in name convetion for the lax method, names chaneged,
included some papers related to the lax method. The method is working but still
need to clearly define de convergence criteria.
Cleaning up the house. Renaming folders an structure.
Done:
+ Imp2DLuWave : + as a code convention example for the rest
                + perfectly working created a cookbook example: Wave2DAcousticSample

+ WaveAnim    : + refactored more coherent and simple
                + added support for 2d wave fields

Todo:
+ finish documentation for implicit 2d wave equation
+ find convergence criteria for 1d lax and document it
+ port segypy or create a library for saving/writing segy files in python
+ code convention for all libraries
+ Sample for cookbook Wave2DAcousticSample created.
+ Modifications on WaveAnim better documented
+ Created base class for 2d fields BaseWave2DField and rewrite Exp2DFourier
using it as base class.
+ Errors with Exp2DFourier exponential growth.
+ Tried to use the segypy library to complicated so far many errors.
+ After started to create my own SegyHeaders.py is part of the effort.
Including the C files original headers from SU for segy/su format.

For now its abandoned since I can plot traces as bitmap using imshow.

+ come back in the future or use other library, a copy of the segypy is
being commited togueter with some might be usefull data (su and segy) some
*.sh too.
Using pydev with eclipse.
Many improvments and change in packages and modules: Filters, CookBook.
Todo: Finish RickerWavelet, BaseWave1DField, Cookbook examples, Seismic