Linear system solver in D?

[BCS]

I am going to have a system of equations like this

a_11*x_1 + a_12*x2 + ... a_1n*x_n = y_1
a_21*x_1 + a_22*x2 + ... a_2n*x_n = y_2
.
.
.
a_m1*x_1 + a_m2*x2 + ... a_mn*x_n = y_m

y_* and a_* known, I need to find x_*

What is the best available solver for such a system that works under D?

C bindings would work, D code would be better and I'd rather stay away from 
non portable (uses __ asm and has no port to other system).

If no one knows of a good lib that is ready to use, what is a good C lib 
to do bindings for?


p.s. I'm going to be putting this in a non-linear root finder, has someone 
already written on of those for D.