Fast Causal Discovery by Approximate Kernel-based Generalized Score Functions with Linear Computational Complexity