s_bootstrap.C File Reference

Command line program to do paleomag s data processing into xyz. More...

#include <gsl/gsl_rng.h>
#include <gsl/gsl_randist.h>
#include <iostream>
#include <iomanip>
#include <fstream>
#include <string>
#include <sstream>
#include "kdsPmagL.H"
#include "SiteSigma.H"
#include "Bootstrap.H"
#include "VecAngle.H"
#include "debug.H"
#include "s_bootstrap_cmd.h"
#include "Eigs.H"

Enumerations

enum  BootTypeEnum { BAD_PARAMETRIC, SITE_PARAMETRIC, SAMPLE_PARAMETRIC }
 What kind of boot strap to do? More...
enum  FormatEnum { BAD_FORMAT, XYZ_FORMAT, TPR_FORMAT, S_FORMAT }

Functions

bool LoadS (const string filename, vector< SVec > &s, vector< float > &sigmas)
 load ascii whitespace delimited text into vectors.
bool GetFiles (const size_t numArg, char **in_arg, vector< string > &inFiles)
BootTypeEnum GetParametricType (const int site_given, const int sample_given)
FormatEnum GetFormat (const char *format_arg)
bool DoS_Bootstrap (const vector< string > &inFiles, ofstream &out1Max, ofstream &out2Int, ofstream &out3Min, const int numout_arg, const FormatEnum format, const BootTypeEnum type, const int draw)
 Actually do the boot strap.
int main (const int argc, char *argv[])

Variables

int debug_level
 The main program needs to define this as a non-extern.
static const UNUSED char * RCSid = "@(#) $Id: s_bootstrap.C,v 1.21 2004/11/04 18:26:51 schwehr Exp $"
 Let the debugger find out which version is being used.

Detailed Description

Command line program to do paleomag s data processing into xyz.

Enumeration Type Documentation

enum BootTypeEnum
 

What kind of boot strap to do?

Site parametic looks at all the data and uses the Hext method to create one sigma that is applied to all samples.

Sample parametric uses the sigma at the end of each sample (possition 7) to make a new sample that particular measurement.

Enumerator:
BAD_PARAMETRIC 
SITE_PARAMETRIC 
SAMPLE_PARAMETRIC 

enum FormatEnum
 

Enumerator:
BAD_FORMAT 
XYZ_FORMAT 
TPR_FORMAT 
S_FORMAT 

Function Documentation

bool DoS_Bootstrap (  const vector< string > &  inFiles,
ofstream &  out1Max,
ofstream &  out2Int,
ofstream &  out3Min,
const int  numout_arg,
const FormatEnum  format,
const BootTypeEnum  type,
const int  draw
) 
 

Actually do the boot strap.

Parameters:
inFiles vector of files to read in and bootstrap
out1Max,out2Int,out3Min Each of the streams to write to. vmax, vint, vmin. Make them all the same and set oneFile true to get one file will all 9 parameters
numout_arg 1 if out1, out2, and out3 are all the same file, otherwise should be 3
format How do we want the output to look. (S, XYZ, other some other day)
type PARAMETRIC_SITE or PARAMETRIC_SAMPLE
draw How many sample to draw out of the magic hat

bool GetFiles (  const size_t  numArg,
char **  in_arg,
vector< string > &  inFiles
) 
 

FormatEnum GetFormat (  const char *  format_arg  ) 
 

BootTypeEnum GetParametricType (  const int  site_given,
const int  sample_given
) 
 

bool LoadS (  const string  filename,
vector< SVec > &  s,
vector< float > &  sigmas
) 
 

load ascii whitespace delimited text into vectors.

Returns:
true if all went well. false if trouble of any kind
Parameters:
filename File to open and read data from
s Return vector of data. The s diagonalized matrix parameters. See s_eigs
sigmas Return vector of sigma errors
Unlike Lisa's code, this one does NOT alter the sigmas on loading which is what the adread subroutine did.` You must call SiteSigma if doing a Site based Parametric Bootstrap

int main (  const int  argc,
char *  argv[]
) 
 

Variable Documentation

int debug_level
 

The main program needs to define this as a non-extern.

const UNUSED char* RCSid = "@(#) $Id: s_bootstrap.C,v 1.21 2004/11/04 18:26:51 schwehr Exp $" [static]
 

Let the debugger find out which version is being used.

Generated on Thu Mar 2 17:55:25 2006 for density by  doxygen 1.4.6