fixed matrix operation, next step=accelerated parallelism