352   register double *wavep;  in freqwave()  local
358   wavep = wave;  in freqwave()
362       *wavep = C/(*freqp);  in freqwave()
370     wavep += swave;  in freqwave()
390   register const double *wavep;  in wavefreq()  local
396   wavep = wave;  in wavefreq()
400     if (*wavep != 0.0) {  in wavefreq()
401       *freqp = C/(*wavep);  in wavefreq()
408     wavep += swave;  in wavefreq()
550   register const double *wavep;  in waveawav()  local
556   wavep = wave;  in waveawav()
560     if (*wavep != 0.0) {  in waveawav()
563         s  = n/(*wavep);  in waveawav()
570       *awavp = (*wavep)/n;  in waveawav()
577     wavep += swave;  in waveawav()
600   register double *wavep;  in awavwave()  local
606   wavep = wave;  in awavwave()
615       *wavep = (*awavp)*n;  in awavwave()
623     wavep += swave;  in awavwave()
645   register const double *wavep;  in wavevelo()  local
650   wavep = wave;  in wavevelo()
654     s = *wavep * *wavep;  in wavevelo()
658     wavep += swave;  in wavevelo()
681   register double *wavep;  in velowave()  local
684   wavep = wave;  in velowave()
689       *wavep = restwav*sqrt((C + *velop)/s);  in velowave()
697     wavep += swave;  in velowave()
1037   register const double *wavep;  in wavevopt()  local
1045   wavep = wave;  in wavevopt()
1049     *voptp = r*(*wavep) - C;  in wavevopt()
1051     wavep += swave;  in wavevopt()
1073   register double *wavep;  in voptwave()  local
1078   wavep = wave;  in voptwave()
1081     *wavep = r*(C + *voptp);  in voptwave()
1084     wavep += swave;  in voptwave()
1106   register const double *wavep;  in wavezopt()  local
1114   wavep = wave;  in wavezopt()
1118     *zoptp = r*(*wavep) - 1.0;  in wavezopt()
1120     wavep += swave;  in wavezopt()
1141   register double *wavep;  in zoptwave()  local
1144   wavep = wave;  in zoptwave()
1147     *wavep = restwav*(1.0 + *zoptp);  in zoptwave()
1150     wavep += swave;  in zoptwave()