A partial hybrid CPU/GPU solver done.

Now we have: 1. multi-partitioned ILUT (on CPU); 2. calculate
spikes (on GPU); 3. sparse sweeps (on CPU).