fix the coherent partial envelop function for TEM

Ivanlh20 · Jan 10, 2023 · 3f55109 · 3f55109
1 parent be8903f
commit 3f55109
Show file tree

Hide file tree

Showing 47 changed files with 48 additions and 48 deletions.
diff --git a/matlab_functions/ilm_read_ap_cif.m b/matlab_functions/ilm_read_ap_cif.m
@@ -1,5 +1,5 @@
 % Read atomic positions from cif file
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 function [atoms, lx, ly, lz] = ilm_read_ap_cif(path, rmsd_3d_0, pbc, na, nb, nc)
     if(nargin<6)
         nc = 1;

diff --git a/matlab_functions/ilm_read_lat_parm_cif.m b/matlab_functions/ilm_read_lat_parm_cif.m
@@ -1,5 +1,5 @@
 % Read lattice parameter from cif file
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 function [a, b, c] = ilm_read_lat_parm_cif(path)
     str_file = fileread(path);
     str_file = strtrim(strsplit(str_file, '\n'));

diff --git a/mex_examples_multem/example_MULTEM_CBED.m b/mex_examples_multem/example_MULTEM_CBED.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % Convergent beam electron diffraction (CBED) simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_MULTEM_CBEI.m b/mex_examples_multem/example_MULTEM_CBEI.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % Convergent beam electron Imaging (CBEI) simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_MULTEM_ED.m b/mex_examples_multem/example_MULTEM_ED.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % Electron diffraction (ED) simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_MULTEM_EELS.m b/mex_examples_multem/example_MULTEM_EELS.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % STEM electron energy loss spectroscopy (EELS) simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_MULTEM_EFTEM.m b/mex_examples_multem/example_MULTEM_EFTEM.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % Energy filtered transmission electron microscopy (EFTEM) simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_MULTEM_EWFS.m b/mex_examples_multem/example_MULTEM_EWFS.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % Exit wave real space (EWRS) simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_MULTEM_EWFS_till_illumination.m b/mex_examples_multem/example_MULTEM_EWFS_till_illumination.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % Exit wave real space (EWRS) simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_MULTEM_EWRS.m b/mex_examples_multem/example_MULTEM_EWRS.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % Exit wave real space (EWRS) simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_MULTEM_EWRS_till_illumination.m b/mex_examples_multem/example_MULTEM_EWRS_till_illumination.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % Exit wave real space (EWRS) simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_MULTEM_HCTEM.m b/mex_examples_multem/example_MULTEM_HCTEM.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % Hollow cone transmission electron microscopy (HCTEM) simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_MULTEM_HRTEM.m b/mex_examples_multem/example_MULTEM_HRTEM.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % High resolution transmission electron microscopy (HRTEM) simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_MULTEM_ISTEM.m b/mex_examples_multem/example_MULTEM_ISTEM.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % Imaging scanning transmission electron microscopy (ISTEM) simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_MULTEM_PED.m b/mex_examples_multem/example_MULTEM_PED.m
@@ -4,7 +4,7 @@
 % 
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
 % 
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_MULTEM_STEM.m b/mex_examples_multem/example_MULTEM_STEM.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % Scanning transmission electron microscopy (STEM) simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_MULTEM_STEM_EELS.m b/mex_examples_multem/example_MULTEM_STEM_EELS.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % STEM electron energy loss spectroscopy (EELS) simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_MULTEM_STEM_matrix_detector.m b/mex_examples_multem/example_MULTEM_STEM_matrix_detector.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % Scanning transmission electron microscopy (STEM) simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_add_amorphous_layer.m b/mex_examples_multem/example_add_amorphous_layer.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_amorphous_specimen.m b/mex_examples_multem/example_amorphous_specimen.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_apply_ctf.m b/mex_examples_multem/example_apply_ctf.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % inclusion of TEM aberrerations into a exitwve simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_beam_self_interaction.m b/mex_examples_multem/example_beam_self_interaction.m
@@ -1,6 +1,6 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_crystal_by_layers.m b/mex_examples_multem/example_crystal_by_layers.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_feg.m b/mex_examples_multem/example_feg.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_fxg.m b/mex_examples_multem/example_fxg.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_incident_wave_1.m b/mex_examples_multem/example_incident_wave_1.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % Incident wave simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_incident_wave_2.m b/mex_examples_multem/example_incident_wave_2.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % Incident wave simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_incident_wave_3.m b/mex_examples_multem/example_incident_wave_3.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % Incident wave simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_incident_wave_4.m b/mex_examples_multem/example_incident_wave_4.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % Incident wave simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) '/mex_bin'])

diff --git a/mex_examples_multem/example_microscope_aberrations.m b/mex_examples_multem/example_microscope_aberrations.m
@@ -1,7 +1,7 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % inclusion of TEM aberrerations into a exitwve simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_pr.m b/mex_examples_multem/example_pr.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_projected_potential.m b/mex_examples_multem/example_projected_potential.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_propagate.m b/mex_examples_multem/example_propagate.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_specimen_creation_from_cif_file.m b/mex_examples_multem/example_specimen_creation_from_cif_file.m
@@ -1,5 +1,5 @@
 % Specimen creation
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear all;clc;
 

diff --git a/mex_examples_multem/example_specimen_planes.m b/mex_examples_multem/example_specimen_planes.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_specimen_rotation.m b/mex_examples_multem/example_specimen_rotation.m
@@ -4,7 +4,7 @@
 % 
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
 % 
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_specimen_slicing_rotation.m b/mex_examples_multem/example_specimen_slicing_rotation.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_specimen_slicing_show.m b/mex_examples_multem/example_specimen_slicing_show.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_specimen_slicing_with_amorphous_layer.m b/mex_examples_multem/example_specimen_slicing_with_amorphous_layer.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_specimen_slicing_zone_axis.m b/mex_examples_multem/example_specimen_slicing_zone_axis.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_specimen_subslicing.m b/mex_examples_multem/example_specimen_subslicing.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_transmission_function.m b/mex_examples_multem/example_transmission_function.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_vp.m b/mex_examples_multem/example_vp.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_vr.m b/mex_examples_multem/example_vr.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_vz.m b/mex_examples_multem/example_vz.m
@@ -1,4 +1,4 @@
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>
 
 clear; clc;
 addpath([fileparts(pwd) filesep 'mex_bin'])

diff --git a/mex_examples_multem/example_wave_function.m b/mex_examples_multem/example_wave_function.m
@@ -1,6 +1,6 @@
 % output_multislice = input_multem.ilc_multem perform TEM simulation
 % All parameters of the input_multem structure are explained in ilm_dflt_input_multem()
-% Copyright 2021 Ivan Lobato <Ivanlh20@gmail.com>clear; clc;
+% Copyright 2023 Ivan Lobato <Ivanlh20@gmail.com>clear; clc;
 
 addpath([fileparts(pwd) filesep 'mex_bin'])
 addpath([fileparts(pwd) filesep 'crystalline_materials'])

diff --git a/src/cgpu_fcns.cuh b/src/cgpu_fcns.cuh
@@ -2019,10 +2019,10 @@ namespace mt
 			{			
 				T_r chi = g2*(lens.c_c_30*g2+lens.c_c_10);
 				T_r c = c_Pi*lens.si_theta_c*lens.ti_iehwgd;
-				T_r u = 1.0 + 2*c*c*g2;
+				T_r u = 1.0 + c*c*g2;
 
 				c = c_Pi*lens.ti_iehwgd*lens.lambda*g2;
-				T_r temp_inc = 0.5*c*c;
+				T_r temp_inc = 0.25*c*c;
 
 				c = c_Pi*lens.si_theta_c*(lens.c_30*lens.lambda2*g2-lens.c_10);
 				T_r spa_inc = c*c*g2;