1 | % ---------------------------------------------------------------------- |
---|
2 | % EXAMPLE: Fermi-Dirac function in Octave. |
---|
3 | % |
---|
4 | % This script represents a newly written application with rappture |
---|
5 | % bindings and interface. |
---|
6 | % |
---|
7 | % ====================================================================== |
---|
8 | % AUTHOR: Michael McLennan, Purdue University |
---|
9 | % AUTHOR: Derrick Kearney, Purdue University |
---|
10 | % Copyright (c) 2004-2005 Purdue Research Foundation |
---|
11 | % |
---|
12 | % See the file "license.terms" for information on usage and |
---|
13 | % redistribution of this file, and for a DISCLAIMER OF ALL WARRANTIES. |
---|
14 | % ====================================================================== |
---|
15 | |
---|
16 | % get input file from the command line |
---|
17 | % open our xml input file. |
---|
18 | lib = rpLib(infile); |
---|
19 | |
---|
20 | % retrieve user specified data out of the input file |
---|
21 | % convert values to correct units. |
---|
22 | Ef = rpLibGetString(lib,"input.number(Ef).current"); |
---|
23 | [Ef,err] = rpUnitsConvertDbl(Ef,"eV"); |
---|
24 | T = rpLibGetString(lib,"input.number(temperature).current"); |
---|
25 | [T,err] = rpUnitsConvertDbl(T,"K"); |
---|
26 | |
---|
27 | % do fermi calculations (science)... |
---|
28 | kT = 8.61734e-5 * T; |
---|
29 | Emin = Ef - 10*kT; |
---|
30 | Emax = Ef + 10*kT; |
---|
31 | |
---|
32 | E = linspace(Emin,Emax,200); |
---|
33 | f = 1.0 ./ (1.0 + exp((E - Ef)/kT)); |
---|
34 | |
---|
35 | % prepare out output section |
---|
36 | % label graphs |
---|
37 | rpLibPutString(lib,"output.curve(f12).about.label","Fermi-Dirac Factor",0); |
---|
38 | rpLibPutString(lib,"output.curve(f12).xaxis.label","Fermi-Dirac Factor",0); |
---|
39 | rpLibPutString(lib,"output.curve(f12).yaxis.label","Energy",0); |
---|
40 | rpLibPutString(lib,"output.curve(f12).yaxis.units","eV",0); |
---|
41 | |
---|
42 | % this is a slow and inefficient method of putting |
---|
43 | % large amounts of data back in to the rappture library object |
---|
44 | %for j=1:200 |
---|
45 | % rpUtilsProgress((j/200*100),'Iterating'); |
---|
46 | % putStr = sprintf('%12g %12g\n', f(j), E(j)); |
---|
47 | % % put the data into the xml file |
---|
48 | % rpLibPutString(lib,"output.curve(f12).component.xy",putStr,1); |
---|
49 | %end |
---|
50 | |
---|
51 | % a better way is to take advantage of octave's vector operations. |
---|
52 | outData = [f;E]; |
---|
53 | putStr = sprintf('%12g %12g\n', outData); |
---|
54 | rpLibPutString(lib,'output.curve(f12).component.xy',putStr,0); |
---|
55 | |
---|
56 | % signal the end of processing |
---|
57 | rpLibResult(lib); |
---|
58 | quit; |
---|