speedup bandpass filtering along multiple channel eeg signal

I have 3 hours long 128 channel dataset sampled at 30Khz 128 x 20000000. I want to filter each channel signal in a specific frequency range say 0-12 Hz. How can I make this computation fast?