Is the new bfast option implemented in 1D?

[andrew-ueltd]

I've modified the refl-angular.py standard python example to use the bfast option (see enclosed zipped .py file), but the output is independent of angle of incidence. I've defined k_point=0 in line with previous discussions since the field is the one with the exp(ik.r) term factored out. Is the bfast method properly implemented in 1D?
bfast_1D_refl-angular.zip

[andrew-ueltd]

I should also say that, based on the refl-angular.py example, if I set the parameter dimensions=3 (the enclosed file has dimensions=1) the execution fails with simulation fields NaN or Inf. Any help appreciated.