1570 int projecteqstate(PetscReal wc[], const PetscReal ueq[], PetscReal lv[][3])  in projecteqstate()
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()
1600       lv[k][j] = 0.;  in eigenvectors()
1607   lv[0][1] = rho * .5;  in eigenvectors()
1608   lv[0][2] = -.5 / csnd;  in eigenvectors()
1609   lv[1][0] = csnd;  in eigenvectors()
1610   lv[1][2] = -1. / csnd;  in eigenvectors()
1611   lv[2][1] = rho * .5;  in eigenvectors()
1612   lv[2][2] = .5 / csnd;  in eigenvectors()
1626   PetscReal lv[3][3];  in initLinearWave()  local
1641   eigenvectors(rv, lv, ueq, gamma);  in initLinearWave()
1642   projecteqstate(wc, ueq, lv);  in initLinearWave()