Add subsystem integration using regularization method like TTL or logH (Slow down will come in an other PR)

Add subsystem integration using regularization method like TTL or logH (Slow down will come in an other PR) #96