Multi-point moment matching reduction of distributed thermal networks