Description of fwd

DESCRIPTION

 FWD_SOLVE: calculate data from a fwd_model object and an image
 
 fwd_solve can be called as
    data= fwd_solve( fwd_model, img)
 or
    data= fwd_solve( img)

 in each case it will call the fwd_model.solve
                        or img.fwd_model.solve method

 For reconstructions on dual meshes, the interpolation matrix
    is defined as fwd_model.coarse2fine. If required, this takes
    coarse2fine * x_coarse = x_fine

 data      is a measurement data structure
 fwd_model is a fwd_model structure
 img       is an img structure

 Options: (not available on all solvers)
    img.fwd_solve.get_all_meas = 1 (data.volt = all FEM nodes, but not CEM)
    img.fwd_solve.get_all_nodes= 1 (data.volt = all nodes, including CEM)

CROSS-REFERENCE INFORMATION

eidors_model_params mdl = eidors_model_params( mdl );

eidors_msg EIDORS_MSG: eidors progress and status messages

eidors_obj EIDORS_OBJ: 'constructor' to create a eidors structure

GN_abs_solve GN_ABS_SOLVER absolute solver using Gauss Newton approximation

aa_calc_jacobian AA_CALC_JACOBIAN: J= aa_calc_jacobian( fwd_model, img)

aa_e_move_jacobian AA_E_MOVE_JACOBIAN: J= aa_e_move_jacobian( fwd_model, img)

calc_all_node_voltages CALC_ALL_NODE_VOLTAGES - calculate voltage on all nodes

calc_noise_figure CALC_NOISE_FIGURE: calculate the noise amplification (NF) of an algorithm

perturb_jacobian PERTURB_JACOBIAN: J= perturb_jacobian( fwd_model, img)

ab_tv_diff_solve AB_TV_DIFF_SOLVE inverse solver for Andrea Borsic's

pdipm_abs PDIPM_ABS inverse solver for absolute data using Primal/Dual interior point method

primaldual_tvrecon_lsearch [rs,x]=primaldual_tvrecon_lsearch(inv_mdl, vmeas, ...

inv_solve_TSVD INV_SOLVE_TSVD: inverse solver based on truncatated SVD

calc_transferimpedance CALC_TRANSFERIMPEDANCE: Calculates transfer impedance matrix

bs_nonlinearGN NONLINEARGN Non-Linear EIT Inverse

calc_move_jacobian CALC_JACOBIAN Computes the Jacobian matrix for conductivity and

calc_meas_icov meas_icov = calc_meas_icov( inv_model )

np_calc_jacobian NP_CALC_JACOBIAN: J= np_calc_jacobian( fwd_model, img)

np_inv_solve NP_INV_SOLVE inverse solver for Nick Polydorides EIDORS3D code

cheating_2d code to simulate inverse crimes in EIT

compare_2d_algs Compare different 2D reconstructions

compare_3d_algs Compare different 3D reconstructions

demo_2d_simdata Example of using EIDORS to simulate 2D data and to

demo_3d_simdata How to make simulation data using EIDORS3D

demo_real [inhomg_img, demo_img] = demo_real;

eidors2d_demo1 EidorsDemo1 Demonstrates the use of 2D EIT Package with linear basis

manchester_tomography Example to show reconstructions from

object_in_tank_2d 2D demo example for reconstruction of object floating inside tank with

resistor_model DEMO to show really simple application of EIDORS framework

show_current SHOW_CURRENT: show current or other quantity defined

mk_GREIT_model MK_GREIT_MODEL: make EIDORS inverse models using the GREIT approach

simulate_2d_movement SIMULATE_2D_MOVEMENT simulate rotational movement in 2D

simulate_3d_movement SIMULATE_3D_MOVEMENT simulate rotational movement in 2D

simulate_movement SIMULATE_MOVEMENT simulate small conductivity perturbations

stim_meas_list STIM_MEAS_LIST: mk stimulation pattern from list of electrodes

test_GREIT_model

calc_data_prior_test Verify dataprior:

calc_jacobian_test Verify Jacobian Calculation by small derivative from forward problem

demo_real_test2 Perform tests based on the demo_real function with new structs

demo_real_test3 Perform tests based on the demo_real function with new structs

perturb_jacobian_test Perturbation Jacobians

test_2d_resistor Create 2D model of a cylindrical resistor

test_3d_resistor Create 3D model of a Rectangular resistor

test_colour_direction Test colour mapping

test_performance TEST_PERFORMANCE: test of difference reconstruction algorithms

test_performance_img TEST_PERFORMANCE: test of difference reconstruction algorithms

SOURCE CODE

0001 function data = fwd_solve( fwd_model, img)
0002 % FWD_SOLVE: calculate data from a fwd_model object and an image
0003 %
0004 % fwd_solve can be called as
0005 %    data= fwd_solve( fwd_model, img)
0006 % or
0007 %    data= fwd_solve( img)
0008 %
0009 % in each case it will call the fwd_model.solve
0010 %                        or img.fwd_model.solve method
0011 %
0012 % For reconstructions on dual meshes, the interpolation matrix
0013 %    is defined as fwd_model.coarse2fine. If required, this takes
0014 %    coarse2fine * x_coarse = x_fine
0015 %
0016 % data      is a measurement data structure
0017 % fwd_model is a fwd_model structure
0018 % img       is an img structure
0019 %
0020 % Options: (not available on all solvers)
0021 %    img.fwd_solve.get_all_meas = 1 (data.volt = all FEM nodes, but not CEM)
0022 %    img.fwd_solve.get_all_nodes= 1 (data.volt = all nodes, including CEM)
0023 
0024 % (C) 2005 Andy Adler. License: GPL version 2 or version 3
0025 % $Id: fwd_solve.html 2819 2011-09-07 16:43:11Z aadler $
0026 
0027 if nargin==1
0028    img= fwd_model;
0029    fwd_model= img.fwd_model;
0030 end
0031 fwd_model= eidors_model_params( fwd_model );
0032 
0033 if isfield(fwd_model,'coarse2fine')
0034    c2f= fwd_model.coarse2fine;
0035    if size(img.elem_data,1)==size(c2f,2)
0036 %     fwd_model data is provided on coarse mesh
0037       img.elem_data = c2f * img.elem_data; 
0038    end
0039 end
0040 
0041 if isfield(fwd_model,'background')
0042     img.elem_data = img.elem_data + fwd_model.background; 
0043 end
0044 
0045 if ~isfield(fwd_model, 'electrode')
0046    error('EIDORS: attempting to solve on model without electrodes');
0047 end
0048 if ~isfield(fwd_model, 'stimulation')
0049    error('EIDORS: attempting to solve on model without stimulation patterns');
0050 end
0051 
0052 data = feval( fwd_model.solve, fwd_model, img);
0053 data= eidors_obj('data',data);  % create data object
0054 
0055 eidors_obj('set-cache', img, 'fwd_solve_data', data);
0056 eidors_msg('fwd_solve: setting cached value',3);

fwd_solve

PURPOSE

SYNOPSIS

DESCRIPTION

CROSS-REFERENCE INFORMATION

SOURCE CODE