pyWafo: Implementing 2D sea state simulation functionality (i.e. SpecData2D.sim())

[stustd]

Are there others interested to in the 2D sea-state simulation
functionality and want to help coding the SpecData2D.sim()? I've made
a preliminary start but I hope I can share the effort...

[Karl Sweitzer]

I would be willing to help.  I have code to create 1D random time histories from given PSD, and have experimented with 2D cross power spectrum techniques.  If you are looking for non Gaussian time histories, then I may be able to help... with some additional assistance from some of my associates.

Let me know...

Regards,

Karl Sweitzer

2011/12/15 stustd <stu...@gmail.com>

Are there others interested to in the 2D sea-state simulation
functionality and want to help coding the SpecData2D.sim()? I've made
a preliminary start but I hope I can share the effort...

--
Du mottar denne meldingen fordi du abonnerer på Google-gruppen «wafo».
Hvis du vil legge inn en melding i denne gruppen, kan du sende e-post til wa...@googlegroups.com.
Hvis du vil avslutte abonnementet på denne gruppen, sender du en e-post til wafo+uns...@googlegroups.com.
 Hvis du vil ha flere alternativer, kan du besøke gruppen på http://groups.google.com/group/wafo?hl=no.

[StUstD]

Thanks Karl, nice you want to help. In SpecData1D.sim()  there are two options to create time histories:

1) Approximate (replicates PSD approximately): i.e. the conventional IFFT way

2) Exact (replicates PSD exactly) through circulant embedding (CE) techniques.

With regard to SpecData2D.sim(), there's CE code available for multivariate Gaussian (random field) simulation  (e.g. www.hermir.org) but most of it is either written in Matlab or in R. Saving the rewrite I think straight linking with some Fortran, C/C++ code would be nicer/faster and quicker to implement. What do you think?

-- Henk