// precomputed the values for some
// polynomial degree
double const *loqv_ptr=0;
- if (degree==2)
- {
- static const double loqv2[1*8]
- ={1/16., 1/16., 1/16., 1/16., 3/16., 3/16., 3/16., 3/16.};
- loqv_ptr=&loqv2[0];
- }
- else if (degree==3)
+ switch (degree)
{
- static const double loqv3[4*12]
- ={80/1053., 1/81., 11/1053., 1/81., 25/117., 44/351.,
- 7/117., 16/351., 7/117., 16/351., 25/117., 44/351.,
- 1/81., 80/1053., 1/81., 11/1053., 44/351., 25/117.,
- 25/117., 44/351., 16/351., 7/117., 7/117., 16/351.,
- 1/81., 11/1053., 1/81., 80/1053., 7/117., 16/351.,
- 16/351., 7/117., 25/117., 44/351., 44/351., 25/117.,
- 11/1053., 1/81., 80/1053., 1/81., 16/351., 7/117.,
- 44/351., 25/117., 44/351., 25/117., 16/351., 7/117.};
-
- loqv_ptr=&loqv3[0];
- }
+ case 2:
+ {
+ static const double loqv2[1*8]
+ ={1/16., 1/16., 1/16., 1/16., 3/16., 3/16., 3/16., 3/16.};
+ loqv_ptr=&loqv2[0];
+
+ break;
+ }
+ case 3:
+ {
+ static const double loqv3[4*12]
+ ={80/1053., 1/81., 11/1053., 1/81., 25/117., 44/351.,
+ 7/117., 16/351., 7/117., 16/351., 25/117., 44/351.,
+ 1/81., 80/1053., 1/81., 11/1053., 44/351., 25/117.,
+ 25/117., 44/351., 16/351., 7/117., 7/117., 16/351.,
+ 1/81., 11/1053., 1/81., 80/1053., 7/117., 16/351.,
+ 16/351., 7/117., 25/117., 44/351., 44/351., 25/117.,
+ 11/1053., 1/81., 80/1053., 1/81., 16/351., 7/117.,
+ 44/351., 25/117., 44/351., 25/117., 16/351., 7/117.};
+ loqv_ptr=&loqv3[0];
+
+ break;
+ }
+
+ // no other cases implemented,
+ // so simply fall through
+ }
+
if (loqv_ptr!=0)
{
// precomputed. copy values to