1570 int projecteqstate(PetscReal wc[], const PetscReal ueq[], PetscReal lv[][3])  in projecteqstate()  argument
1576     for (j = 0; j < 3; ++j) wc[k] += lv[k][j] * ueq[j];  in projecteqstate()
1592 int eigenvectors(PetscReal rv[][3], PetscReal lv[][3], const PetscReal ueq[], PetscReal gamma)  in eigenvectors()  argument
1603   rho = ueq[0];  in eigenvectors()
1605   p0       = ueq[2];  in eigenvectors()
1629   PetscReal eps, ueq[3], rho0, twopi;  in initLinearWave()  local
1637   ueq[0] = rho0;  in initLinearWave()
1638   ueq[1] = u0;  in initLinearWave()
1639   ueq[2] = p0;  in initLinearWave()
1641   eigenvectors(rv, lv, ueq, gamma);  in initLinearWave()
1642   projecteqstate(wc, ueq, lv);  in initLinearWave()