So I solved this using this solution:
Compare the intersect plane with the circle plane.
If they are parallel (i.e. absolute value of dot product of normals equals 1) there are zero solutions, unless they are the same plane, in which case there are infinite solutions.
Otherwise, project* the plane in the circle plane, yielding a 2D-line. Also you may need to project the circle center in the circle plane, yielding a 2D-coordinate for the circle center, depending on where you center the coordinate system.
Then I simply solved this as a 2D line-circle intersection with the 2D center and radius.
This yields 0, 1 or 2 solutions which I projected back into 3D space (in the plane of the circle).
*Project = Take two arbitrary orthogonal axes in the plane (that are orthogonal to the normal). Each 2D dimension is the dot product of the corresponding axis and the 3D coordinate. To "reverse project", multiply each dimension with the corresponding axis and add them together. Also to move the reverse projected plane into the circle plane (restoring the information that was lost when 2D projecting), add the dot product of the circle center and the circle normal.