src/inversion/Model.h

Go to the documentation of this file.

/*
 * Model.h
 *
 *  Created on: Aug 18, 2010
 *      Author: Brenda Javornik (e-mail: brendaj (at) nwra.com)
 */

#ifndef MODEL_H_
#define MODEL_H_

#include "MapDepthMatrix.h"
#include "FITSParameters.h"
#include "BaseDirectory.h"

typedef struct {
        MAP_DEPTH_MATRIX w;                     
        MAP_DEPTH_MATRIX q;                     
        MAP_DEPTH_MATRIX AvgKern;       
        MAP_DEPTH_MATRIX noise;         
        MAP_DEPTH_MATRIX dTilde;        
        MAP_DEPTH_MATRIX bFull;         
        MATRIX kxGridValues;            
        MATRIX kyGridValues;            
        bool kxkyDone;                      
        MATRIX lambdaGridValues;        

        bool convergence;                   
        bool noMaps;                    

        bool xt;                                        

        MATRIX bVECTORS;                
        MATRIX A1d;                             
        MATRIX Cov1d;                   

        Float mu;               
        Float sigma;   
        Float nu;      
        FITSParameters parms; 
} Model;

void initializeModel(Model *model, int X, int Y, int M, int Z);

void initializeOutputModel(Model *model, MAP_DEPTH_MATRIX kernelTdm);

void initializeForwardOutputModel(Model *model, MAP_DEPTH_MATRIX kernelTdm);

void initializeCalcQOutputModel(Model *model, MAP_DEPTH_MATRIX w);

void useNoMaps(Model *model);

void setParameters(Model *model, Float mu, Float sigma);

void setParameters2(Model *model, Float *values, int nValues);

void clearParameters(Model *model);

void saveW(Model *model, MATRIX w, int kx, int ky);

void saveB(Model *model, MATRIX b, int nStep);

void saveGridValues(Model *m, int r, int c, Float kx, Float ky, Float lambda);

void calculateQValues(Model *model, MAP_DEPTH_MATRIX tau, bool DValueShift);

void calculateDerivedValues(Model *model, MAP_DEPTH_MATRIX K, MAP_DEPTH_MATRIX tau, MAP_DEPTH_MATRIX Lambda,
                bool DValueShift);

void calculateDerivedValuesNoQ(Model *model, MAP_DEPTH_MATRIX K, MAP_DEPTH_MATRIX Lambda,
                bool DValueShift);

void writeQ(Model *model, BaseDirectory *baseName);

void writeW(Model *model, BaseDirectory *baseName);

void writeAvgKern(Model *model, BaseDirectory *baseName);

void writeBVectors1d(Model *model, BaseDirectory *baseName);

void writeAVectors1d(Model *model, BaseDirectory *baseName, Float lambda);

void writeCovVectors1d(Model *model, BaseDirectory *baseName, Float lambda);

void writeNoise(Model *model, BaseDirectory *baseName);

void writeDTildaValues3d(Model *model, BaseDirectory *baseName, bool useFT);

void writeAvgKernValues3d(Model *model, BaseDirectory *baseName);

void writeGridValues3d(Model *model, BaseDirectory *baseName);

void writeBVectors3d(Model *model, BaseDirectory *baseName);

void freeModel(Model *model);

#endif /* MODEL_H_ */