A fully-antisymmetrized random phase approximation calculation employing the continued fraction technique is performed to study nuclear matter response functions with the finite range Gogny force. The most commonly used parameter sets of this force, as well as some recent generalizations that include the tensor terms are considered and the corresponding response functions are shown. The calculations are performed at the first and second order in the continued fraction expansion and the explicit expressions for the second order tensor contributions are given. Comparison between first and second order continued fraction expansion results are provided. The differences between the responses obtained at the two orders turn to be more pronounced for the forces including tensor terms than for the standard Gogny ones. In the vector channels the responses calculated with Gogny forces including tensor terms are characterized by a large heterogeneity, reflecting the different choices for the tensor part of the interaction. For sake of comparison the response functions obtained considering a G-matrix based nuclear interaction are also shown. As first application of the present calculation, the possible existence of spurious finite-size instabilities of the Gogny forces with or without tensor terms has been investigated. The positive conclusion is that all the Gogny forces, but the GT2 one, are free of spurious finite-size instabilities. In perspective, the tool developed in the present paper can be inserted in the fitting procedure to construct new Gogny-type forces.