periodic bc implemented.
before editing restriction.
outputting different matrix and vector for differnt level.
global constraint good.
finite difference matrix and load vector have to be initilized before use.
restriction verified to work.
not working...
seems working!
now working. trying to parallelize as much as possible.
i think i parallelize it as much as possible and seems working properly.
parallels assembly of sparse fd matrix almost complete except global constraint.
sparse matrix assembly verified.
sparse jacobi method implemented.
merge sort implemented.
works great on stampede too
main function modified.
before merging