Hi, I'm looking at using Tpetra for some problems where I'd like higher precision (e.g., quad). I'm currently using lapack's dgesv for a small (4x4 or 5x5) dense linear solve on the inside of my algorithm. Can anyone recommend a replacement for dgesv that supports templated scalar types? Thanks in advance ~ Pat