I have C++ code which solves a partial differential equation numerically. Everything works fine, but after the first 10 or so iterations, the code starts to run about 3 times as slow.
Some possibly relevant information:
- I am using the package Eigen for a lot of calculations
- I am storing a large amount of data in an array (about 100000 double precision values every iteration)
I don't know very much about how computers work or about coding in c++, and was wondering if anybody might be able to point me in the right direction so I can possibly fix this problem.
I can also post the code if it would be helpful, but it's very long.