removed seasonality design matrix, using Lagrange multipliers instead
finally got seasonality working correctly