From 7e9efcd2e4446f6c51f457640d1a34f52d4b1ab3 Mon Sep 17 00:00:00 2001 From: Wojtek Kosior Date: Wed, 24 Jul 2019 13:45:58 +0200 Subject: correct polynomial order for gauss --- src/quadratures.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) (limited to 'src/quadratures.f90') diff --git a/src/quadratures.f90 b/src/quadratures.f90 index 6fbb5cb..0fcaa2c 100644 --- a/src/quadratures.f90 +++ b/src/quadratures.f90 @@ -198,11 +198,11 @@ CONTAINS procedure(quadrature), pointer :: quad SELECT CASE (p) - CASE (1) + CASE (0) quad => gauss_n1 - CASE (2) + CASE (1) quad => gauss_n2 - CASE (3) + CASE (2) quad => gauss_n3 CASE default val = ieee_value(val, ieee_quiet_nan) -- cgit v1.2.3