Computing a theoretical dispersion curve - Revision history

Admin: /* Theoretical considerations */

2013-01-11T14:41:42Z

Theoretical considerations

Admin: /* Programming interface */

2010-02-22T18:38:01Z

Programming interface

Admin: /* Programming interface */

2010-02-22T18:29:39Z

Programming interface

Admin: /* Computing the dispersion curves */

2010-02-22T18:27:00Z

Computing the dispersion curves

Admin: /* Saving the dispersion curves */

2010-02-22T18:26:34Z

Saving the dispersion curves

Admin: /* Plotting the dispersion curves */

2010-02-22T18:26:09Z

Plotting the dispersion curves

Admin: /* Saving the dispersion curves */

2010-02-22T18:25:26Z

Saving the dispersion curves

Admin: /* Plotting the dispersion curves */

2010-02-22T18:25:16Z

Plotting the dispersion curves

Admin: /* Creating a model file */

2010-02-22T18:23:16Z

Creating a model file

Marc: /* Plotting the dispersion curves */

2009-07-24T12:21:33Z

Plotting the dispersion curves

← Older revision		Revision as of 14:41, 11 January 2013
Line 183:		Line 183:
	=== Theoretical considerations ===		=== Theoretical considerations ===

−	The equations and methods implemented in ''~~libqtbwave~~.so'' are fully described from a theoretical point of view in [http://marc.geopsy.org/publications.html Marc Wathelet's PhD thesis] (chapter 3).	+	The equations and methods implemented in ''libQGpCoreWave.so'' are fully described from a theoretical point of view in [http://marc.geopsy.org/publications.html Marc Wathelet's PhD thesis] (chapter 3).

@@ Line 73: / Line 73: @@
 === Programming interface ===
-The computation of dispersion is implemented in ''libQGpCoreWave.so'' (or equivalent for Windows and Mac). You can access it through various languages: C++, C or Fortran
+The computation of dispersion is implemented in ''libQGpCoreWave.so'' (or equivalent for Windows and Mac). You can access it through various languages: C++, C or Fortran. Complete examples are provided for C and Fortran in the source distribution archive. Download ''gpdc'' to have the simplest archive.
-A C++ example:
+==== A C++ example ====
   #include <QGpCoreWave.h>
@@ Line 119: / Line 119: @@
   g++ [...] -lQGpCoreTools -lQGpCoreWave
-A Fortran example:
+==== Fortran interface ====
-  call dispersion_curve_init
+  call dispersion_curve_init(verbose)
-  call dispersion_curve_rayleigh(nLayers, h, vp, vs, rho, nSamples, omega, nModes, slowness );
+  call dispersion_curve_rayleigh(nLayers, h, vp, vs, rho, nSamples, omega, nModes, slowness, group);
-  call dispersion_curve_love(nLayers, h, vp, vs, rho, nSamples, omega, nModes, slowness );
+  call dispersion_curve_love(nLayers, h, vp, vs, rho, nSamples, omega, nModes, slowness, group);
 {| border="1" cellpadding="5" cellspacing="0"
@@ Line 130: / Line 130: @@
   ! Type
   ! Description
   |-
   | nLayers
@@ Line 166: / Line 170: @@
   | Vector of floats 8 bytes (nSamples*nModes, elements)
   | Output of slowness values
   |}
@@ Line 171: / Line 179: @@
 To compile:
-  gfortran [...] -lqtbtools -lqtbwave
+  gfortran [...] -lQGpCoreTools -lQGpCoreWave
 === Theoretical considerations ===
 The equations and methods implemented in ''libqtbwave.so'' are fully described from a theoretical point of view in [http://marc.geopsy.org/publications.html Marc Wathelet's PhD thesis] (chapter 3).

@@ Line 73: / Line 73: @@
 === Programming interface ===
-The computation of dispersion is implemented in ''libqtbwave.so'' (or equivalent for Windows and Mac). You can access it through various languages: C++, C or Fortran
+The computation of dispersion is implemented in ''libQGpCoreWave.so'' (or equivalent for Windows and Mac). You can access it through various languages: C++, C or Fortran
 A C++ example:
-  #include "qtblayeredmodel.h"
+  #include <QGpCoreWave.h>
   [...]
   // 3-layer model initialization
-  QtbLayeredModel model(3);
+  LayeredModel model(3);
   model.h()[0]=7.5;
   model.h()[1]=25.0;
@@ Line 101: / Line 98: @@
   [...]
   x[49]=...;
-  // Or use QtbCurve to generate linear or log series
+  // Or use Curve to generate linear or log series
-  QtbCurve<QtbPoint1D> c;
+  Curve<Point1D> c;
   c.line( minFreq, 0.0, maxFreq, 0.0 );
-  c.resample( nSamples, minFreq, maxFreq, Qtb::Log | Qtb::Function );
+  c.resample( nSamples, minFreq, maxFreq, Log | Function );
   c.multiply( 2*M_PI ); // convert to angular frequency
   QVector<double> x = c.xVector();
   // Dispersion computation, fundamental mode only
-  QtbRayleigh rayleigh( &model );
+  Rayleigh rayleigh( &model );
-  //QtbLove love( &model );
+  //Love love( &model );
-  QtbDispersion dispersion( 1, &x);
+  Dispersion dispersion( 1, &x);
   dispersion.calculate( modelRayleigh, 0 ); // 0 is to avoid ellipticity computation
   // Access computed values
-  QtbCurve<QtbPoint2D> * dc = dispersion.curve(0);
+  Curve<Point2D> * dc = dispersion.curve(0);
-  // See QtbCurve header for more information
+  // See Curve header for more information
   // Or directly, fundamental mode, ith sample value.
   dispersion.mode(0)[i]
@@ Line 120: / Line 117: @@
 To compile:
-  g++ [...] -lqtbtools -lqtbwave
+  g++ [...] -lQGpCoreTools -lQGpCoreWave
 A Fortran example:

@@ Line 23: / Line 23: @@
 To get the dispersion curve of Rayleigh fundamental mode:
-  gpdc test.model
+  gpdc Test.model
 You get as output:
@@ Line 43: / Line 43: @@
 To compute Rayleigh higher modes (4 in this case) as well as the fundamental mode:
-  gpdc test.model -R 5
+  gpdc Test.model -R 5
 To compute Love higher modes (4 in this case) as well as the fundamental mode:
-  gpdc test.model -R 0 -L 5
+  gpdc Test.model -R 0 -L 5
 To compute Love modes, you must set explicitly the number of Rayleigh to 0, it is one by default.

@@ Line 53: / Line 53: @@
 == Saving the dispersion curves ==
-  gpdc Test.model -R 2 > test_Rayeigh_2modes.disp
+  gpdc Test.model -R 2 > Test_Rayeigh_2modes.disp
 == Plotting the dispersion curves ==

@@ Line 57: / Line 57: @@
 == Plotting the dispersion curves ==
-[[Image:gpdc_dispersion_plot.png|thumb|right|600px|Dispersion curve plot with figue]]
+[[Image:gpdc_dispersion_plot.png|thumb|right|300px|Dispersion curve plot with figue]]
 We use '''figue''' to plot couples of values (X, Y), additionally we provide a '''make-up''' file to get the correct type and labeling of axes. Before running the next command you can download [[Media:Dc.mkup|Dc.mkup]].

@@ Line 59: / Line 59: @@
 [[Image:gpdc_dispersion_plot.png|thumb|right|600px|Dispersion curve plot with figue]]
-We use '''figue''' to plot couples of values (X, Y), additionally we provide a '''make-up''' file to get the correct type and labeling of axes. Before running the next command you can download [[Media:Dc.mkup|dc.mkup]].
+We use '''figue''' to plot couples of values (X, Y), additionally we provide a '''make-up''' file to get the correct type and labeling of axes. Before running the next command you can download [[Media:Dc.mkup|Dc.mkup]].
-  gpdc test.model -R 2 | figue -c -m dc.mkup
+  gpdc Test.model -R 2 | figue -c -m Dc.mkup
 == Going further ==

@@ Line 5: / Line 5: @@
 == Creating a model file ==
-Crete a text file (e.g. "test.model"):
+Crete a text file (e.g. "[[Media:Test.model|test.model]]"):
   # My first model: two layers over a half-space