| 461 | #endif |
| 462 | |
| 463 | void wNorm(int *degw, int *lpol, int npol, double *rel) |
| 464 | { |
| 465 | int i, j, ecu, ec; |
| 466 | int *ex; |
| 467 | double *r; |
| 468 | |
| 469 | ex = degw; |
| 470 | r = rel; |
| 471 | for (i = 0; i < npol; i++) |
| 472 | { |
| 473 | ecu = *ex++; |
| 474 | for (j = lpol[i] - 1; j!=0 ; j--) |
| 475 | { |
| 476 | ec = *ex++; |
| 477 | if (ec > ecu) |
| 478 | ecu = ec; |
| 479 | } |
| 480 | *r = (double)1.0 / (double)(ecu * ecu); |
| 481 | r++; |
| 482 | } |
| 483 | } |