Classes
class	AbstractComplexPointOptions
	Brief description of class still missing. More...
class	AbstractFunction2
	Abstract 2D function. More...
class	AbstractNumericalCache
	Brief description of class still missing. More...
class	AbstractNumericalKey
	Brief description of class still missing. More...
class	AbstractParameters
	Brief description of class still missing. More...
class	AbstractProgress
	Brief description of class still missing. More...
class	AbstractStream
	Brief description of class still missing. More...
class	AllocatedCacheItem
class	Angle
class	App
class	ApplicationClock
	Brief description of class still missing. More...
class	ApplicationHelp
	Brief description of class still missing. More...
class	ArgumentStdinReader
	Brief description of class still missing. More...
class	AsciiLineParser
	A fast low level parser for ASCII contents. More...
class	AtomicBoolean
	A boolean value tha can can accessed in an atomic way. More...
class	Cache
	Cache space with "inteligent" swap mechanism. More...
class	CacheItem
	Brief description of class still missing. More...
class	CacheProcess
	A list of CacheItem to be processed. More...
class	Circle
class	ColumnTextDescription
	Storage for the decription of column in a ColumnTextParser. More...
class	ColumnTextIterator
	Brief description of class still missing. More...
class	ColumnTextParser
	Brief description of class still missing. More...
class	Complex
class	ComplexExponential
	Fast computation of complex exponential with a rational argument. More...
class	ComplexMatrix
	Brief description of class still missing. More...
class	ComplexPointOptions
	Brief description of class still missing. More...
class	ComplexStatisticalPoint
class	ComplexStatisticalValue
class	ConcentricSearch
	Searches for maxima across a grid in concentric way. More...
class	ConsoleProgress
	Brief description of class still missing. More...
class	CoreApplication
	Brief description of class still missing. More...
class	CoreApplicationPrivate
	Brief description of class still missing. More...
class	Covariance
	Brief description of class still missing. More...
class	Curve
	Template class for a curve It is a vector of "points" that may eventually contain extra data. The "point" type must at least have the following functions: More...
class	CurvePointOptions
	Storage for curve point options. More...
class	DateTime
	A high precision time class without limitation (e.g. 1970-01-01) More...
class	DateTimeData
class	DateTimeFromString
class	DateTimeParser
class	DateTimeToString
class	DirectionalSearch
	Search a grid for maxima (global or local) along one direction only. More...
class	DoubleMatrix
	Brief description of class still missing. More...
class	DynamicBuffer
class	Ellipse
	Brief description of class still missing. More...
class	EventRestrictor
	Restricts the number of events per unit of time. More...
class	ExpressionAction
class	ExpressionActions
	Brief description of class still missing. More...
class	ExpressionBracket
class	ExpressionContext
	Brief description of class still missing. More...
class	ExpressionFunction
class	ExpressionOperator
class	ExpressionParser
class	ExpressionStorage
	Brief description of class still missing. More...
class	ExpressionString
	Brief description of class still missing. More...
class	ExpressionValue
class	ExpressionVariable
class	ExpressionVariableStorage
class	File
class	FileStream
	Brief description of class still missing. More...
class	FletcherChecksum
	Brief description of class still missing. More...
class	Function2Search
	Brief description of class still missing. More...
class	Function2SearchMaximum
class	GaussDistribution
class	GlobalObject
class	GoogleEarthKML
	Brief description of class still missing. More...
class	Grid2D
	A 2D grid with regular cell sizes. More...
class	Grid3D
	A 3D grid with regular cell sizes. More...
class	GridSearch
	Search a grid for maxima (global or local) More...
class	HashNumericalKey
class	Histogram
	Brief description of class still missing. More...
class	IncreaseStorage
	Brief description of class still missing. More...
class	IrregularGrid2D
	A 2D grid where the coordinates of the rectangular cell are not regularly distributed. More...
class	IrregularGrid2DData
	Brief description of class still missing. More...
class	Leds
	Brief description of class still missing. More...
class	Line2D
class	LineParser
	Brief description of class still missing. More...
class	LocalPointer
	Automatically deletes a pointer upon exit of a function. More...
class	LoopTask
	Task of a parallel loop. More...
class	Matrix
class	Matrix2x2
class	Matrix3x3
	Brief description of class still missing. More...
class	Matrix4x4
class	MatrixColumnIterator
class	MatrixData
	Brief description of class still missing. More...
class	MatrixIterator
	Brief description of class still missing. More...
class	MatrixMultiply
	Brief description of class still missing. More...
class	MatrixMultiplyKey
class	MatrixRowIterator
class	Message
	Show messages to user If not graphical user interface is running this class send messages to and get information from the user through the console. In graphical user interface, a GuiMessage class is also available (in library qtbguitools), which transparently (for the developer) send the message to dialog boxes. Those two classes have eaxctly the same API. So always use this class for send messages and automatically the messages will be displayed in the appropriate way. More...
class	MessageContext
	Define a context for recording and automatic answer to user messages. More...
class	MutableMatrixColumnIterator
class	MutableMatrixIterator
class	MutableMatrixRowIterator
class	NamedPoint
	Brief description of class still missing. More...
class	Number
class	PackageInfo
	Brief description of class still missing. More...
class	Parallelepiped
	Brief description of class still missing. More...
class	ParallelLoop
	Class that helps to parallelize loops. More...
class	PathTranslator
	Utility to ease access to file which path is stored as absolute. It can translate obsolete paths into new ones, for instance, if the user change the organisation of his file tree. The new path is only asked one time. For all other requests, the obsolete path is searched in a map and if found, the file path is modified. More...
class	PermutationVector
	To perfom permutations in a matrix. More...
struct	PitsaTime
class	Plane
	A geometrical plane. More...
class	PluginCoreApplication
	Brief description of class still missing. More...
class	Point
class	Point1D
	Brief description of class still missing. More...
class	Point2D
class	PointLocate
	Brief description of class still missing. More...
class	PointND
	Brief description of class still missing. More...
class	PointNDOptions
	Brief description of class still missing. More...
class	ProcessStatistics
class	Random
	Random generator for intensive random computations. More...
class	RealStatisticalPoint
class	RealStatisticalValue
class	RealValue
class	Rect
	Storage for a rectangle with coordinates in double precision A rectangle is defined by top-left corner (x1,y1) and bottom-right corner (x2,y2). All function that changes the coordinates of these two points check for validity. More...
class	RootSolver
class	RootSolverTemplate
	Template class to search roots of a 1D function. More...
class	SamplingParameters
	An axis sampling definition. More...
class	Segment
	Brief description of class still missing. More...
class	Segment2D
	A 2D segment made of two points. More...
class	SharedObject
	A simple thread-safe shared object. More...
class	SmoothingParameters
	Brief description of class still missing. More...
class	SortedVector
	A vector of sorted elements. More...
class	SparseMatrix
	Sparse matrix storage based on rows. More...
class	StandardStream
	A stream directed to stdout or stderr. More...
class	StatComparator
class	Stationarity
class	StatisticalPoint
	Brief description of class still missing. More...
class	StatisticalValue
	Brief description of class still missing. More...
class	Statistics
	Brief description of class still missing. More...
class	StreamRedirection
	Begin and end stream redirection. More...
class	StringSection
	Section of a QString. More...
class	StringStream
	A stream directed to a string. More...
class	Tar
	Brief description of class still missing. More...
struct	TarHeader
class	Thread
	Thread with an internal stack to record trace points. More...
class	ThreadTimer
class	TraceBug
	Brief description of class still missing. More...
class	TraceInfo
	Abstract class for all types of trace information. More...
class	TraceLog
	Brief description of class still missing. More...
class	TraceStamp
	Storage for trace stamps. More...
class	TreeContainer
	Brief description of class still missing. More...
class	TreeItem
	Brief description of class still missing. More...
class	Value
	Storage for any double precision value with a flag for validitidy. More...
class	Version
	Brief description of class still missing. More...
class	XMLByteArrayStream
	Brief description of class still missing. More...
class	XMLClass
	XMLClass is a parser for a simplified XML syntax. More...
class	XMLClassCreator
class	XMLClassFactory
	XMLClassFactory provides creation of instances of classes from their XML tag. More...
class	XMLContext
	Brief description of class still missing. More...
class	XMLErrorReport
	Brief description of class still missing. More...
class	XMLGenericItem
	Brief description of class still missing. More...
class	XMLHeader
	Brief description of class still missing. More...
class	XMLMap
	Save/Restore a QMap to a XML stream. More...
class	XMLMember
	XMLMember is used by XMLClass to receive the type of tag to parse. More...
class	XMLParser
	Brief description of class still missing. More...
class	XMLSaveAttribute
class	XMLSaveAttributes
	Brief description of class still missing. More...
class	XMLStream
	Brief description of class still missing. More...
class	XMLStringStream
	Brief description of class still missing. More...
class	XMLTarStream
	Brief description of class still missing. More...
Typedefs
typedef double v2df	__attribute__ ((vector_size(16)))
typedef AbstractStream &(*	AbstractStreamFunction )(AbstractStream &)
typedef Value< Complex >	ComplexValue
typedef QList < Function2SearchMaximum > ::iterator	Function2MaximaIterator
typedef struct QGpCoreTools::PitsaTime	TIME
typedef QHash< StringSection, StringSection > ::const_iterator	XMLRestoreAttributeIterator
typedef QHash< StringSection, StringSection >	XMLRestoreAttributes
Enumerations
enum	AxisType { XAxis, YAxis, ZAxis }
enum	MisfitType { L1, L1_Normalized, L1_LogNormalized, L1_NormalizedBySigmaOnly, L2, L2_Normalized, L2_LogNormalized, L2_NormalizedBySigmaOnly, Akaike, AkaikeFewSamples }
enum	SamplingOption { LinearScale = 0, LogScale = 1, InversedScale = 2, Interpole = 4, Function = 8 }
Functions
double	abs (double a)
double	abs (const Complex &c)
double	abs2 (double a)
double	abs2 (const Complex &c)
void	amoeba (double *p, double y[], int ndim, double ftol, double(funk)(double[]), int *nfunk)
double	amotry (double *p, double y[], double psum[], int ndim, double(funk)(double[]), int ihi, double fac)
Complex	asin (const Complex &c)
Complex	conjugate (const Complex &c)
Complex	cos (const Complex &c)
int	dgesvd_ (char jobu, char jobvt, int m, int n, double a, int lda, double s, double u, int ldu, double vt, int ldvt, double work, int lwork, int info)
int	dgetrf_ (int m, int n, double a, int lda, int ipiv, int info)
int	dgetri_ (int n, double a, int lda, int ipiv, double work, int lwork, int *info)
QString	encodeToHtml (QString str)
AbstractStream &	endl (AbstractStream &s)
Complex	exp (const Complex &c)
AbstractStream &	flush (AbstractStream &s)
Complex	inverse (const Complex &c)
Complex	log (const Complex &c)
Complex	operator* (const double &d, const Complex &c)
Complex	operator+ (const double &d, const Complex &c)
Complex	operator- (const double &d, const Complex &c)
Complex	operator/ (const double &d, const Complex &c)
QTextStream &	operator<< (QTextStream &s, const ComplexMatrix &m)
AbstractStream &	operator<< (AbstractStream &s, const DoubleMatrix &m)
QTextStream &	operator<< (QTextStream &s, const Complex &c)
template<class numberType >
QGPCORETOOLS_EXPORT QDataStream &	operator<< (QDataStream &s, const StatisticalValue< numberType > &p)
AbstractStream &	operator<< (AbstractStream &s, AbstractStreamFunction f)
QDataStream &	operator<< (QDataStream &s, const Complex &c)
template<class ValueType >
QGPCORETOOLS_EXPORT QTextStream &	operator<< (QTextStream &s, Grid2D< ValueType > &grd)
template<class ValueType >
QGPCORETOOLS_EXPORT QTextStream &	operator<< (QTextStream &s, Grid3D< ValueType > &grd)
QTextStream &	operator<< (QTextStream &s, const IrregularGrid2D &grd)
QTextStream &	operator<< (QTextStream &s, const QList< Point > &plist)
QTextStream &	operator<< (QTextStream &s, const Point &p)
QDataStream &	operator<< (QDataStream &s, const Point &p)
template<typename T >
QDataStream &	operator<< (QDataStream &s, const Matrix< T > &m)
template<typename T >
QDataStream &	operator<< (QDataStream &s, const SparseMatrix< T > &m)
template<class ValueType >
QTextStream &	operator<< (QTextStream &s, const Grid2D< double > &grd)
template<class ValueType >
QTextStream &	operator<< (QTextStream &s, const Grid3D< double > &grd)
QTextStream &	operator<< (QTextStream &s, const Point2D &p)
QDataStream &	operator<< (QDataStream &s, const Point2D &p)
QTextStream &	operator<< (QTextStream &s, const QList< Point2D > &plist)
QTextStream &	operator>> (QTextStream &s, IrregularGrid2D &grd)
QTextStream &	operator>> (QTextStream &s, Complex &c)
template<class numberType >
QGPCORETOOLS_EXPORT QDataStream &	operator>> (QDataStream &s, StatisticalValue< numberType > &p)
QDataStream &	operator>> (QDataStream &s, Complex &c)
template<class ValueType >
QGPCORETOOLS_EXPORT QTextStream &	operator>> (QTextStream &s, Grid2D< ValueType > &grd)
template<class ValueType >
QGPCORETOOLS_EXPORT QTextStream &	operator>> (QTextStream &s, Grid3D< ValueType > &grd)
QTextStream &	operator>> (QTextStream &s, QList< Point > &plist)
QTextStream &	operator>> (QTextStream &s, Point &p)
QDataStream &	operator>> (QDataStream &s, Point &p)
template<typename T >
QDataStream &	operator>> (QDataStream &s, Matrix< T > &m)
template<typename T >
QDataStream &	operator>> (QDataStream &s, SparseMatrix< T > &m)
QTextStream &	operator>> (QTextStream &s, Point2D &p)
QDataStream &	operator>> (QDataStream &s, Point2D &p)
QTextStream &	operator>> (QTextStream &s, QList< Point2D > &plist)
	PACKAGE_INFO (QGpCoreTools, QGPCORETOOLS)
uint	qHash (StringSection key)
int	qHash (const HashNumericalKey &key)
uint	qHash (Point2D key)
Complex	sin (const Complex &c)
Complex	sqrt (const Complex &c)
Complex	tan (const Complex &c)
QString	tr (const char s, const char comments=0)
template<typename T >
void	unique (QList< T > &list)
template<typename T , typename Equal >
void	unique (QList< T > &list, Equal equal)
template<typename T >
void	unique (QVector< T > &v)
template<typename T , typename Equal >
void	unique (QVector< T > &v, Equal equal)
template<typename T >
void	uniqueDeleteAll (QList< T > &list)
template<typename T , typename Equal >
void	uniqueDeleteAll (QList< T > &list, Equal equal)
template<typename T >
void	uniqueDeleteAll (QVector< T > &v)
template<typename T , typename Equal >
void	uniqueDeleteAll (QVector< T > &v, Equal equal)

Detailed Description

QGpCoreTools.h

Typedef Documentation

typedef double v2df QGpCoreTools::__attribute__((vector_size(16)))

typedef AbstractStream&(* QGpCoreTools::AbstractStreamFunction)(AbstractStream &)

typedef Value<Complex> QGpCoreTools::ComplexValue

typedef QList<Function2SearchMaximum>::iterator QGpCoreTools::Function2MaximaIterator

typedef struct QGpCoreTools::PitsaTime QGpCoreTools::TIME

typedef QHash<StringSection, StringSection>::const_iterator QGpCoreTools::XMLRestoreAttributeIterator

typedef QHash<StringSection, StringSection> QGpCoreTools::XMLRestoreAttributes

Enumeration Type Documentation

enum QGpCoreTools::AxisType

Enumerator:

XAxis
YAxis
ZAxis

{XAxis, YAxis, ZAxis};

enum QGpCoreTools::MisfitType

Enumerator:

L1
L1_Normalized
L1_LogNormalized
L1_NormalizedBySigmaOnly
L2
L2_Normalized
L2_LogNormalized
L2_NormalizedBySigmaOnly
Akaike
AkaikeFewSamples

                {L1, L1_Normalized, L1_LogNormalized, L1_NormalizedBySigmaOnly,
                 L2, L2_Normalized, L2_LogNormalized, L2_NormalizedBySigmaOnly,
                 Akaike, AkaikeFewSamples};

enum QGpCoreTools::SamplingOption

Enumerator:

LinearScale
LogScale
InversedScale
Interpole
Function

{LinearScale=0, LogScale=1, InversedScale=2, Interpole=4, Function=8};

Function Documentation

double QGpCoreTools::abs ( double a ) [inline]

{
  return ::fabs(a);
}

double QGpCoreTools::abs ( const Complex & c ) [inline]

References QGpCoreTools::Complex::abs().

Referenced by PPSBuffer::bytesAvailable(), GeopsyCore::Gse::checksum(), TFAResults::compute(), SciFigs::ImageLayer::grayFilter(), SciFigs::MouseTracking::isEmptyRectangle(), Station::lessThan(), Event::lessThan(), QGpCoreTools::StatisticalValue< numberType >::misfit(), SciFigs::GraphicObject::mouseMoveEvent(), SciFigs::CircleMask::paintData(), SciFigs::CircleViewer::paintData(), GeopsyCore::FilterParameters::setCausal(), GeopsyCore::FilterParameters::setOrder(), QGpCoreTools::Complex::setPhase(), and DinverCore::Neighborhood::setVolumes().

{
  return c.abs();
}

double QGpCoreTools::abs2 ( double a ) [inline]

{
  return a*a;
}

double QGpCoreTools::abs2 ( const Complex & c ) [inline]

References QGpCoreTools::Complex::abs2().

Referenced by QGpCoreTools::Complex::abs(), and QGpCoreTools::StatisticalValue< numberType >::average().

{
  return c.abs2();
}

QGPCORETOOLS_EXPORT void QGpCoreTools::amoeba	(	double **	p,
		double	y[],
		int	ndim,
		double	ftol,
		double(*)(double[])	funk,
		int *	nfunk
	)

Simplex minimization

Multidimensional minimization of the function funk(x) where x[1..ndim] is a vector in ndim dimensions, by the downhill simplex method of Nelder and Mead. The matrix p[1..ndim+1] [1..ndim] is input. Its ndim+1 rows are ndim-dimensional vectors which are the vertices of the starting simplex. Also input is the vector y[1..ndim+1], whose components must be preinitialized to the values of funk evaluated at the ndim+1 vertices (rows) of p; and ftol the fractional convergence tolerance to be achieved in the function value (n.b.!). On output, p and y will have been reset to ndim+1 new points all within ftol of a minimum function value, and nfunk gives the number of function evaluations taken.

References amotry(), NMAX, SWAP, TINY, and TRACE.

Referenced by DampingResults::compute().

{
  TRACE;
  double amotry(double **p, double y[], double psum[], int ndim,
                 double ( *funk) (double [] ), int ihi, double fac);
  int i, ihi, ilo, inhi, j, mpts=ndim + 1;
  double rtol, sum, swap, ysave, ytry, *psum;

  psum=new double [ ndim ];
  *nfunk=0;
  // GET_PSUM
  for(j=0;j < ndim;j++ )
  {
    for(sum=0.0, i=0;i < mpts;i++ )
      sum += p[ i ][ j ];
    psum[ j ]=sum;
  }

  for( ;; )
  {
    ilo=0;
    /* First we must determine which point is the highest (worst),
       next-highest, and lowest (best), by looping over the points
       in the simplex. */
    ihi=y[ 0 ] > y[ 1 ] ? (inhi=1, 0) : (inhi=0, 1);
    for(i=0;i < mpts;i++ ) {
      if(y[ i ] <= y[ ilo ] )
        ilo=i;
      if(y[ i ] > y[ ihi ] ) {
        inhi=ihi;
        ihi=i;
      } else if(y[ i ] > y[ inhi ] && i!=ihi)
        inhi=i;
    }
    rtol=2.0 * fabs(y[ ihi ] - y[ ilo ] )/(fabs(y[ ihi ] ) + fabs(y[ ilo ] ) + TINY);
    /* Compute the fractional range from highest to lowest and return if
       satisfactory. */
    if(rtol < ftol) {
      // If returning, put best point and value in slot 1.
      SWAP(y[ 0 ], y[ ilo ] );
      for(i=0;i < ndim;i++ )
        SWAP(p[ 0 ][ i ], p[ ilo ][ i ] );
      break;
    }
    if( *nfunk >= NMAX) {
      fprintf(stderr, "NMAX exceeded\n" );
      return ;
    }
    *nfunk += 2;
    /* Begin a new iteration. First extrapolate by a factor  1 through the
       face of the simplex across from the high point, i.e., reflect the
       simplex from the high point. */
    ytry=amotry(p, y, psum, ndim, funk, ihi, -1.0);
    if(ytry <= y[ ilo ] )
      /* Gives a result better than the best point, so try an additional
         extrapolation by a factor 2. */
      ytry=amotry(p, y, psum, ndim, funk, ihi, 2.0);
    else if(ytry >= y[ inhi ] ) {
      /* The reflected point is worse than the second-highest, so look for an
         intermediate lower point, i.e., do a one-dimensional contraction. */
      ysave=y[ ihi ];
      ytry=amotry(p, y, psum, ndim, funk, ihi, 0.5);
      if(ytry >= ysave) {
        /* Can t seem to get rid of that high point. Better contract around the
           lowest (best) point. */
        for(i=0;i < mpts;i++ ) {
          if(i!=ilo) {
            for(j=0;j < ndim;j++ )
              p[ i ][ j ]=psum[ j ]=0.5 * (p[ i ][ j ] + p[ ilo ][ j ] );
            y[ i ]=( *funk) (psum);
          }
        }
        *nfunk += ndim;
        // Keep track of function evaluations. GET_PSUM Recompute psum.
      }
    } else
      --( *nfunk); // Correct the evaluation count
  }  // Go back for the test of doneness and the next iteration.
  delete [] psum;
}

double QGpCoreTools::amotry	(	double **	p,
		double	y[],
		double	psum[],
		int	ndim,
		double(*)(double[])	funk,
		int	ihi,
		double	fac
	)

References TRACE.

Referenced by amoeba().

{
  TRACE;
  int j;
  double fac1, fac2, ytry, *ptry;
  ptry=new double [ ndim ];
  fac1=(1.0 - fac)/ndim;
  fac2=fac1 - fac;
  for(j=0;j < ndim;j++ )
    ptry[ j ]=psum[ j ] * fac1 - p[ ihi ][ j ] * fac2;
  ytry=( *funk) (ptry);
  // Evaluate the function at the trial point.
  if(ytry < y[ ihi ] )
  {
    // If it s better than the highest, then replace the highest.
    y[ ihi ]=ytry;
    for(j=0;j < ndim;j++ ) {
      psum[ j ] += ptry[ j ] - p[ ihi ][ j ];
      p[ ihi ][ j ]=ptry[ j ];
    }
  }
  delete [] ptry;
  return ytry;
}

Complex QGpCoreTools::asin ( const Complex & c ) [inline]

asin(C)=-i log (sqrt(1-C^2)+i*C), where C is this

References log(), and sqrt().

Referenced by QGpCoreWave::RefractionDippingModel::begin(), and QGpCoreTools::Point2D::geographicalDistanceTo().

{
  return Complex(0.0,-1.0)*log(sqrt(Complex(1.0)-c*c)+Complex(0.0,1.0)*c);
}

Complex QGpCoreTools::conjugate ( const Complex & c ) [inline]

References QGpCoreTools::Complex::im, and QGpCoreTools::Complex::re.

Referenced by QGpCoreWave::MagnetoTelluricFactory::calculate(), GeopsyCore::DoubleSignal::complexSpectrum(), ArrayCore::HRFK::crossCorrelationMatrix(), PhaseShifter::initGrid(), ArrayCore::HRFK::initOperator(), QGpCoreTools::Complex::operator/(), operator/(), QGpCoreTools::Complex::operator/=(), ArrayCore::HRFK::value(), and PhaseShifter::value().

{
  return Complex(c.re(),-c.im());
}

Complex QGpCoreTools::cos ( const Complex & c ) [inline]

References QGpCoreTools::Complex::im, QGpCoreTools::Complex::re, and sin().

Referenced by GeopsyCore::TaperDelegate::apply(), QGpCoreWave::RefractionDippingModel::begin(), GeopsyCore::DoubleSignal::discreteFourierTransform(), Simulator::displaySource(), EllipseItem::EllipseItem(), SignalReader::exec(), exp(), QGpCoreTools::ExpressionContext::functionValue(), QGpCoreTools::Point2D::geographicalAzimuthTo(), QGpCoreTools::Point2D::geographicalDistanceTo(), QGpCoreTools::Point2D::geographicalToRectangular(), QGpCoreTools::Point2D::geographicalToUtm(), HVStationSignals::horizontal(), DrawScene::moveTo(), SignalReader::parse(), QGpCoreTools::PointLocate::position(), QGpCoreWave::RefractionDippingModel::propagate(), QGpCoreTools::Point2D::rectangularToGeographical(), QGpCoreTools::DirectionalSearch::setAzimuth(), SourceParameters::setAzimuth(), SourceParameters::setDistance(), QGpCoreTools::Complex::setExp(), FKArrayMap::setSourceSignals(), QGpCoreTools::Complex::setUnitExp(), GeopsyCore::DoubleSignal::shift(), sin(), sqrt(), tan(), QGpCoreTools::Point2D::utmToGeographical(), QGpCoreWave::TheoreticalFK::value(), QGpCoreTools::Ellipse::xRadius(), QGpCoreWave::Love::y(), QGpCoreWave::RayleighTemplate< RealType >::y(), and QGpCoreTools::Ellipse::yRadius().

{
  return Complex( ::cos(c.re()) *::cosh(c.im()), -::sin(c.re()) * sinh(c.im()));
}

int QGpCoreTools::dgesvd_	(	char *	jobu,
		char *	jobvt,
		int *	m,
		int *	n,
		double *	a,
		int *	lda,
		double *	s,
		double *	u,
		int *	ldu,
		double *	vt,
		int *	ldvt,
		double *	work,
		int *	lwork,
		int *	info
	)

Referenced by QGpCoreTools::DoubleMatrix::eigenVectors().

int QGpCoreTools::dgetrf_	(	int *	m,
		int *	n,
		double *	a,
		int *	lda,
		int *	ipiv,
		int *	info
	)

Referenced by QGpCoreTools::DoubleMatrix::invert().

int QGpCoreTools::dgetri_	(	int *	n,
		double *	a,
		int *	lda,
		int *	ipiv,
		double *	work,
		int *	lwork,
		int *	info
	)

Referenced by QGpCoreTools::DoubleMatrix::invert().

QString QGpCoreTools::encodeToHtml ( QString str )

References str.

{
  str.replace("&","&amp;");
  str.replace("'","&apos;");
  str.replace("\"","&quot;");
  str.replace("<","&lt;");
  str.replace(">","&gt;");
  return str;
}

QGPCORETOOLS_EXPORT AbstractStream & QGpCoreTools::endl ( AbstractStream & s )

Not a flush like in QTextStream. This is the mandatory item to terminate a line.

References QGpCoreTools::AbstractStream::endl().

  {
    s.endl();
    return s;
  }

Complex QGpCoreTools::exp ( const Complex & c ) [inline]

References cos(), QGpCoreTools::Complex::im, QGpCoreTools::Complex::re, and sin().

Referenced by GeopsyCore::TaperDelegate::apply(), GeopsyCore::AsciiSignalFormat::assign(), QGpCoreTools::StatComparator::betai(), GaussianFrequencyBand::calculate(), QGpCoreWave::MagnetoTelluricFactory::calculate(), DampingResults::compute(), QGpCompatibility::CompatFunction::getModes(), GeopsyCore::SEGYTraceHeader::IBMFloat2ieee(), DinverCore::ParameterGrid::init(), DampingResults::misfitFunk(), GeopsyCore::DoubleSignal::morletWavelet(), DinverCore::PdfPoint::probability(), QGpCoreTools::GaussDistribution::rho(), GroupView::select(), ArrayCore::FK::setFrequencyBand(), Event::setId(), FKArrayMap::setSourceSignals(), QGpCoreTools::IrregularGrid2DData::smoothInvalidCells(), QGpCoreTools::Number::toDouble(), QGpCoreWave::Love::y(), and QGpCoreWave::RayleighTemplate< RealType >::y().

{
  double r=::exp(c.re());
  return Complex(r*::cos(c.im()), r*::sin(c.im()));
}

QGPCORETOOLS_EXPORT AbstractStream & QGpCoreTools::flush ( AbstractStream & s )

References QGpCoreTools::AbstractStream::flush().

Referenced by GeopsyCore::GeopsyCoreEngine::askLoadFormat(), GeopsyCore::GeopsyCoreEngine::askTimeReference(), QGpCoreTools::ApplicationHelp::exec(), SignalReader::exec(), QGpCoreWave::Seismic1DModel::exportHerrmann(), CurveReader::generateFunction(), QGpCoreTools::Message::getExistingDirectoryInternal(), QGpCoreTools::Message::getOpenFileNameInternal(), QGpCoreTools::Message::getOpenFileNamesInternal(), QGpCoreTools::Message::getSaveFileNameInternal(), DinverDCCore::TargetList2D::humanInfo(), DinverDCCore::Param3DGroundModel::humanInfo(), main(), QGpCoreTools::Message::message(), operator<<(), QGpCoreTools::CoreApplicationPrivate::printArgumentLists(), replace(), DinverDCCore::ModalStorageReader::toStream(), QGpCoreWave::MagnetoTelluricFactory::toStream(), QGpCoreWave::Profile::toStream(), QGpCoreWave::RefractionFactory::toStream(), QGpCoreWave::RefractionDippingModel::toStream(), and QGpCoreWave::Seismic1DModel::toStream().

  {
    s.flush();
    return s;
  }

Complex QGpCoreTools::inverse ( const Complex & c ) [inline]

References QGpCoreTools::Complex::abs2(), QGpCoreTools::Complex::im, and QGpCoreTools::Complex::re.

{
  double m2=1.0/c.abs2();
  return Complex(c.re() * m2, -m2 * c.im());
}

Complex QGpCoreTools::log ( const Complex & c ) [inline]

References QGpCoreTools::Complex::abs(), and QGpCoreTools::Complex::phase().

{
  return Complex(::log(c.abs()), c.phase());
}

Complex QGpCoreTools::operator*	(	const double &	d,
		const Complex &	c
	)		`[inline]`

References QGpCoreTools::Complex::im, and QGpCoreTools::Complex::re.

Referenced by QGpCoreTools::Matrix2x2::operator*(), QGpCoreTools::Matrix3x3::operator*(), and QGpCoreTools::Matrix4x4::operator*().

{
  return Complex(d*c.re(), d*c.im());
}

Complex QGpCoreTools::operator+	(	const double &	d,
		const Complex &	c
	)		`[inline]`

{
  return c+d;
}

Complex QGpCoreTools::operator-	(	const double &	d,
		const Complex &	c
	)		`[inline]`

References QGpCoreTools::Complex::im, and QGpCoreTools::Complex::re.

{
  return Complex(d-c.re(), -c.im());
}

Complex QGpCoreTools::operator/	(	const double &	d,
		const Complex &	c
	)		`[inline]`

References QGpCoreTools::Complex::abs2(), and conjugate().

{
  double m2=c.abs2();
  Complex r(d);
  r*=conjugate(c);
  r/=m2;
  return r;
}

QGPCORETOOLS_EXPORT QTextStream & QGpCoreTools::operator<<	(	QTextStream &	s,
		const ComplexMatrix &	m
	)

References QGpCoreTools::Matrix< T >::columnCount(), flush(), QGpCoreTools::MatrixIterator< T >::hasNext(), QGpCoreTools::MatrixIterator< T >::next(), QGpCoreTools::Complex::re, QGpCoreTools::Matrix< T >::rowCount(), and TRACE.

{
  TRACE;
  s << "\t";
  for(int j=0;j < m.columnCount();j++ ) s << j << "\t";
  s << "\n";
  for(int i=0;i < m.rowCount();i++ ) {
    MatrixRowIterator<Complex> it(m, i);
    s << i << ":\t";
    while(it.hasNext()) {
      it.next();
      s << (*it).re() << "+i*" << (*it).im() << "\t";
    }
    s << "\n";
  }
  s << flush;
  return s;
}

QGPCORETOOLS_EXPORT AbstractStream & QGpCoreTools::operator<<	(	AbstractStream &	s,
		const DoubleMatrix &	m
	)

References QGpCoreTools::Matrix< T >::columnCount(), QGpCoreTools::MatrixIterator< T >::hasNext(), QGpCoreTools::MatrixIterator< T >::next(), QGpCoreTools::Matrix< T >::rowCount(), and TRACE.

  {
    TRACE;
    s << "\t";
    for(int j=0;j < m.columnCount();j++) s << QString::number(j) << "\t";
    s << "\n";
    for(int i=0;i < m.rowCount();i++ ) {
      MatrixRowIterator<double> it(m, i);
      s << QString::number(i) << ":\t";
      while(it.hasNext()) {
        it.next();
        s << QString::number(*it) << "\t";
      }
      s << "\n";
    }
    return s;
  }

QGPCORETOOLS_EXPORT QTextStream & QGpCoreTools::operator<<	(	QTextStream &	s,
		const Complex &	c
	)

References QGpCoreTools::Complex::im, QGpCoreTools::Complex::re, and TRACE.

{
  TRACE;
  s.setRealNumberNotation(QTextStream::SmartNotation);
  s.setRealNumberPrecision(20);
  s << c.re() << " " << c.im();
  return s;
}

template<class numberType >

QDataStream & QGpCoreTools::operator<<	(	QDataStream &	s,
		const StatisticalValue< numberType > &	p
	)

References TRACE.

{
  TRACE;
  s << p.mean() << p.stddev() << p.weight();
  return s;
}

AbstractStream& QGpCoreTools::operator<<	(	AbstractStream &	s,
		AbstractStreamFunction	f
	)		`[inline]`

  {
    return (*f)(s);
  }

QGPCORETOOLS_EXPORT QDataStream & QGpCoreTools::operator<<	(	QDataStream &	s,
		const Complex &	c
	)

References QGpCoreTools::Complex::im, QGpCoreTools::Complex::re, and TRACE.

{
  TRACE;
  s << c.re() << c.im();
  return s;
}

template<class ValueType >

QGPCORETOOLS_EXPORT QTextStream& QGpCoreTools::operator<<	(	QTextStream &	s,
		Grid2D< ValueType > &	grd
	)

template<class ValueType >

QGPCORETOOLS_EXPORT QTextStream& QGpCoreTools::operator<<	(	QTextStream &	s,
		Grid3D< ValueType > &	grd
	)

QGPCORETOOLS_EXPORT QTextStream & QGpCoreTools::operator<<	(	QTextStream &	s,
		const IrregularGrid2D &	grd
	)

References endl(), QGpCoreTools::IrregularGrid2D::nx(), QGpCoreTools::IrregularGrid2D::ny(), TRACE, QGpCoreTools::IrregularGrid2D::valuePointer(), QGpCoreTools::IrregularGrid2D::x(), and QGpCoreTools::IrregularGrid2D::y().

{
  TRACE;
  s << "x y val" << ::endl;
  const double * ptr=grd.valuePointer(0, 0);
  for(int iy=0;iy < grd.ny();iy++ ) {
    for(int ix=0;ix < grd.nx();ix++ )
      s << grd.x(ix) << " " << grd.y(iy) << " " << *(ptr++ ) << ::endl;
  }
#if 0
  // Old format for grid, not practical for processing with matlab or gnuplot for instance
  s << "nx=" << grd.nx() << endl;
  s << "ny=" << grd.ny() << endl;
  s << "x" << endl;
  int i;
  for(i=0;i < grd.nx();i++ ) s << grd.x(i) << endl;
  s << "y" << endl;
  for(i=0;i < grd.ny();i++ ) s << grd.y(i) << endl;
  s << "values" << endl;
  double * ptr=grd.valueAddress(0, 0);
  for(i=0;i < grd.ny();i++ ) {
    for(int ix=0;ix < grd.nx();ix++ )
      s << *(ptr++ ) << " " << ix << " " << i << endl;
  }
#endif
  return s;
}

QGPCORETOOLS_EXPORT QTextStream & QGpCoreTools::operator<<	(	QTextStream &	s,
		const QList< Point > &	plist
	)

References endl(), and TRACE.

  {
    TRACE;
    s << "%x      y      z"<< ::endl;
    QList<Point>::const_iterator it;
    for(it=plist.begin();it!=plist.end();++it) s << *it << ::endl;
    return s;
  }

QGPCORETOOLS_EXPORT QTextStream & QGpCoreTools::operator<<	(	QTextStream &	s,
		const Point &	p
	)

References TRACE, QGpCoreTools::Point2D::x(), QGpCoreTools::Point2D::y(), and QGpCoreTools::Point::z().

  {
    TRACE;
    s.setRealNumberNotation(QTextStream::SmartNotation);
    s.setRealNumberPrecision(20);
    s << p.x() << " " << p.y() << " " << p.z();
    return s;
  }

QGPCORETOOLS_EXPORT QDataStream & QGpCoreTools::operator<<	(	QDataStream &	s,
		const Point &	p
	)

References TRACE, QGpCoreTools::Point2D::x(), QGpCoreTools::Point2D::y(), and QGpCoreTools::Point::z().

  {
    TRACE;
    s << p.x() << p.y() << p.z();
    return s;
  }

template<typename T >

QDataStream& QGpCoreTools::operator<<	(	QDataStream &	s,
		const Matrix< T > &	m
	)

  {
    s << m.rowCount() << m.columnCount();
    for(int i=m.rowCount()*m.columnCount()-1; i>=0; i--) {
      s << m.data()[i];
    }
    return s;
  }

template<typename T >

QDataStream& QGpCoreTools::operator<<	(	QDataStream &	s,
		const SparseMatrix< T > &	m
	)

References QGpCoreTools::SparseMatrix< T >::begin(), and QGpCoreTools::SparseMatrix< T >::end().

  {
    s << m.rowCount() << m.columnCount();
    s << m.realRowCount();
    for(typename SparseMatrix<T>::ConstRowIterator itr=m.begin(); itr!=m.end(); itr++) {
      const typename SparseMatrix<T>::Row& rowHash=itr.value();
      s << itr.key() << rowHash.count();
      for(typename SparseMatrix<T>::ConstColumnIterator itc=rowHash.begin(); itc!=rowHash.end(); itc++) {
        s<< itc.key() << itc.value();
      }
    }
    return s;
  }

template<class ValueType >

QTextStream& QGpCoreTools::operator<<	(	QTextStream &	s,
		const Grid2D< double > &	grd
	)

References endl(), and TRACE.

{
  TRACE;
  s << "x y val" << ::endl;
  const double * ptr=grd.valuePointer(0, 0);
  for(int iy=0;iy < grd.ny();iy++ ) {
    for(int ix=0;ix < grd.nx();ix++ )
      s << grd.x(ix) << " " << grd.y(iy) << " " << *(ptr++ ) << ::endl;
  }
  return s;
}

template<class ValueType >

QTextStream& QGpCoreTools::operator<<	(	QTextStream &	s,
		const Grid3D< double > &	grd
	)

References endl(), and TRACE.

{
  TRACE;
  s << "x y z val" << ::endl;
  const double * ptr=grd.valuePointer(0, 0, 0);
  for(int iz=0;iz < grd.nz();iz++ ) {
    for(int iy=0;iy < grd.ny();iy++ ) {
      for(int ix=0;ix < grd.nx();ix++ )
        s << grd.x(ix) << " " << grd.y(iy) << " " << grd.z(iz) << " " << *(ptr++ ) << ::endl;
    }
  }
  return s;
}

QGPCORETOOLS_EXPORT QTextStream & QGpCoreTools::operator<<	(	QTextStream &	s,
		const Point2D &	p
	)

References TRACE, QGpCoreTools::Point2D::x(), and QGpCoreTools::Point2D::y().

  {
    TRACE;
    s.setRealNumberNotation(QTextStream::SmartNotation);
    s.setRealNumberPrecision(20);
    s << p.x() << " " << p.y();
    return s;
  }

QGPCORETOOLS_EXPORT QDataStream & QGpCoreTools::operator<<	(	QDataStream &	s,
		const Point2D &	p
	)

References TRACE, QGpCoreTools::Point2D::x(), and QGpCoreTools::Point2D::y().

  {
    TRACE;
    s << p.x() << p.y();
    return s;
  }

QGPCORETOOLS_EXPORT QTextStream & QGpCoreTools::operator<<	(	QTextStream &	s,
		const QList< Point2D > &	plist
	)

References endl(), and TRACE.

  {
    TRACE;
    s << "%x      y"<< ::endl;
    QList<Point2D>::const_iterator it;
    for(it=plist.begin();it!=plist.end();++it) s << *it << ::endl;
    return s;
  }

QGPCORETOOLS_EXPORT QTextStream & QGpCoreTools::operator>>	(	QTextStream &	s,
		IrregularGrid2D &	grd
	)

References endl(), QGpCoreTools::Point::fromString(), QGpCoreTools::IrregularGrid2D::init(), QGpCoreTools::IrregularGrid2D::nx(), QGpCoreTools::IrregularGrid2D::ny(), QGpCoreTools::IrregularGrid2D::setValue(), QGpCoreTools::IrregularGrid2D::setX(), QGpCoreTools::IrregularGrid2D::setY(), str, QGpCoreTools::App::stream(), tr(), TRACE, and QGpCoreTools::IrregularGrid2D::valuePointer().

{
  TRACE;
  QString str;
  bool ok;
  int i;
  int line=0;
  str=s.readLine().trimmed(); line++;
  while(str.isEmpty() || str[0]=='#') {  // Accept blank lines and comments only before grid definition
    str=s.readLine().trimmed(); line++;
  }
  if(str.startsWith("nx=")) { // Historical format for grids
    int nx=str.section( "=", 1, 1).toInt(&ok);
    if( !ok || nx<=0) {
      return s;
    }
    str=s.readLine(); line++;
    int ny=str.section( "=", 1, 1).toInt(&ok);
    if( !ok || ny<=0) {
      return s;
    }
    grd.init(nx, ny);
    //  X coordinates of nodes
    str=s.readLine(); line++;
    if(str!="x" ) {
        App::stream() << IrregularGrid2D::tr("Wrong format at line 3, expected: x") << endl;
        return s;
    }
    for(i=0;i < grd.nx();i++ ) {
      str=s.readLine();
      grd.setX(i, str.section( " ", 0, 0).toDouble(&ok));
      if( !ok) {
        App::stream() << IrregularGrid2D::tr("Wrong format at line %1, expected: float value").arg(line) << endl;
        return s;
      }
    }
    // Y coordinates of nodes
    str=s.readLine(); line++;
    if(str!="y" ) {
        App::stream() << IrregularGrid2D::tr("Wrong format at line %1, expected: y").arg(line) << endl;
        return s;
    }
    for(i=0;i < grd.ny();i++ ) {
      str=s.readLine(); line++;
      grd.setY(i, str.section( " ", 0, 0).toDouble(&ok));
      if( !ok) {
        App::stream() << IrregularGrid2D::tr("Wrong format at line %1, expected: float value").arg(line) << endl;
        return s;
      }
    }
    // Values
    int n=grd.nx() * grd.ny();
    str=s.readLine().trimmed(); line++;
    if(str=="values" ) {
      double * ptr=grd.valuePointer(0, 0);
      for(i=0;i < n;i++ ) {
        str=s.readLine(); line++;
        *(ptr++ )=str.section( " ", 0, 0).toDouble(&ok);
        if( !ok) {
          App::stream() << IrregularGrid2D::tr("Wrong format at line %1, expected: float value").arg(line) << endl;
          return s;
        }
      }
    } else if(str=="binary values" ) { // No longer exported like this but kept for compatibility
      QDataStream sb(s.device());
      sb.readRawData ((char * ) grd.valuePointer(0, 0), n * sizeof(double));
      str=s.readLine(); line++;
    } else {
      App::stream() << IrregularGrid2D::tr("Wrong format at line %1, expected: values").arg(line) << endl;
    }
  } else if(str=="x y val" || // More practical format for grids
              str=="x y z" ) {  // Compatibility with snapshot release
    QList<Point> values;
    Point p;
    while( !s.atEnd()) {
      str=s.readLine().trimmed(); line++;
      if(str.isEmpty() || str[0]=='#') break;
      StringSection pStr(str);
      if( !p.fromString(pStr)) {
        App::stream() << IrregularGrid2D::tr("Wrong format at line %1, expected: float value").arg(line) << endl;
        return s;
      }
      values.append(p);
    }
    // Get nx and ny
    QMap<double,int> xMap, yMap;
    QMap<double,int>::iterator itMap;
    for(QList<Point>::iterator it=values.begin();it!=values.end(); it++) {
      itMap=xMap.find(it->x());
      if(itMap==xMap.end()) {
        xMap.insert(it->x(),0);
      }
      itMap=yMap.find(it->y());
      if(itMap==yMap.end()) {
        yMap.insert(it->y(),0);
      }
    }
    if(xMap.isEmpty() || yMap.isEmpty()) {
      return s;
    }
    grd.init(xMap.count(), yMap.count(), 0.0);
    // Set index in x and y maps
    int i;
    i=0;
    for(QMap<double,int>::iterator it=xMap.begin();it!=xMap.end(); it++, i++) {
      grd.setX(i, it.key());
      it.value()=i;
      //printf("x[%i]=%lf\n",i,it.key());
    }
    i=0;
    for(QMap<double,int>::iterator it=yMap.begin();it!=yMap.end(); it++, i++) {
      grd.setY(i, it.key());
      it.value()=i;
      //printf("y[%i]=%lf\n",i,it.key());
    }
    // Set values
    for(QList<Point>::iterator it=values.begin();it!=values.end(); it++) {
      grd.setValue(xMap[it->x()], yMap[it->y()], it->z());
    }
  } else if(str=="binary values" ) { // No longer exported like this but kept for compatibility
    QDataStream sb(s.device());
    sb.readRawData ((char * ) grd.valuePointer(0, 0), grd.nx() * grd.ny() * sizeof(double));
    str=s.readLine(); line++;
  } else {
    App::stream() << IrregularGrid2D::tr("Wrong format at line %1, expected: values").arg(line) << endl;
  }
  return s;
}

QGPCORETOOLS_EXPORT QTextStream & QGpCoreTools::operator>>	(	QTextStream &	s,
		Complex &	c
	)

References QGpCoreTools::Complex::set(), and TRACE.

{
  TRACE;
  double re, im;
  s >> re >> im;
  c.set(re, im);
  return s;
}

template<class numberType >

QDataStream & QGpCoreTools::operator>>	(	QDataStream &	s,
		StatisticalValue< numberType > &	p
	)

References QGpCoreTools::StatisticalValue< numberType >::setMean(), QGpCoreTools::StatisticalValue< numberType >::setStddev(), QGpCoreTools::StatisticalValue< numberType >::setWeight(), and TRACE.

{
  TRACE;
  numberType m;
  s >> m;
  p.setMean(m);
  double v;
  s >> v;
  p.setStddev(v);
  s >> v;
  p.setWeight(v);
  return s;
}

QGPCORETOOLS_EXPORT QDataStream & QGpCoreTools::operator>>	(	QDataStream &	s,
		Complex &	c
	)

References QGpCoreTools::Complex::set(), and TRACE.

{
  TRACE;
  double re, im;
  s >> re >> im;
  c.set(re, im);
  return s;
}

template<class ValueType >

QTextStream & QGpCoreTools::operator>>	(	QTextStream &	s,
		Grid2D< ValueType > &	grd
	)

References endl(), QGpCoreTools::Grid2D< ValueType >::init(), QGpCoreTools::StringSection::isValid(), QGpCoreTools::StringSection::nextField(), QGpCoreTools::StringSection::readLine(), QGpCoreTools::Grid2D< ValueType >::setDeltaX(), QGpCoreTools::Grid2D< ValueType >::setDeltaY(), QGpCoreTools::Grid2D< ValueType >::setOrigin(), QGpCoreTools::Grid2D< ValueType >::setValue(), QGpCoreTools::Point2D::setX(), QGpCoreTools::Point2D::setY(), str, QGpCoreTools::App::stream(), QGpCoreTools::StringSection::toDouble(), TRACE, QGpCoreTools::StringSection::trimmed(), QGpCoreTools::Grid2D< ValueType >::value(), QGpCoreTools::Point2D::x(), and QGpCoreTools::Point2D::y().

{
  TRACE;
  QString str;
  bool ok;
  int i;
  int line=0;
  str=s.readLine().trimmed(); line++;
  while(str.isEmpty() || str[0]=='#') {  // Accept blank lines and comments only before grid definition
    str=s.readLine().trimmed(); line++;
  }
  if(str=="x y val" ) {
    QVector<Point2D> points;
    QVector<ValueType> values;
    Point2D p;
    ValueType val;
    const QChar * ptr;
    StringSection f;
    bool ok=true;
    while( !s.atEnd()) {
      str=s.readLine().trimmed(); line++;
      if(str.isEmpty() || str[0]=='#') break;
      StringSection pStr(str);
      ptr=0;
      f=pStr.nextField(ptr);
      if(f.isValid()) p.setX(f.toDouble()); else ok=false;
      f=pStr.nextField(ptr);
      if(f.isValid()) p.setY(f.toDouble()); else ok=false;
      f=pStr.nextField(ptr);
      if(f.isValid()) val=f.toDouble(); else ok=false;
      if( !ok) {
        App::stream() << QCoreApplication::translate( "Grid2D", "Wrong format at line %1, expected: float value" ).arg(line) << endl;
        return s;
      }
      points.append(p);
      values.append(val);
    }
    // Deduce nx, ny, delta and origin from x and y vectors
    QMap<double,int> xMap, yMap;
    QMap<double,int>::iterator itMap;
    for(QVector<Point2D>::iterator it=values.begin();it!=values.end(); it++) {
      itMap=xMap.find(it->x());
      if(itMap==xMap.end()) {
        xMap.insert(it->x(),0);
      }
      itMap=yMap.find(it->y());
      if(itMap==yMap.end()) {
        yMap.insert(it->y(),0);
      }
    }
    if(xMap.isEmpty() || yMap.isEmpty()) {
      return s;
    }
    grd.setOrigin(Point2D(xMap.begin().value(), yMap.begin().value()));
    grd.setDeltaX(((--xMap.end()).value() - xMap.begin().value())/xMap.count());
    grd.setDeltaY(((--yMap.end()).value() - yMap.begin().value())/yMap.count());
    grd.init(xMap.count(), yMap.count(), 0.0);
    // Set index in x and y maps
    int i;
    i=0;
    for(QMap<double,int>::iterator it=xMap.begin();it!=xMap.end(); it++, i++) {
      grd.setX(i, it.key());
      it.value()=i;
      //printf("x[%i]=%lf\n",i,it.key());
    }
    i=0;
    for(QMap<double,int>::iterator it=yMap.begin();it!=yMap.end(); it++, i++) {
      grd.setY(i, it.key());
      it.value()=i;
      //printf("y[%i]=%lf\n",i,it.key());
    }
    // Set values
    for(int i=points.count()-1; i>=0; i--) {
      Point2D& p=points.at(i);
      grd.setValue(xMap[p.x()], yMap[p.y()], values.at(i));
    }
  } else {
    App::stream() << QCoreApplication::translate( "Grid2D", "Wrong format at line %1, expected: values" ).arg(line) << endl;
  }
  return s;
}

template<class ValueType >

QTextStream & QGpCoreTools::operator>>	(	QTextStream &	s,
		Grid3D< ValueType > &	grd
	)

References endl(), QGpCoreTools::Grid3D< ValueType >::init(), QGpCoreTools::StringSection::isValid(), QGpCoreTools::StringSection::nextField(), QGpCoreTools::StringSection::readLine(), QGpCoreTools::Grid3D< ValueType >::setDeltaX(), QGpCoreTools::Grid3D< ValueType >::setDeltaY(), QGpCoreTools::Grid3D< ValueType >::setDeltaZ(), QGpCoreTools::Grid3D< ValueType >::setOrigin(), QGpCoreTools::Grid3D< ValueType >::setValue(), QGpCoreTools::Point2D::setX(), QGpCoreTools::Point2D::setY(), QGpCoreTools::Point::setZ(), str, QGpCoreTools::App::stream(), QGpCoreTools::StringSection::toDouble(), TRACE, QGpCoreTools::StringSection::trimmed(), QGpCoreTools::Grid3D< ValueType >::value(), QGpCoreTools::Point2D::x(), QGpCoreTools::Point2D::y(), and QGpCoreTools::Point::z().

{
  TRACE;
  QString str;
  bool ok;
  int i;
  int line=0;
  str=s.readLine().trimmed(); line++;
  while(str.isEmpty() || str[0]=='#') {  // Accept blank lines and comments only before grid definition
    str=s.readLine().trimmed(); line++;
  }
  if(str=="x y z val" ) {
    QVector<Point> points;
    QVector<ValueType> values;
    Point p;
    ValueType val;
    const QChar * ptr;
    StringSection f;
    bool ok=true;
    while( !s.atEnd()) {
      str=s.readLine().trimmed(); line++;
      if(str.isEmpty() || str[0]=='#') break;
      StringSection pStr(str);
      ptr=0;
      f=pStr.nextField(ptr);
      if(f.isValid()) p.setX(f.toDouble()); else ok=false;
      f=pStr.nextField(ptr);
      if(f.isValid()) p.setY(f.toDouble()); else ok=false;
      f=pStr.nextField(ptr);
      if(f.isValid()) p.setZ(f.toDouble()); else ok=false;
      f=pStr.nextField(ptr);
      if(f.isValid()) val=f.toDouble(); else ok=false;
      if( !ok) {
        App::stream() << QCoreApplication::translate( "Grid3D", "Wrong format at line %1, expected: float value" ).arg(line) << endl;
        return s;
      }
      points.append(p);
      values.append(val);
    }
    // Deduce nx, ny, nz, delta and origin from x, y and z vectors
    QMap<double,int> xMap, yMap, zMap;
    QMap<double,int>::iterator itMap;
    for(QVector<Point>::iterator it=points.begin(); it!=points.end(); it++ ) {
      itMap=xMap.find(it->x());
      if(itMap==xMap.end()) {
        xMap.insert(it->x(),0);
      }
      itMap=yMap.find(it->y());
      if(itMap==yMap.end()) {
        yMap.insert(it->y(),0);
      }
      itMap=zMap.find(it->z());
      if(itMap==zMap.end()) {
        zMap.insert(it->z(),0);
      }
    }
    if(xMap.isEmpty() || yMap.isEmpty() || zMap.isEmpty()) {
      return s;
    }
    grd.setOrigin(Point(xMap.begin().value(), yMap.begin().value(), zMap.begin().value()));
    grd.setDeltaX(((--xMap.end()).value() - xMap.begin().value())/xMap.count());
    grd.setDeltaY(((--yMap.end()).value() - yMap.begin().value())/yMap.count());
    grd.setDeltaZ(((--zMap.end()).value() - zMap.begin().value())/zMap.count());
    grd.init(xMap.count(), yMap.count(), zMap.count(), 0.0);
    // Set index in x, y and z maps
    int i;
    i=0;
    for(QMap<double,int>::iterator it=xMap.begin();it!=xMap.end(); it++, i++) {
      grd.setX(i, it.key());
      it.value()=i;
      //printf("x[%i]=%lf\n",i,it.key());
    }
    i=0;
    for(QMap<double,int>::iterator it=yMap.begin();it!=yMap.end(); it++, i++) {
      grd.setY(i, it.key());
      it.value()=i;
      //printf("y[%i]=%lf\n",i,it.key());
    }
    i=0;
    for(QMap<double,int>::iterator it=zMap.begin();it!=zMap.end(); it++, i++) {
      grd.setZ(i, it.key());
      it.value()=i;
      //printf("z[%i]=%lf\n",i,it.key());
    }
    // Set values
    for(int i=points.count()-1; i>=0; i--) {
      Point& p=points.at(i);
      grd.setValue(xMap[p.x()], yMap[p.y()], zMap[p.z()], values.at(i));
    }
  } else {
    App::stream() << QCoreApplication::translate( "Grid3D", "Wrong format at line %1, expected: values" ).arg(line) << endl;
  }
  return s;
}

QGPCORETOOLS_EXPORT QTextStream & QGpCoreTools::operator>>	(	QTextStream &	s,
		QList< Point > &	plist
	)

References str, and TRACE.

  {
    TRACE;
    QString str;
    plist.clear();
    while(!s.atEnd()) {
      str=s.readLine();
      while(str[0]!='%') str=s.readLine();
      QTextStream currentLine(&str, QIODevice::ReadOnly);
      Point * p=new Point;
      currentLine >> *p;
      plist.append(*p);
    }
    return s;
  }

QGPCORETOOLS_EXPORT QTextStream & QGpCoreTools::operator>>	(	QTextStream &	s,
		Point &	p
	)

References QGpCoreTools::Point2D::setX(), QGpCoreTools::Point2D::setY(), QGpCoreTools::Point::setZ(), and TRACE.

  {
    TRACE;
    double x, y, z;
    s >> x >> y >> z;
    p.setX(x);
    p.setY(y);
    p.setZ(z);
    return s;
  }

QGPCORETOOLS_EXPORT QDataStream & QGpCoreTools::operator>>	(	QDataStream &	s,
		Point &	p
	)

References QGpCoreTools::Point2D::setX(), QGpCoreTools::Point2D::setY(), QGpCoreTools::Point::setZ(), and TRACE.

  {
    TRACE;
    double x, y, z;
    s >> x >> y >> z;
    p.setX(x);
    p.setY(y);
    p.setZ(z);
    return s;
  }

template<typename T >

QDataStream& QGpCoreTools::operator>>	(	QDataStream &	s,
		Matrix< T > &	m
	)

References QGpCoreTools::Matrix< T >::data(), and QGpCoreTools::Matrix< T >::resize().

  {
    int nrow, ncol;
    s >> nrow >> ncol;
    m.resize(nrow, ncol);
    for(int i=nrow*ncol-1; i>=0; i--) {
      s >> m.data()[i];
    }
    return s;
  }

template<typename T >

QDataStream& QGpCoreTools::operator>>	(	QDataStream &	s,
		SparseMatrix< T > &	m
	)

References QGpCoreTools::SparseMatrix< T >::addRow(), QGpCoreTools::SparseMatrix< T >::clear(), QGpCoreTools::SparseMatrix< T >::setColumnCount(), and QGpCoreTools::SparseMatrix< T >::setRowCount().

  {
    m.clear();
    int nrow, ncol;
    s >> nrow >> ncol;
    m.setRowCount(nrow);
    m.setColumnCount(ncol);
    s >> nrow;
    int index;
    T value;
    for(int ir=0; ir<nrow; ir++) {
      s >> index >> ncol;
      typename SparseMatrix<T>::Row& rowHash=m.addRow(index);
      for(int ic=0; ic<ncol; ic++) {
        s >> index >> value;
        rowHash.insert(index, value);
      }
    }
    return s;
  }

QGPCORETOOLS_EXPORT QTextStream & QGpCoreTools::operator>>	(	QTextStream &	s,
		Point2D &	p
	)

References QGpCoreTools::Point2D::setX(), QGpCoreTools::Point2D::setY(), and TRACE.

  {
    TRACE;
    double x, y;
    s >> x >> y;
    p.setX(x);
    p.setY(y);
    return s;
  }

QGPCORETOOLS_EXPORT QDataStream & QGpCoreTools::operator>>	(	QDataStream &	s,
		Point2D &	p
	)

References QGpCoreTools::Point2D::setX(), QGpCoreTools::Point2D::setY(), and TRACE.

  {
    TRACE;
    double x, y;
    s >> x >> y;
    p.setX(x);
    p.setY(y);
    return s;
  }

QGPCORETOOLS_EXPORT QTextStream & QGpCoreTools::operator>>	(	QTextStream &	s,
		QList< Point2D > &	plist
	)

References str, and TRACE.

  {
    TRACE;
    QString str;
    plist.clear();
    Point2D p;
    while(!s.atEnd()) {
      str=s.readLine();
      while(str[0]!='%') str=s.readLine();
      QTextStream currentLine(&str, QIODevice::ReadOnly);
      currentLine >> p;
      plist.append(p);
    }
    return s;
  }

QGpCoreTools::PACKAGE_INFO	(	QGpCoreTools	,
		QGPCORETOOLS
	)

QGPCORETOOLS_EXPORT uint QGpCoreTools::qHash ( StringSection key )

Global Hash function to store StringSection into QHash

This code was copied from qHash(QString key)

References QGpCoreTools::StringSection::data(), QGpCoreTools::StringSection::size(), and TRACE.

  {
    TRACE;
    const QChar * p=key.data();
    int n=key.size();
    uint h=0;
    uint g;

    while(n-- ) {
      h=(h << 4) + ( *p++ ).unicode();
      if((g=(h & 0xf0000000))!=0)
        h ^= g >> 23;
      h &= ~g;
    }
    return h;
  }

int QGpCoreTools::qHash ( const HashNumericalKey & key )

References QGpCoreTools::HashNumericalKey::hash().

Referenced by QGpCoreTools::Message::messageId(), and qHash().

{
  return key.hash();
}

QGPCORETOOLS_EXPORT uint QGpCoreTools::qHash ( Point2D key )

Global Hash function to store Point2D into QHash

Copied from uint qHash(const QPair<T1, T2> &key)

References qHash(), QGpCoreTools::Point2D::x(), and QGpCoreTools::Point2D::y().

  {
    uint h1=::qHash(quint64(key.x()));
    uint h2=::qHash(quint64(key.y()));
    return ((h1 << 16) | (h1 >> 16)) ^ h2;
  }

Complex QGpCoreTools::sin ( const Complex & c ) [inline]

References cos(), QGpCoreTools::Complex::im, and QGpCoreTools::Complex::re.

Referenced by GeopsyCore::TaperDelegate::apply(), QGpCoreWave::RefractionDippingModel::begin(), DampingResults::compute(), cos(), GeopsyCore::DoubleSignal::discreteFourierTransform(), Simulator::displaySource(), EllipseItem::EllipseItem(), exp(), QGpCoreTools::ExpressionContext::functionValue(), QGpCoreTools::Point2D::geographicalAzimuthTo(), QGpCoreTools::Point2D::geographicalDistanceTo(), QGpCoreTools::Point2D::geographicalToRectangular(), QGpCoreTools::Point2D::geographicalToUtm(), HVStationSignals::horizontal(), DampingResults::misfitFunk(), DrawScene::moveTo(), QGpCoreTools::PointLocate::position(), QGpCoreWave::RefractionDippingModel::propagate(), QGpCoreTools::Point2D::rectangularToGeographical(), QGpCoreTools::DirectionalSearch::setAzimuth(), SourceParameters::setAzimuth(), SourceParameters::setDistance(), QGpCoreTools::Complex::setExp(), QGpCoreTools::Complex::setUnitExp(), GeopsyCore::DoubleSignal::shift(), sqrt(), tan(), QGpCoreTools::Point2D::utmToGeographical(), QGpCoreWave::TheoreticalFK::value(), QGpCoreTools::Ellipse::xRadius(), QGpCoreWave::Love::y(), QGpCoreWave::RayleighTemplate< RealType >::y(), and QGpCoreTools::Ellipse::yRadius().

{
  return Complex( ::sin(c.re()) *::cosh(c.im()), ::cos(c.re()) * sinh(c.im()));
}

Complex QGpCoreTools::sqrt ( const Complex & c ) [inline]

References QGpCoreTools::Complex::abs(), cos(), QGpCoreTools::Complex::phase(), and sin().

Referenced by QGpCoreTools::Complex::abs(), GeopsyCore::DoubleSignal::agc(), asin(), QGpCompatibility::CompatFunctionList::average(), QGpCoreTools::StatisticalValue< numberType >::average(), ArrayCore::KminSolver::calculate(), GpCoreTools::Average::calculate(), QGpCoreWave::MagnetoTelluricFactory::calculate(), QGpCompatibility::CompatDispersionData::closestModeMisfit(), QGpGuiWave::MagnetoTelluricProxy::columnValue(), DampingResults::compute(), QGpCoreTools::ProcessStatistics::current(), QGpCoreTools::ProcessStatistics::currentLog(), GeopsyCore::MorletParameters::deltaF(), GeopsyCore::MorletParameters::deltaT(), QGpCoreTools::Point::distanceTo(), QGpCoreTools::Point2D::distanceTo(), QGpCoreTools::Point::distanceToSegment(), QGpCoreTools::Point::elevationTo(), QGpCoreTools::ExpressionContext::functionValue(), QGpCoreTools::Point2D::geographicalDistanceTo(), QGpCoreTools::Point2D::geographicalToUtm(), DinverDCCore::PoissonCondition::getLimits(), QGpCoreWave::Seismic1DModel::halfSpaceRayleighSlowness(), ArrayCore::StationCoupleProcess::horizontalAutocorr(), ArrayCore::HRFK::initOperator(), QGpCoreTools::Circle::intersect(), DinverDCCore::PoissonCondition::isOk(), main(), GeopsyCore::DoubleSignal::maximumAmplitude(), GeopsyCore::DoubleSignal::maximumAmplitudeAt(), PhaseShifter::maximumSlowness(), ArrayCore::FK::maximumSlowness(), QGpCoreTools::ConcentricSearch::maxRadius(), Histogram2D::meanCurve(), Histogram2D::medianCurve(), Spac3CForward::misfit(), QGpCompatibility::CompatDispersionData::misfit(), QGpCompatibility::CompatModalEllipticity::misfit(), Histogram2D::modeCurve(), LinearFKActiveStationSignals::normalize(), GeopsyCore::DoubleSignal::normalizedCorrelation(), StatGridAnalyser::on_freqScroll_valueChanged(), QGpCoreTools::PointLocate::position(), QGpCoreTools::Curve< pointType >::project(), QGpCoreWave::RefractionDippingModel::propagate(), QGpCoreTools::Curve< pointType >::resample(), QGpCompatibility::CompatFunction::resample_ds(), ArrayCore::FKRadial::rotationFactors(), ArrayCore::FKTransverse::rotationFactors(), SPACLoopTask::run(), QGpCoreTools::Angle::set(), MonoStation::StationResults::setAverageLayer(), MonoStation::AbstractSummary::setAverageValues(), GeopsyCore::ComplexSignal::sqrt(), GeopsyCore::DoubleSignal::sqrt(), QGpCompatibility::CompatInversionReport::statDispersion(), GpCoreTools::LinearRegression::stddev(), GpCoreTools::Statistics::stddev(), QGpCoreTools::Statistics::stddev(), QGpCoreTools::ProcessStatistics::stddev(), QGpCompatibility::CompatFunction::stddev(), QGpCoreTools::Covariance::stddev2D(), GeopsyCore::DoubleSignal::stddevClip(), QGpCoreTools::ProcessStatistics::stddevLog(), MonoStation::StatisticResults::studentTest(), DinverDCCore::DCReportBlock::surfaceWaveModel(), CoordReader::terminate(), GeopsyCore::DoubleSignal::timeCorrelation(), QGpCoreWave::MagnetoTelluricPointOptions::toAbsoluteValue(), GeopsyCore::DoubleSignal::unglitch(), QGpCoreTools::Point2D::utmToGeographical(), ArrayCore::StationCoupleProcess::verticalAutocorr(), QGpCoreWave::Love::y(), QGpCoreWave::RayleighTemplate< RealType >::y(), and QGpCoreTools::Curve< pointType >::ySqrt().

{
  double r2=::sqrt(c.abs());
  double phi2=c.phase()/2;
  return Complex(r2*::cos(phi2), r2*::sin(phi2));
}

Complex QGpCoreTools::tan ( const Complex & c ) [inline]

References cos(), and sin().

Referenced by QGpCoreTools::ExpressionContext::functionValue(), QGpCoreTools::Point2D::geographicalToRectangular(), QGpCoreTools::Line2D::Line2D(), and QGpCoreTools::Point2D::rectangularToGeographical().

{
  return sin(c)/cos(c);
}

QString QGpCoreTools::tr	(	const char *	s,
		const char *	comments = `0`
	)		`[inline]`

{
  return QCoreApplication::translate("global", s, comments);
}

template<typename T >

QGpCoreTools::unique ( QList< T > & list )

Removes repeated entries. The container must be sorted.

Referenced by GeopsyCore::SubSignalPool::availableTimePicks(), QGpCoreTools::Curve< pointType >::average(), RecordList::contracts(), DinverDCGui::GroundModelViewer::exportModels(), QGpGuiTools::XMLItem::getSelection(), ToolRefra::initStations(), DinverDCCore::TargetList2D::misfit(), outputDCModel(), MainEnvironment::projects(), MainWindow::projects(), DinverDCCore::DCReportBlock::readSeismicModel(), MonoStation::WinResults::setCategories(), MonoStation::AbstractSummary::setStations(), QGpCoreWave::MagnetoTelluricFactory::setX(), QGpCoreWave::ModalFactory::setX(), QGpCoreWave::RefractionFactory::setX(), QGpCompatibility::CompatFunction::sort(), QGpCoreTools::Curve< pointType >::unique(), and GraphicWindow::updateLabels().

  {
    for(int i=1; i<list.count(); ) {
      if(list.at(i-1)==list.at(i)) {
        list.removeAt(i);
      } else {
        i++;
      }
    }
  }

template<typename T , typename Equal >

void QGpCoreTools::unique	(	QList< T > &	list,
		Equal	equal
	)

  {
    for(int i=1; i<list.count(); ) {
      if(equal(list.at(i-1), list.at(i))) {
        list.removeAt(i);
      } else {
        i++;
      }
    }
  }

template<typename T >

QGpCoreTools::unique ( QVector< T > & v )

Removes repeated entries. The container must be sorted.

  {
    for(int i=1; i<v.count(); ) {
      if(v.at(i-1)==v.at(i)) {
        v.remove(i);
      } else {
        i++;
      }
    }
  }

template<typename T , typename Equal >

void QGpCoreTools::unique	(	QVector< T > &	v,
		Equal	equal
	)

  {
    for(int i=1; i<v.count(); ) {
      if(equal(v.at(i-1), v.at(i))) {
        v.remove(i);
      } else {
        i++;
      }
    }
  }

template<typename T >

void QGpCoreTools::uniqueDeleteAll ( QList< T > & list )

Referenced by TapePositioningSystem::Triangulator::aggregate().

  {
    for(int i=1; i<list.count(); ) {
      if(list.at(i-1)==list.at(i)) {
        delete list.takeAt(i);
      } else {
        i++;
      }
    }
  }

template<typename T , typename Equal >

void QGpCoreTools::uniqueDeleteAll	(	QList< T > &	list,
		Equal	equal
	)

  {
    for(int i=1; i<list.count(); ) {
      if(equal(list.at(i-1), list.at(i))) {
        delete list.takeAt(i);
      } else {
        i++;
      }
    }
  }

template<typename T >

void QGpCoreTools::uniqueDeleteAll ( QVector< T > & v )

  {
    for(int i=1; i<v.count(); ) {
      if(v.at(i-1)==v.at(i)) {
        delete v.at(i);
        v.remove(i);
      } else {
        i++;
      }
    }
  }

template<typename T , typename Equal >

void QGpCoreTools::uniqueDeleteAll	(	QVector< T > &	v,
		Equal	equal
	)

  {
    for(int i=1; i<v.count(); ) {
      if(equal(v.at(i-1), v.at(i))) {
        delete v.at(i);
        v.remove(i);
      } else {
        i++;
      }
    }
  }

Classes

Typedefs

Enumerations

Functions

Detailed Description

Typedef Documentation

Enumeration Type Documentation

Function Documentation