fix description on power flow in doc/theory.pdf
change interface of Network.loss() and Network.optimize()
- loss() returns detailed information like lower bound with an option
- optimize() returns only an optimal configuration
update README
- modify it according to the two proceeding commits
- how to select configurations common to multiple load profiles
- modify description about service restoration
- several small modifications