initial version

Ignore-this: 88f46035a701c216435a2c4208fbbd83
The changes I have: 1) combine myinterp and myinterpslab into myinterp; use Mz to denote the thickness in z direction while Mzslab=0/1 to indicate whether it's the slab case;

2)add addmuabsorption, add ignore_zero_entries options, modify MOperator (introducing an arrary)

darcs-hash:20111210232740-88e91-bfcd50754621677a3206e5e9a97be24076641777.gz
add command line options, merge three functions to one MLDOS

Ignore-this: 5d3bec0bb8e86be73cc47cef0ba1111b
In this version, 1)I change the previous passing parameters to main by taking them from commandline options. In this way, I do not need to worry about the orders. 2) Merge AnalyzeStructure, FindOpt, and ComputeQ into one single file MLDOS. Hope it will be easy to use.

darcs-hash:20111212195222-88e91-b8261a6a487dada7405f922c00bb294f3e87f591.gz
add MLDOS_sample to illustrate usage of MLDOS

Ignore-this: 2e2bca23309e26c92cbced15d59fb857

darcs-hash:20111212195732-88e91-68dad7b8efec4136dd7ee5707cd58b6a8dd76f02.gz
use VecToArray and ArrayToVec in ResonantSolver.c

Ignore-this: 5e7dac8515bff7714a4babbae642e2bc

darcs-hash:20111212202001-88e91-66dd686f0a74e9e9ec3562ea3e66908f2d660fd6.gz
bug fix for Mxyz

Ignore-this: 48f6d6c55606c6b4d560d2e02b422882
Now Mzslab = 1 for slab case and Mz=0 for nonslab case.

darcs-hash:20111217235234-88e91-5a5db042b67e83cc68b28268ea8bcafba48cf6b0.gz
create rountines for mutiple RHS and fix mldos for full domian ( without symmetry)

Ignore-this: 759fa5d3c493b1f158f102ac1633173
main changes: 1) modification on the current rounites so that it automatically adjust the position of currents based on LowerPML. 2) create routine for muliltple RHS currents from single currents. (only copied, not started yet).

darcs-hash:20111226201036-88e91-788fa92c092b9fa4591ae15fdb5589dc897d9b5f.gz
implemented multirhs rountine and API interface change for SourceSingleSetX[YZ]

Ignore-this: ed54005734c8adb759234018ed270e7f

1) I implemented the mldosrhs rountine to solve ldos at different postions.
2) API interface to SourceSingleSetX[YZ] are changed. passing *J-> passing J.
3) In mldos.c, I add cx, cy and cz to SourceSingleSet in order to indicate the position of currents ( instead of 0,0,0 used before)

darcs-hash:20120116213140-88e91-f2dbcc8e5afd7fab08654eeb2a1dbfd1fd3df7fd.gz
create MLDOSPOLXY files

Ignore-this: b99f9217954ad0e503930aed71871f1c

Create MLDOSPOLXY files to compute the ldos at single source with two polarizations. To be implemented.

darcs-hash:20120116232648-88e91-b67184349e99d5754095eca8978812701f7d6e15.gz
implemented Polarization in XY routines

Ignore-this: 617e8237b0737c4b88e1e23ab45937df

darcs-hash:20120120035345-88e91-876219a93aedfce7d1dc48e5c8678e754384192d.gz
output weight when job is analyzing

Ignore-this: 78463347f01210e493d9b4af1ace4504

darcs-hash:20120122022108-88e91-75db575134104f49ce9b83dd824a95a6e7163929.gz
allow current in y direction in mldos.c

Ignore-this: 445fda16557a78c5107a95292644627b

darcs-hash:20120123170221-88e91-48985efa04cd443597a6925ad6235b6aecac3a38.gz
change "Job" to global variable

Ignore-this: 98a375ab102250f9a053fe1558e0a2dc

darcs-hash:20120125171118-88e91-a893cf55e41355ef0f9e25d86854e095e2a69f28.gz
add MLDOSSolar (initial)

Ignore-this: e05d4124a67203cebdf7810760f7cf81

darcs-hash:20120125174119-88e91-779e28ae6606a1d48047435ca3c1ca378ca62db7.gz
implemented MLDOSSolar

Ignore-this: ba2006e91528d85f6db5e792a442fff9

darcs-hash:20120129021246-88e91-cb6d051296fcfeda66d78b1d24d351510e264f54.gz
fixed epsilon=epshi at the center

Ignore-this: f5a5a7a71733c5f6a5a34a06d7d9e26c

Following the popular definition of Purcell factor, I fixed the dielectric constant at the center.

darcs-hash:20120214032852-88e91-56873a2f807f82c1a54bf5486b218a2bb3c1bdcc.gz
relocate the codes for how to fix the center dielectric

Ignore-this: b3c668e447aa77a956e70fb92cc45baa

Instead of fix the dielectric at center in the main body, I set it only for Job=1 style.

darcs-hash:20120216163002-88e91-b92985862bc33a479b4079f36daa252e140a2696.gz
add ResonatorSolverSolar.c

Ignore-this: f8af4e0a6395f6f6b837f8aa1546a8cf
Actually, I have this file long ago, but forgot to add to darcs until I need to modify it to handle E*eps*J case.

darcs-hash:20120218234420-88e91-54a0f763e40afd3fe4c24bc3caf598ee5c2ca64a.gz
implement the new derivative in ResonatorSolverSolar.c

Ignore-this: a37d6a8e3172e7ab882bc345d945a158
I implmented the -E*J*hxyz*epscenter derivative and tested it.

darcs-hash:20120219223405-88e91-434b9d5a2c3df46c97c0d65a20dcc937952c67e3.gz
change prop def in mldossolar.c

Ignore-this: d139a5c58e9c8e1fa22f39cca8ff139f
change omega=omega*prop instead of previous omega/prop; also use tab to format the code

darcs-hash:20120219230054-88e91-e098441dae725e4b114ab4388bf96748765573a7.gz
create MoperatorGeneralBloch.c to be implemented

Ignore-this: ea748bebdc670e6abf79c50f62772c

In order to incorporate Bloch Period Boundary Condition, I create this file from MoperatorGeneral; To be implemented

darcs-hash:20120222031804-88e91-2eb54bb94603787fc473d5db02753ac10eb7df87.gz