def default_residential_outputs_1(request, gsod_722880_2012_2014_weather_source): elec_model_params, gas_model_params, \ elec_annualized_usage, gas_annualized_usage, \ elec_gross_savings, gas_gross_savings, \ elec_rmse, gas_rmse, elec_r_squared, gas_r_squared, \ temp_unit, cdd_tmy, hdd_tmy, total_cdd, total_hdd \ = request.param period = Period(datetime(2012,1,1,tzinfo=pytz.utc), datetime(2014,12,31,tzinfo=pytz.utc)) retrofit_start_date = datetime(2013,6,1,tzinfo=pytz.utc) retrofit_completion_date = datetime(2013,8,1,tzinfo=pytz.utc) datetimes = generate_monthly_billing_datetimes(period,randint(30,31)) # generate electricity consumption elec_model = AverageDailyTemperatureSensitivityModel(cooling=True,heating=True) elec_params = { "base_daily_consumption": elec_model_params[0], "heating_balance_temperature": elec_model_params[1], "heating_slope": elec_model_params[2], "cooling_balance_temperature": elec_model_params[3], "cooling_slope": elec_model_params[4], } elec_gen = MonthlyBillingConsumptionGenerator("electricity", "kWh", temp_unit, elec_model, elec_params) elec_consumption_data = elec_gen.generate(gsod_722880_2012_2014_weather_source, datetimes) elec_consumption_kWh_per_day, elec_consumption_n_days = \ elec_consumption_data.average_daily_consumptions() elec_params = elec_model.param_type(elec_params) # generate natural_gas consumption gas_model = AverageDailyTemperatureSensitivityModel(cooling=False,heating=True) gas_params = { "base_daily_consumption": gas_model_params[0], "heating_balance_temperature": gas_model_params[1], "heating_slope": gas_model_params[2], } gas_gen = MonthlyBillingConsumptionGenerator("natural_gas", "therm", temp_unit, gas_model, gas_params) gas_consumption_data = gas_gen.generate(gsod_722880_2012_2014_weather_source, datetimes) gas_consumption_kWh_per_day, gas_consumption_n_days = \ gas_consumption_data.average_daily_consumptions() gas_params = gas_model.param_type(gas_params) fixture = elec_consumption_data, gas_consumption_data, \ elec_params, gas_params, \ elec_annualized_usage, gas_annualized_usage, \ elec_gross_savings, gas_gross_savings, \ elec_rmse, gas_rmse, \ elec_r_squared, gas_r_squared, \ elec_consumption_kWh_per_day, gas_consumption_kWh_per_day, \ elec_consumption_n_days, gas_consumption_n_days, \ temp_unit, retrofit_start_date, retrofit_completion_date, \ cdd_tmy, hdd_tmy, total_cdd, total_hdd return fixture
def test_model_weather_input_not_np_array(): model = AverageDailyTemperatureSensitivityModel(heating=False, cooling=False) params = model.param_type([1]) observed_temps = [[70], [65, 60]] usages = model.transform(observed_temps, params) assert_allclose(usages, [1, 1], rtol=1e-2, atol=1e-2)
def test_average_daily_baseload_heating_cooling_consumption_model(): initial_params = { "base_daily_consumption": 0, "heating_slope": 0, "cooling_slope": 0, "heating_balance_temperature": 55, "cooling_balance_temperature": 57, } param_bounds = { "base_daily_consumption": [0, 100], "heating_slope": [0, 100], "cooling_slope": [0, 100], "heating_balance_temperature": [50, 60], "cooling_balance_temperature": [52, 72], } model = AverageDailyTemperatureSensitivityModel( cooling=True, heating=True, initial_params=initial_params, param_bounds=param_bounds) params = model.param_type([1, 60, 1, 65, 1]) observed_temps = np.array([[i] for i in range(50, 70)]) usages = model.transform(observed_temps, params) assert_allclose(usages[8:18], [3, 2, 1, 1, 1, 1, 1, 1, 2, 3], rtol=1e-2, atol=1e-2) opt_params = model.fit(observed_temps, usages) assert_allclose(params.to_list(), opt_params.to_list(), rtol=1e-2, atol=1e-2)
def consumption_generator_2(request): model_params, fuel_type, consumption_unit_name, temperature_unit_name = \ request.param model = AverageDailyTemperatureSensitivityModel(cooling=True, heating=True) params = { "base_daily_consumption": model_params[0], "heating_slope": model_params[1], "heating_balance_temperature": model_params[2], "cooling_slope": model_params[3], "cooling_balance_temperature": model_params[4] } generator = MonthlyBillingConsumptionGenerator(fuel_type, consumption_unit_name, temperature_unit_name, model, params) params = model.param_type(params) return generator, params
def test_TemperatureSensitivityModel_with_cooling(): initial_params = { "base_daily_consumption": 0, "cooling_slope": 0, "cooling_balance_temperature": 57, } param_bounds = { "base_daily_consumption": [0,100], "cooling_slope": [0,100], "cooling_balance_temperature": [52,72], } model = AverageDailyTemperatureSensitivityModel(heating=False,cooling=True,initial_params=initial_params,param_bounds=param_bounds) params = model.param_type([1,60,1]) observed_temps = np.array([[i] for i in range(50,70)]) usages = model.transform(observed_temps, params) assert_allclose(usages[8:13],[1,1,1,2,3]) opt_params = model.fit(observed_temps, usages) assert_allclose(params.to_list(), opt_params.to_list(), rtol=1e-2, atol=1e-2)
def test_model_weather_input_not_np_array(): model = AverageDailyTemperatureSensitivityModel(heating=False,cooling=False) params = model.param_type([1]) observed_temps = [[70],[65,60]] usages = model.transform(observed_temps, params) assert_allclose(usages, [1,1], rtol=1e-2, atol=1e-2)