diff options
| author | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-14 21:23:57 +0000 |
|---|---|---|
| committer | Eugeniy Mikhailov <evgmik@gmail.com> | 2009-12-14 21:23:57 +0000 |
| commit | f585ad0d4d495312d41d759e8695c46d1c7162ef (patch) | |
| tree | 1a1d8aa4ade5d38a152da44ab01ed7f973918117 /useful_functions.m | |
| parent | cd1f205ca2b52bd90d2c523cfa76ebd4a15505de (diff) | |
| download | multi_mode_eit-f585ad0d4d495312d41d759e8695c46d1c7162ef.tar.gz multi_mode_eit-f585ad0d4d495312d41d759e8695c46d1c7162ef.zip | |
now I use parcellfun for loop evaluation, and storage of the global variables in file, so algorithm is less memory hungry.
Diffstat (limited to 'useful_functions.m')
| -rw-r--r-- | useful_functions.m | 9 |
1 files changed, 6 insertions, 3 deletions
diff --git a/useful_functions.m b/useful_functions.m index e0c9e8d..bc8f709 100644 --- a/useful_functions.m +++ b/useful_functions.m @@ -236,13 +236,16 @@ endfunction function xi=susceptibility_steady_state_at_freq( atom_field_problem) % find steady state susceptibility at particular modulation frequency element % at given E_field - L0m = atom_field_problem.L0m ; - polarizability_m = atom_field_problem.polarizability_m ; - dipole_elements = atom_field_problem.dipole_elements ; + %L0m = atom_field_problem.L0m ; + %polarizability_m = atom_field_problem.polarizability_m ; + %dipole_elements = atom_field_problem.dipole_elements ; E_field = atom_field_problem.E_field ; modulation_freq = atom_field_problem.modulation_freq ; freq_index = atom_field_problem.freq_index ; + % load from file values of 'L0m', 'polarizability_m', 'dipole_elements'; + load( atom_field_problem.atom_properties_fname ); + rhoLiouville=rhoLiouville_steady_state(L0m, polarizability_m, E_field, modulation_freq); xi=susceptibility(freq_index, rhoLiouville, dipole_elements); endfunction |