add test: masked data should not affect fit
ignore masked values in error function
note on fork, speed