openstudio::TimeSeries ScheduleVariableInterval_Impl::timeSeries() const { unsigned numExtensibleGroups = this->numExtensibleGroups(); if (numExtensibleGroups == 0){ return TimeSeries(Date(MonthOfYear::Jan, 1), 0, Vector(), ""); } DateTimeVector dateTimes; Vector values(numExtensibleGroups); unsigned i = 0; for (const ModelExtensibleGroup& group : castVector<ModelExtensibleGroup>(extensibleGroups())) { OptionalInt month = group.getInt(0); OptionalInt day = group.getInt(1); OptionalInt hour = group.getInt(2); OptionalInt minute = group.getInt(3); OptionalDouble x = group.getDouble(4); OS_ASSERT(month); OS_ASSERT(day); OS_ASSERT(hour); OS_ASSERT(minute); OS_ASSERT(x); dateTimes.push_back(DateTime(Date(MonthOfYear(*month), *day), Time(0, *hour, *minute))); values[i] = *x; ++i; } TimeSeries result(dateTimes, values, ""); result.setOutOfRangeValue(this->outOfRangeValue()); return result; }
TEST(Filetypes, EpwFile_International_Data) { try{ path p = resourcesPath() / toPath("utilities/Filetypes/CHN_Guangdong.Shaoguan.590820_CSWD.epw"); EpwFile epwFile(p,true); EXPECT_EQ(p, epwFile.path()); EXPECT_EQ("B68C068B", epwFile.checksum()); EXPECT_EQ("Shaoguan", epwFile.city()); EXPECT_EQ("Guangdong", epwFile.stateProvinceRegion()); EXPECT_EQ("CHN", epwFile.country()); EXPECT_EQ("CSWD", epwFile.dataSource()); EXPECT_EQ("590820", epwFile.wmoNumber()); EXPECT_EQ(24.68, epwFile.latitude()); EXPECT_EQ(113.6, epwFile.longitude()); EXPECT_EQ(8, epwFile.timeZone()); EXPECT_EQ(61, epwFile.elevation()); EXPECT_EQ(Time(0,1,0,0), epwFile.timeStep()); EXPECT_EQ(DayOfWeek(DayOfWeek::Sunday), epwFile.startDayOfWeek()); EXPECT_EQ(Date(MonthOfYear::Jan, 1), epwFile.startDate()); EXPECT_EQ(Date(MonthOfYear::Dec, 31), epwFile.endDate()); // Up to here, everything should be the same as the first test. Now ask for the data std::vector<EpwDataPoint> data = epwFile.data(); EXPECT_EQ(8760,data.size()); // The last data point check EXPECT_EQ(14.7,data[8759].dryBulbTemperature().get()); EXPECT_EQ(101100,data[8759].atmosphericStationPressure().get()); // Try out the alternate access functions, dew point temperature should be -1C EXPECT_EQ(11.7,data[8759].fieldByName("Dew Point Temperature").get()); EXPECT_EQ(11.7,data[8759].field(EpwDataField("Dew Point Temperature")).get()); // The last data point should not have a liquid precipitation depth EXPECT_FALSE(data[8759].fieldByName("Liquid Precipitation Depth")); // Get a time series boost::optional<openstudio::TimeSeries> series = epwFile.getTimeSeries("Wind Speed"); ASSERT_TRUE(series); ASSERT_EQ(8760,series->values().size()); DateTimeVector seriesTimes = series->dateTimes(); ASSERT_EQ(8760,seriesTimes.size()); // Check the times in the data and the time series DateTime current(Date(1,1,1999),Time(0,1)); // Use 1999 to avoid leap years Time delta(0,1); for(unsigned i=0;i<8760;i++) { // This is a lot more complicated that it probably should be to avoid the year being a problem DateTime datatime = data[i].dateTime(); EXPECT_EQ(datatime.date().monthOfYear(), current.date().monthOfYear()); EXPECT_EQ(datatime.date().dayOfMonth(), current.date().dayOfMonth()); EXPECT_EQ(datatime.time().hours(), current.time().hours()); EXPECT_EQ(datatime.time().minutes(), current.time().minutes()); DateTime seriestime = seriesTimes[i]; EXPECT_EQ(seriestime.date().monthOfYear(), current.date().monthOfYear()); EXPECT_EQ(seriestime.date().dayOfMonth(), current.date().dayOfMonth()); EXPECT_EQ(seriestime.time().hours(), current.time().hours()); EXPECT_EQ(seriestime.time().minutes(), current.time().minutes()); current += delta; } // No need to redo the original tests here since the data should have been loaded in the constructor }catch(...){ ASSERT_TRUE(false); } }
TEST(Filetypes, EpwFile_IWEC_Data) { try { path p = resourcesPath() / toPath("utilities/Filetypes/TUN_Tunis.607150_IWEC.epw"); EpwFile epwFile(p, true); EXPECT_EQ(p, epwFile.path()); EXPECT_EQ("FEAB878E", epwFile.checksum()); EXPECT_EQ("TUNIS", epwFile.city()); EXPECT_EQ("-", epwFile.stateProvinceRegion()); EXPECT_EQ("TUN", epwFile.country()); EXPECT_EQ("IWEC Data", epwFile.dataSource()); EXPECT_EQ("607150", epwFile.wmoNumber()); EXPECT_EQ(36.83, epwFile.latitude()); EXPECT_EQ(10.23, epwFile.longitude()); EXPECT_EQ(1, epwFile.timeZone()); EXPECT_EQ(4, epwFile.elevation()); EXPECT_EQ(Time(0, 1, 0, 0), epwFile.timeStep()); EXPECT_EQ(DayOfWeek(DayOfWeek::Sunday), epwFile.startDayOfWeek()); EXPECT_EQ(Date(MonthOfYear::Jan, 1), epwFile.startDate()); EXPECT_EQ(Date(MonthOfYear::Dec, 31), epwFile.endDate()); // Up to here, everything should be the same as the first test. Now ask for the data std::vector<EpwDataPoint> data = epwFile.data(); EXPECT_EQ(8760, data.size()); // The last data point check EXPECT_EQ(11.3, data[8759].dryBulbTemperature().get()); EXPECT_EQ(102400, data[8759].atmosphericStationPressure().get()); // Try out the alternate access functions EXPECT_EQ(9.8, data[8759].getFieldByName("Dew Point Temperature").get()); EXPECT_EQ(9.8, data[8759].getField(EpwDataField("Dew Point Temperature")).get()); // Get a time series boost::optional<openstudio::TimeSeries> series = epwFile.getTimeSeries("Wind Speed"); ASSERT_TRUE(series); ASSERT_EQ(8760, series->values().size()); DateTimeVector seriesTimes = series->dateTimes(); ASSERT_EQ(8760, seriesTimes.size()); // Check the times in the data and the time series DateTime current(Date(1, 1, 1999), Time(0, 1)); // Use 1999 to avoid leap years Time delta(0, 1); for (unsigned i = 0; i<8760; i++) { // This is a lot more complicated that it probably should be to avoid the year being a problem DateTime datatime = data[i].dateTime(); EXPECT_EQ(datatime.date().monthOfYear(), current.date().monthOfYear()); EXPECT_EQ(datatime.date().dayOfMonth(), current.date().dayOfMonth()); EXPECT_EQ(datatime.time().hours(), current.time().hours()); EXPECT_EQ(datatime.time().minutes(), current.time().minutes()); DateTime seriestime = seriesTimes[i]; EXPECT_EQ(seriestime.date().monthOfYear(), current.date().monthOfYear()); EXPECT_EQ(seriestime.date().dayOfMonth(), current.date().dayOfMonth()); EXPECT_EQ(seriestime.time().hours(), current.time().hours()); EXPECT_EQ(seriestime.time().minutes(), current.time().minutes()); current += delta; } // No need to redo the original tests here since the data should have been loaded in the constructor } catch (...) { ASSERT_TRUE(false); } }
TEST_F(DataFixture,TimeSeries_AddSubtractSameTimePeriod) { std::string units = "W"; Date startDate(Date(MonthOfYear(MonthOfYear::Feb),21)); DateTime startDateTime(startDate, Time(0,1,0,0)); // interval Time interval = Time(0,1,0,0); Vector intervalValues(3); intervalValues(0) = 0; intervalValues(1) = 1; intervalValues(2) = 2; TimeSeries intervalTimeSeries(startDateTime, interval, intervalValues, units); ASSERT_TRUE(!intervalTimeSeries.values().empty()); // detailed DateTimeVector dateTimes; dateTimes.push_back(startDateTime + Time(0,0,0,0)); dateTimes.push_back(startDateTime + Time(0,0,30,0)); dateTimes.push_back(startDateTime + Time(0,1,0,0)); dateTimes.push_back(startDateTime + Time(0,1,30,0)); dateTimes.push_back(startDateTime + Time(0,2,0,0)); Vector detailedValues(5); detailedValues(0) = 0.0; // 1:00 detailedValues(1) = 0.5; // 1:30 detailedValues(2) = 1.0; // 2:00 detailedValues(3) = 1.5; // 2:30 detailedValues(4) = 2.0; // 3:00 TimeSeries detailedTimeSeries(dateTimes, detailedValues, units); ASSERT_TRUE(!detailedTimeSeries.values().empty()); // sum and difference TimeSeries sum = intervalTimeSeries + detailedTimeSeries; TimeSeries diff1 = intervalTimeSeries - detailedTimeSeries; TimeSeries diff2 = detailedTimeSeries - intervalTimeSeries; ASSERT_TRUE(!sum.values().empty()); ASSERT_TRUE(!diff1.values().empty()); ASSERT_TRUE(!diff2.values().empty()); // EXPECT_EQ((unsigned)5, sum.dateTimes().size()); // EXPECT_EQ((unsigned)5, diff1.dateTimes().size()); // EXPECT_EQ((unsigned)5, diff2.dateTimes().size()); EXPECT_EQ((unsigned)5, sum.daysFromFirstReport().size()); EXPECT_EQ((unsigned)5, diff1.daysFromFirstReport().size()); EXPECT_EQ((unsigned)5, diff2.daysFromFirstReport().size()); // EXPECT_EQ(startDateTime, sum.dateTimes().front()); // EXPECT_EQ(startDateTime, diff1.dateTimes().front()); // EXPECT_EQ(startDateTime, diff2.dateTimes().front()); EXPECT_EQ(startDateTime, sum.firstReportDateTime()); EXPECT_EQ(startDateTime, diff1.firstReportDateTime()); EXPECT_EQ(startDateTime, diff2.firstReportDateTime()); DateTime endDateTime = startDateTime + Time(0,2,0,0); // EXPECT_EQ(endDateTime, sum.dateTimes().back()); // EXPECT_EQ(endDateTime, diff1.dateTimes().back()); // EXPECT_EQ(endDateTime, diff2.dateTimes().back()); EXPECT_EQ(endDateTime, sum.firstReportDateTime() + Time(sum.daysFromFirstReport(sum.daysFromFirstReport().size()-1))); EXPECT_EQ(endDateTime, diff1.firstReportDateTime() + Time(diff1.daysFromFirstReport(diff1.daysFromFirstReport().size()-1))); EXPECT_EQ(endDateTime, diff2.firstReportDateTime() + Time(diff2.daysFromFirstReport(diff2.daysFromFirstReport().size()-1))); // 1:00 EXPECT_EQ(0, sum.value(Time(0,0,0,0))); EXPECT_EQ(0, diff1.value(Time(0,0,0,0))); EXPECT_EQ(0, diff2.value(Time(0,0,0,0))); // 1:30 EXPECT_EQ(1.5, sum.value(Time(0,0,30,0))); EXPECT_EQ(0.5, diff1.value(Time(0,0,30,0))); EXPECT_EQ(-0.5, diff2.value(Time(0,0,30,0))); // 2:00 EXPECT_EQ(2, sum.value(Time(0,1,0,0))); EXPECT_EQ(0.0, diff1.value(Time(0,1,0,0))); EXPECT_EQ(0.0, diff2.value(Time(0,1,0,0))); // 2:30 EXPECT_EQ(3.5, sum.value(Time(0,1,30,0))); EXPECT_EQ(0.5, diff1.value(Time(0,1,30,0))); EXPECT_EQ(-0.5, diff2.value(Time(0,1,30,0))); // Test helper function for summing a vector. TimeSeriesVector sumAndDiffs; sumAndDiffs.push_back(sum); sumAndDiffs.push_back(diff1); sumAndDiffs.push_back(diff2); TimeSeries ans = openstudio::sum(sumAndDiffs); EXPECT_FALSE(ans.values().empty()); // 1:00 EXPECT_DOUBLE_EQ(0, ans.value(Time(0,0,0,0))); // 1:30 EXPECT_DOUBLE_EQ(1.5, ans.value(Time(0,0,30,0))); // 2:00 EXPECT_DOUBLE_EQ(2.0, ans.value(Time(0,1,0,0))); // 2:30 EXPECT_DOUBLE_EQ(3.5, ans.value(Time(0,1,30,0))); // Test multiplication and division with a scalar sumAndDiffs.push_back(sum/2.0); sumAndDiffs.push_back(3.0*diff1); ans = openstudio::sum(sumAndDiffs); EXPECT_FALSE(ans.values().empty()); // 1:00 EXPECT_DOUBLE_EQ(0, ans.value(Time(0,0,0,0))); // 1:30 EXPECT_DOUBLE_EQ(3.75, ans.value(Time(0,0,30,0))); // 2:00 EXPECT_DOUBLE_EQ(3.0, ans.value(Time(0,1,0,0))); // 2:30 EXPECT_DOUBLE_EQ(6.75, ans.value(Time(0,1,30,0))); }
TEST_F(DataFixture,TimeSeries_DetailedConstructor) { std::string units = "W"; Date startDate(Date(MonthOfYear(MonthOfYear::Feb),21)); DateTime startDateTime(startDate, Time(0,1,0,0)); // fill vector with 3 hours of data Vector values(3); DateTimeVector dateTimes; for (unsigned i = 0; i < 3; ++i){ values(i) = i; dateTimes.push_back(startDateTime + Time(0,i,0,0)); } unsigned numValues = values.size(); // create detailed timeSeries TimeSeries timeSeries(dateTimes, values, units); ASSERT_TRUE(!timeSeries.values().empty()); // check interval EXPECT_FALSE(timeSeries.intervalLength()); // check start date and time // DateTime firstDateTime = timeSeries.dateTimes().front(); DateTime firstDateTime = timeSeries.firstReportDateTime(); EXPECT_EQ(DateTime(Date(MonthOfYear(MonthOfYear::Feb), 21), Time(0,1,0,0)), firstDateTime); // check end date and time // DateTime endDateTime = timeSeries.dateTimes().back(); DateTime endDateTime = timeSeries.firstReportDateTime() + Time(timeSeries.daysFromFirstReport(timeSeries.daysFromFirstReport().size()-1)); EXPECT_EQ(DateTime(Date(MonthOfYear(MonthOfYear::Feb), 21), Time(0,3,0,0)), endDateTime); // check out of range timeSeries.setOutOfRangeValue(-99); EXPECT_EQ(-99, timeSeries.value(startDateTime + Time(0,0,-61,0))); // out of range EXPECT_EQ(-99, timeSeries.value(startDateTime + Time(0,0,-60,0))); // out of range EXPECT_EQ(-99, timeSeries.value(startDateTime + Time(0,0,-59,0))); // out of range EXPECT_EQ(-99, timeSeries.value(startDateTime + Time(0,0,-1,0))); // out of range EXPECT_EQ(0, timeSeries.value(startDateTime)); // in range EXPECT_EQ(2, timeSeries.value(endDateTime)); // in range EXPECT_EQ(-99, timeSeries.value(endDateTime + Time(0,1,0,0))); // out of range // check values for (unsigned i = 0; i < numValues; ++i){ double numPeriods = (double)i; Time interval(0,1,0,0); double lastPeriodEnd = (numPeriods-1.0)*interval.totalDays(); double periodBegin = (numPeriods-0.99)*interval.totalDays(); double periodMiddle = (numPeriods-0.5)*interval.totalDays(); double periodEnd = (numPeriods)*interval.totalDays(); if (i > 0){ EXPECT_EQ(numPeriods-1, timeSeries.value(lastPeriodEnd)); EXPECT_EQ(numPeriods, timeSeries.value(periodBegin)); EXPECT_EQ(numPeriods, timeSeries.value(periodMiddle)); } EXPECT_EQ(numPeriods, timeSeries.value(periodEnd)); } }
TEST(Filetypes, EpwFile_Data) { try{ path p = resourcesPath() / toPath("runmanager/USA_CO_Golden-NREL.724666_TMY3.epw"); EpwFile epwFile(p); EXPECT_EQ(p, epwFile.path()); EXPECT_EQ("F188656D", epwFile.checksum()); EXPECT_EQ("Denver Centennial Golden Nr", epwFile.city()); EXPECT_EQ("CO", epwFile.stateProvinceRegion()); EXPECT_EQ("USA", epwFile.country()); EXPECT_EQ("TMY3", epwFile.dataSource()); EXPECT_EQ("724666", epwFile.wmoNumber()); EXPECT_EQ(39.74, epwFile.latitude()); EXPECT_EQ(-105.18, epwFile.longitude()); EXPECT_EQ(-7, epwFile.timeZone()); EXPECT_EQ(1829, epwFile.elevation()); EXPECT_EQ(Time(0,1,0,0), epwFile.timeStep()); EXPECT_EQ(DayOfWeek(DayOfWeek::Sunday), epwFile.startDayOfWeek()); EXPECT_EQ(Date(MonthOfYear::Jan, 1), epwFile.startDate()); EXPECT_EQ(Date(MonthOfYear::Dec, 31), epwFile.endDate()); // Up to here, everything should be the same as the first test. Now ask for the data std::vector<EpwDataPoint> data = epwFile.data(); EXPECT_EQ(8760,data.size()); // The last data point should be on 12/31/1996, with a dry bulb temp of 4C and presure 81100 EXPECT_EQ(4.0,data[8759].dryBulbTemperature().get()); EXPECT_EQ(81100,data[8759].atmosphericStationPressure().get()); // Try out the alternate access functions, dew point temperature should be -1C EXPECT_EQ(-1.0,data[8759].fieldByName("Dew Point Temperature").get()); EXPECT_EQ(-1.0,data[8759].field(EpwDataField("Dew Point Temperature")).get()); // The last data point should not have a liquid precipitation depth EXPECT_FALSE(data[8759].fieldByName("Liquid Precipitation Depth")); // Get a time series boost::optional<openstudio::TimeSeries> series = epwFile.getTimeSeries("Wind Speed"); ASSERT_TRUE(series); ASSERT_EQ(8760,series->values().size()); DateTimeVector seriesTimes = series->dateTimes(); ASSERT_EQ(8760,seriesTimes.size()); // Check the times in the data and the time series DateTime current(Date(1,1,1999),Time(0,1)); // Use 1999 to avoid leap years Time delta(0,1); for(unsigned i=0;i<8760;i++) { // This is a lot more complicated that it probably should be to avoid the year being a problem DateTime datatime = data[i].dateTime(); EXPECT_EQ(datatime.date().monthOfYear(), current.date().monthOfYear()); EXPECT_EQ(datatime.date().dayOfMonth(), current.date().dayOfMonth()); EXPECT_EQ(datatime.time().hours(), current.time().hours()); EXPECT_EQ(datatime.time().minutes(), current.time().minutes()); DateTime seriestime = seriesTimes[i]; EXPECT_EQ(seriestime.date().monthOfYear(), current.date().monthOfYear()); EXPECT_EQ(seriestime.date().dayOfMonth(), current.date().dayOfMonth()); EXPECT_EQ(seriestime.time().hours(), current.time().hours()); EXPECT_EQ(seriestime.time().minutes(), current.time().minutes()); current += delta; } // We should redo the original tests because we have reparsed the entire file EXPECT_EQ(p, epwFile.path()); EXPECT_EQ("F188656D", epwFile.checksum()); EXPECT_EQ("Denver Centennial Golden Nr", epwFile.city()); EXPECT_EQ("CO", epwFile.stateProvinceRegion()); EXPECT_EQ("USA", epwFile.country()); EXPECT_EQ("TMY3", epwFile.dataSource()); EXPECT_EQ("724666", epwFile.wmoNumber()); EXPECT_EQ(39.74, epwFile.latitude()); EXPECT_EQ(-105.18, epwFile.longitude()); EXPECT_EQ(-7, epwFile.timeZone()); EXPECT_EQ(1829, epwFile.elevation()); EXPECT_EQ(Time(0,1,0,0), epwFile.timeStep()); EXPECT_EQ(DayOfWeek(DayOfWeek::Sunday), epwFile.startDayOfWeek()); EXPECT_EQ(Date(MonthOfYear::Jan, 1), epwFile.startDate()); EXPECT_EQ(Date(MonthOfYear::Dec, 31), epwFile.endDate()); }catch(...){ ASSERT_TRUE(false); } }
TEST_F(DataFixture,TimeSeries_DetailedConstructor_WrapAroundDates) { std::string units = "W"; Date startDate(MonthOfYear(MonthOfYear::Jan),1); // fill vector with 365 days of data Vector values(365); DateTimeVector dateTimes; for (unsigned i = 0; i < 365; ++i){ values(i) = i; dateTimes.push_back(startDate + Time(i+1,0,0,0)); } unsigned numValues = values.size(); ASSERT_EQ(365u, numValues); ASSERT_EQ(365u, dateTimes.size()); EXPECT_EQ(DateTime(Date(MonthOfYear(MonthOfYear::Jan), 1), Time(1,0,0,0)), dateTimes[0].date()); EXPECT_EQ(DateTime(Date(MonthOfYear(MonthOfYear::Dec), 31), Time(1,0,0,0)), dateTimes[364].date()); // create detailed timeSeries TimeSeries timeSeries(dateTimes, values, units); ASSERT_EQ(365u, timeSeries.values().size()); // check interval EXPECT_FALSE(timeSeries.intervalLength()); // check start date and time DateTime firstDateTime = timeSeries.firstReportDateTime(); EXPECT_EQ(DateTime(Date(MonthOfYear(MonthOfYear::Jan), 1), Time(1,0,0,0)), firstDateTime); // check values for (unsigned i = 0; i < numValues; ++i){ DateTime dateTime = dateTimes[i]; EXPECT_EQ(i, timeSeries.value(dateTime)) << dateTime; } Vector test = timeSeries.values(DateTime(Date(MonthOfYear(MonthOfYear::Apr), 1), Time(1,0,0,0)), DateTime(Date(MonthOfYear(MonthOfYear::Apr), 30), Time(1,0,0,0))); ASSERT_EQ(30, test.size()); for (unsigned i = 0; i < 30; ++i){ EXPECT_EQ(i + 90, test[i]); } // now rearrange date times to make a wrap around year Vector wrappedValues(365); DateTimeVector wrappedDateTimes; unsigned j = 0; for (unsigned i = 100; i < 365; ++i, ++j){ wrappedValues[j] = values[i]; wrappedDateTimes.push_back(dateTimes[i]); } for (unsigned i = 0; i < 100; ++i, ++j){ wrappedValues[j] = values[i]; wrappedDateTimes.push_back(dateTimes[i]); } ASSERT_EQ(365u, wrappedValues.size()); ASSERT_EQ(365u, wrappedDateTimes.size()); EXPECT_EQ(100, wrappedValues[0]); EXPECT_EQ(99, wrappedValues[364]); EXPECT_EQ(DateTime(Date(MonthOfYear(MonthOfYear::Apr), 11), Time(1,0,0,0)), wrappedDateTimes[0].date()); EXPECT_EQ(DateTime(Date(MonthOfYear(MonthOfYear::Apr), 10), Time(1,0,0,0)), wrappedDateTimes[364].date()); // create detailed timeSeries TimeSeries wrappedTimeSeries(wrappedDateTimes, wrappedValues, units); ASSERT_EQ(365u, wrappedTimeSeries.values().size()); // check interval EXPECT_FALSE(wrappedTimeSeries.intervalLength()); // check start date and time firstDateTime = wrappedTimeSeries.firstReportDateTime(); EXPECT_EQ(DateTime(Date(MonthOfYear(MonthOfYear::Apr), 11), Time(1,0,0,0)), firstDateTime); // check values for (unsigned i = 0; i < numValues; ++i){ DateTime dateTime = dateTimes[i]; EXPECT_EQ(i, wrappedTimeSeries.value(dateTime)) << dateTime; } test = timeSeries.values(DateTime(Date(MonthOfYear(MonthOfYear::Apr), 1), Time(1,0,0,0)), DateTime(Date(MonthOfYear(MonthOfYear::Apr), 30), Time(1,0,0,0))); ASSERT_EQ(30, test.size()); for (unsigned i = 0; i < 30; ++i){ EXPECT_EQ(i + 90, test[i]); } }