JWST computation is now double precision, changed type of variables in JWST_data_processor to float64, complex128

Compared result to matlab result:

After 1 iteration:

Maximum norm of difference:
2.25868816441e-12
Frobenius norm of difference:
1.9427347826e-11

After 500 iterations:

Maximum norm of difference:
8.35619946335e-05
Frobenius norm of difference:
0.00514996012436

So JWST seems to be working now...
Also epsilon in magproj was changed back to 1e-30 (no overflow with double precision)