Commit Graph

148 Commits

Author SHA1 Message Date
39c05d2cab no reason to demand float32 yet 2019-07-22 00:27:48 -07:00
89976647f2 test (and fix tests) for constant non-1 dxes
Still need to look at non-constant dxes
2019-07-22 00:27:32 -07:00
7092c13088 better error messages when tests fail 2019-07-22 00:26:34 -07:00
f1fc308d25 Add JdotE test 2019-07-21 22:06:57 -07:00
7f8a326114 Loosen tolerances on tests 2019-07-21 22:06:24 -07:00
30ddeb7b73 fix typo in fdfd.vec() 2019-07-21 22:05:40 -07:00
223b202d03 More test cases 2019-07-19 00:19:47 -07:00
fb3c88a78d add test_poynting_planes 2019-07-19 00:19:32 -07:00
950e70213a Consolidate variables in test case setups 2019-07-18 00:03:32 -07:00
f858cb8bbb Fix poynting e2h test 2019-07-17 23:48:04 -07:00
2cec4fabaf Account for dxes 2019-07-17 23:47:45 -07:00
a528effd89 add some more tests 2019-07-17 00:51:28 -07:00
935b2c9a80 remove extra dt 2019-07-17 00:51:13 -07:00
79e14af4db poynting divergence doesn't use dt, and can have default dxes 2019-07-17 00:50:49 -07:00
dd4e6f294f update fdtd and add some fdtd tests 2019-07-15 01:21:12 -07:00
ecaf9fa3d0 Test code for cylindrical wg modesolver 2019-07-09 20:21:14 -07:00
099966f291 Add poynting vector and divergence 2019-07-09 20:20:49 -07:00
a8a5a69e1a Eliminate iterations over lists (assume ndarray instead of list of ndarrays) 2019-07-09 20:20:05 -07:00
557e748356 Reduce number of allocations during maxwell curls 2019-07-09 20:19:35 -07:00
9d1d8fe869 Improve wisdom management 2019-07-09 20:13:49 -07:00
8e634e35df Add experimental source types 2019-07-09 20:13:31 -07:00
4c2035c882 Add m2j() function 2019-07-09 20:13:07 -07:00
d462ae9412 unvec to (3, *shape) rather than list-of-ndarrays 2019-07-09 20:12:48 -07:00
2acbda4764 Force slices to be a tuple 2019-07-09 20:12:03 -07:00
3a5d75cde4 fix typo 2019-07-09 20:11:45 -07:00
2b3a74b737 Fix waveguide source computation for different polarities etc. 2019-07-09 20:11:32 -07:00
5dd26915fc wavenumber correction must take dx into account 2019-07-09 20:09:12 -07:00
c3f248a73c Clarify beta=wavenumber 2019-07-09 20:08:44 -07:00
001c32a2e0 Partially fix arbitrary mode phase 2019-07-09 20:08:33 -07:00
41cd94fe48 More detailed logging 2019-07-09 20:07:44 -07:00
c7d4c4a8e6 Add callback for block mode solve progress 2019-07-09 20:07:13 -07:00
jan
1f9a9949c0 Clarify memo and cleanup 2018-01-15 22:44:59 -08:00
jan
323bcf88ad Propagate mu correctly 2018-01-15 22:44:26 -08:00
jan
ee9abb77d9 Fix approx_inverse operator 2018-01-15 22:44:14 -08:00
jan
c1f65f61c1 Use pyfftw if available 2018-01-15 22:43:59 -08:00
jan
e8f836c908 Cleanup 2018-01-15 22:43:33 -08:00
jan
0e47fdd5fb randomize imaginary part of starting vector 2018-01-09 00:00:58 -08:00
jan
e02040c709 fixes and clarification 2018-01-08 23:33:22 -08:00
jan
c4cbdff751 cleanup 2018-01-08 23:28:57 -08:00
jan
4067766478 use own CG implementation 2018-01-08 17:01:44 -08:00
jan
47dd0df8bc fix operator test 2018-01-08 11:59:00 -08:00
jan
66712efd49 scipy L-BFGS silently converts to float, so view as floats when dealing with it.' 2017-12-27 01:44:45 -08:00
jan
a70687f5e3 add bloch example 2017-12-21 20:11:42 -08:00
jan
85030448c3 Use L-BFGS instead of CG, and remove rayleigh iteration refinement
scipy CG doesn't seem to converge as well as L-BFGS... worth looking
into later
2017-12-21 20:11:30 -08:00
jan
16f97d7f6b Add ability to set bounds for find_k 2017-12-18 00:13:29 -08:00
jan
f312d73503 Return real part of the gradient 2017-12-17 22:55:55 -08:00
jan
39979edc44 implement eigenvalue algorithm from Johnson paper. Could also use arpack + refinement, but that's also slow. 2017-12-17 21:33:53 -08:00
jan
4a9596921f rename search_direction to direction 2017-12-17 21:32:59 -08:00
jan
000cfabd78 switch fft, ifft 2017-12-17 21:32:29 -08:00
jan
d09eff990f Update Rayleigh quotient iteration to allow arbitrary linear operators 2017-12-17 20:51:34 -08:00