num_elems

PURPOSE

SYNOPSIS

DESCRIPTION

 NUM_ELEMS: number of elemnts in a (fwd or inv model or image)
 num_elems = num_elems( mdl );

CROSS-REFERENCE INFORMATION

• mk_common_model MK_COMMON_MODEL: make common EIT models
• mk_image MK_IMAGE: create eidors image object
• num_elems NUM_ELEMS: number of elemnts in a (fwd or inv model or image)
• unit_test_cmp UNIT_TEST_CMP: compare matrices in eidors output

• show_3d_slices show_3d_slices(img, z_cuts, x_cuts, y_cuts, any_cuts)
• show_fem SHOW_FEM: show the EIDORS3D finite element model
• show_fem_enhanced SHOW_FEM_ENHANCED: show the EIDORS3D finite element model
• mat_idx_to_electrode MAT_IDX_TO_ELECTRODE: create electrodes from mat_idx values
• ng_write_opt NG_WRITE_OPT Write an ng.opt file in current directory
• calc_model_reduction calc_model_reduction: calculate the fields for a reduced model
• fix_model FIX_MODEL: Add useful fields to a model
• get_elem_volume GET_ELEM_VOLUME: VOL = get_elem_volume(fwd_model, map_node )
• interp_mesh INTERP_MESH: calculate interpolation points onto mdl elements
• linear_reorder [fwd_model] = linear_reorder(fwd_model,ccw)
• mk_c2f_circ_mapping MK_C2F_CIRC_MAPPING: create a mapping matrix from circles/spheres to FEM
• mk_grid_c2f MK_GRID_C2F - calculate a coarse2fine mapping for grid coarse models.
• mk_image MK_IMAGE: create eidors image object
• num_elems NUM_ELEMS: number of elemnts in a (fwd or inv model or image)
• point_in_tet POINT_IN_TET test for points contained in elements
• remove_unused_nodes REMOVE_UNUSED_NODES: identify and remove unused nodes in model
• jacobian_adjoint JACOBIAN_ADJOINT: J= jacobian_adjoint( img )
• jacobian_adjoint_2p5d_1st_order JACOBIAN_ADJOINT_2P5D: J= jacobian_adjoint_2p5d_1st_order( img )
• jacobian_movement_2p5d_1st_order JACOBIAN_MOVEMENT_2P5D: J = jacobian_movement_2p5d_1st_order( img )
• system_mat_1st_order SYSTEM_MAT_1ST_ORDER: SS= system_mat_1st_order( fwd_model, img)
• system_mat_2p5d_1st_order SYSTEM_MAT_2P5D_1ST_ORDER: 2.5D system matrix
• system_mat_fields SYSTEM_MAT_FIELDS: fields (elem to nodes) fraction of system mat
• GREIT_desired_img_FEMmesh GREIT_DESIRED_IMG_FEMmesh GREIT onto a FEM mesh
• calc_solution_error CALC_SOLUTION_ERROR Calculate residuals for a solution
• prior_gaussian_HPF PRIOR_GAUSSIAN_HPF calculate image prior
• test_2d_resistor Create 2D model of a cylindrical resistor

SUBFUNCTIONS

• function do_unit_test

SOURCE CODE

0001 function num = num_elems( mdl );
0002 % NUM_ELEMS: number of elemnts in a (fwd or inv model or image)
0003 % num_elems = num_elems( mdl );
0004 
0005 % $Id: num_elems.m 5112 2015-06-14 13:00:41Z aadler $
0006 
0007 if ischar(mdl) && strcmp(mdl,'UNIT_TEST'); do_unit_test; return; end
0008 
0009 if ~isfield(mdl,'type') && isfield(mdl,'elems') % Set default for this case
0010    mdl.type = 'fwd_model';
0011 end
0012 
0013 switch mdl.type
0014   case 'image';      fmdl = mdl.fwd_model;
0015   case 'inv_model';  fmdl = mdl.fwd_model;
0016   case 'fwd_model';  fmdl = mdl;
0017   otherwise;
0018       error('can''t process model of type %s', mdl.type );
0019 end
0020 
0021 num = size(fmdl.elems,1);
0022 
0023 
0024 function do_unit_test
0025    mdl = mk_common_model('a2c2',8);
0026    ne = num_elems( mdl );
0027    unit_test_cmp('test1',ne, 64);
0028 
0029    ne = num_elems( mdl.fwd_model );
0030    unit_test_cmp('test2',ne, 64);
0031 
0032    ne = num_elems( mk_image( mdl ));
0033    unit_test_cmp('test3',ne, 64);
0034 
0035    mdl = mk_common_model('n3r2',[16,2]);
0036    ne = num_elems( mk_image( mdl ));
0037    unit_test_cmp('test4',ne, 828);