1. Use actual Solve Op to infer output dtype as CholSolve outputs a different dtype than basic Solve in Scipy==1.15

2. Tweaked test related to https://github.com/pymc-devs/pytensor/issues/1152

3. Tweak tolerage