A numerical algorithm for carrying out "Whole building hourly analyses" on PCM containing buildings