Commit 0ff70335 authored by Anne Reinarz's avatar Anne Reinarz

fix EV for alpha=0

parent 8026c92c
......@@ -204,7 +204,7 @@ void PDEncp(const double* const Q,const double* const gradQ,double* const BgradQ
}
void PDEEigenvalues(const double* const Q,const int direction,double* const L){
if(Q[5]<1.e-10){
if(Q[5]<1.e-4){
L[0] = 1.0;
return;
}
......
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment