Quite a number of changes:
- templates to provide operations for fixed matrix sizes
- sorting of nets to increase convergence 
- successive over relaxation - for popeye a parameter of 1.05 works best.
- The code commented out may provide benefits on different architectures like true
  vector architectures.
  
Popeye run peaks at 1800% (up 20%) and pong has a slight performance increase as well.
Disable printf logging of stats