1 /* FCE Ultra - NES/Famicom Emulator
3 * Copyright notice for this file:
4 * Copyright (C) 2002 Ben Parnell
6 * This program is free software; you can redistribute it and/or modify
7 * it under the terms of the GNU General Public License as published by
8 * the Free Software Foundation; either version 2 of the License, or
9 * (at your option) any later version.
11 * This program is distributed in the hope that it will be useful,
12 * but WITHOUT ANY WARRANTY; without even the implied warranty of
13 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14 * GNU General Public License for more details.
16 * You should have received a copy of the GNU General Public License
17 * along with this program; if not, write to the Free Software
18 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
21 /********************************************************/
24 /******* Sound emulation code and waveform synthesis */
25 /******* routines. A few ideas were inspired */
26 /******* by code from Marat Fayzullin's EMUlib */
28 /********************************************************/
46 int32 WaveFinal[2048];
47 int16 WaveFinalMono[2048];
49 EXPSOUND GameExpSound={0,0,0};
58 static int32 count[5];
59 static int64 sqacc[2]={0,0};
71 static const uint8 Slengthtable[0x20]=
73 0x5,0x7f,0xA,0x1,0x14,0x2,0x28,0x3,0x50,0x4,0x1E,0x5,0x7,0x6,0x0E,0x7,
74 0x6,0x08,0xC,0x9,0x18,0xa,0x30,0xb,0x60,0xc,0x24,0xd,0x8,0xe,0x10,0xf
77 static uint32 lengthtable[0x20];
79 static const uint32 SNoiseFreqTable[0x10]=
81 2,4,8,0x10,0x20,0x30,0x40,0x50,0x65,0x7f,0xbe,0xfe,0x17d,0x1fc,0x3f9,0x7f2
83 static uint32 NoiseFreqTable[0x10];
87 static const uint8 NTSCPCMTable[0x10]=
89 0xd6,0xbe,0xaa,0xa0,0x8f,0x7f,0x71,0x6b,
90 0x5f,0x50,0x47,0x40,0x35,0x2a,0x24,0x1b
93 static const uint8 PALPCMTable[0x10]= // These values are just guessed.
95 0xc6,0xb0,0x9d,0x94,0x84,0x75,0x68,0x63,
96 0x58,0x4a,0x41,0x3b,0x31,0x27,0x21,0x19
102 // $4011 - Actual data outputted
103 // $4012 - Address register: $c000 + V*64
104 // $4013 - Size register: Size in bytes = (V+1)*64
107 static int64 PCMacc=0;
111 uint32 PCMAddressIndex=0;
112 int32 PCMSizeIndex=0;
116 static void Dummyfunc(void) {};
118 static void (*DoNoise)(void)=Dummyfunc;
119 static void (*DoTriangle)(void)=Dummyfunc;
120 static void (*DoPCM)(void)=Dummyfunc;
121 static void (*DoSQ1)(void)=Dummyfunc;
122 static void (*DoSQ2)(void)=Dummyfunc;
124 static void CalcDPCMIRQ(void)
131 freq=(PALPCMTable[PSG[0x10]&0xF]<<4);
133 freq=(NTSCPCMTable[PSG[0x10]&0xF]<<4);
135 cycles=(((PSG[0x13]<<4)+1));
137 honk=((PSG[0x13]<<4)+1)*freq;
140 //else honk/=(double)113.66666666;
142 //PCMIRQCount=honk*3; //180;
143 //if(PAL) PCMIRQCount*=.93;
147 static void PrepDPCM()
149 PCMAddressIndex=0x4000+(PSG[0x12]<<6);
150 PCMSizeIndex=(PSG[0x13]<<4)+1;
152 //PCMBuffer=ARead[0x8000+PCMAddressIndex](0x8000+PCMAddressIndex);
154 PCMfreq=PALPCMTable[PSG[0x10]&0xF];
156 PCMfreq=NTSCPCMTable[PSG[0x10]&0xF];
157 PCMacc=(int64)PCMfreq<<50;
160 uint8 sweepon[2]={0,0};
161 int32 curfreq[2]={0,0};
167 uint8 DecCountTo1[3];
175 /* Instantaneous? Maybe the new freq value is being calculated all of the time... */
176 static int FASTAPASS(2) CheckFreq(uint32 cf, uint8 sr)
188 static DECLFW(Write0x11)
194 static uint8 DutyCount[2]={0,0};
196 static DECLFW(Write_PSG)
198 //if((A>=0x4004 && A<=0x4007) || A==0x4015)
199 //printf("$%04x:$%02x, %d\n",A,V,SOUNDTS);
220 lengthcount[0]=lengthtable[(V>>3)&0x1f];
223 sweepon[0]=PSG[1]&0x80;
224 curfreq[0]=PSG[0x2]|((V&7)<<8);
226 DecCountTo1[0]=(PSG[0]&0xF)+1;
227 SweepCount[0]=((PSG[0x1]>>4)&7)+1;
229 sqacc[0]=((int64)curfreq[0]+1)<<50;
249 lengthcount[1]=lengthtable[(V>>3)&0x1f];
252 sweepon[1]=PSG[0x5]&0x80;
253 curfreq[1]=PSG[0x6]|((V&7)<<8);
255 DecCountTo1[1]=(PSG[0x4]&0xF)+1;
256 SweepCount[1]=((PSG[0x5]>>4)&7)+1;
258 sqacc[1]=((int64)curfreq[1]+1)<<50;
271 case 0xa:DoTriangle();
278 lengthcount[2]=lengthtable[(V>>3)&0x1f];
281 tricoop=PSG[0x8]&0x7f;
282 trimode=PSG[0x8]&0x80;
288 case 0xE:DoNoise();break;
294 lengthcount[3]=lengthtable[(V>>3)&0x1f];
297 DecCountTo1[2]=(PSG[0xC]&0xF)+1;
301 X6502_IRQEnd(FCEU_IQDPCM);
320 if(!(PSG[0x15]&0x10))
330 X6502_IRQEnd(FCEU_IQDPCM);
339 X6502_IRQEnd(FCEU_IQFCOUNT);
351 ret=(PSG[0x15]&(sqnon|0x10))|SIRQStat;
353 X6502_IRQEnd(/*FCEU_IQDPCM|*/FCEU_IQFCOUNT);
357 DECLFR(Read_PSGDummy)
361 ret=(PSG[0x15]&sqnon)|SIRQStat;
363 X6502_IRQEnd(/*FCEU_IQDPCM|*/FCEU_IQFCOUNT);
367 static void FASTAPASS(1) FrameSoundStuff(int V)
377 case 1: /* Envelope decay, linear counter, length counter, freq sweep */
378 if(PSG[0x15]&4 && sqnon&4)
384 if(lengthcount[2]<=0)
394 if(PSG[0x15]&(P+1) && sqnon&(P+1))
396 if(!(PSG[P<<2]&0x20))
401 if(lengthcount[P]<=0)
408 /* Frequency Sweep Code Here */
410 /* xxxx = hz. 120/(x+1)*/
415 if(SweepCount[P]>0) SweepCount[P]--;
418 SweepCount[P]=((PSG[(P<<2)+0x1]>>4)&7)+1; //+1;
420 if(PSG[(P<<2)+0x1]&0x8)
422 mod-=(P^1)+((curfreq[P])>>(PSG[(P<<2)+0x1]&7));
424 if(curfreq[P] && (PSG[(P<<2)+0x1]&7)/* && sweepon[P]&0x80*/)
431 mod=curfreq[P]>>(PSG[(P<<2)+0x1]&7);
432 if((mod+curfreq[P])&0x800)
439 if(curfreq[P] && (PSG[(P<<2)+0x1]&7)/* && sweepon[P]&0x80*/)
450 if(PSG[0x15]&0x8 && sqnon&8)
457 if(lengthcount[3]<=0)
465 case 0: /* Envelope decay + linear counter */
471 if(tricoop==1) DoTriangle();
478 if(DecCountTo1[P]>0) DecCountTo1[P]--;
479 if(DecCountTo1[P]<=0)
481 DecCountTo1[P]=(PSG[P<<2]&0xF)+1;
482 if(decvolume[P] || PSG[P<<2]&0x20)
485 /* Step from 0 to full volume seems to take twice as long
486 as the other steps. I don't know if this is the correct
487 way to double its length, though(or if it even matters).
489 if((PSG[P<<2]&0x20) && (decvolume[P]==0))
494 if(!(PSG[P<<2]&0x10))
495 realvolume[P]=decvolume[P];
498 if(DecCountTo1[2]>0) DecCountTo1[2]--;
499 if(DecCountTo1[2]<=0)
501 DecCountTo1[2]=(PSG[0xC]&0xF)+1;
502 if(decvolume[2] || PSG[0xC]&0x20)
505 /* Step from 0 to full volume seems to take twice as long
506 as the other steps. I don't know if this is the correct
507 way to double its length, though(or if it even matters).
509 if((PSG[0xC]&0x20) && (decvolume[2]==0))
515 realvolume[2]=decvolume[2];
522 void FrameSoundUpdate(void)
524 // Linear counter: Bit 0-6 of $4008
525 // Length counter: Bit 4-7 of $4003, $4007, $400b, $400f
531 if(!(PSG[0x17]&0xC0))
534 X6502_IRQBegin(FCEU_IQFCOUNT);
537 //if(SIRQStat&0x40) X6502_IRQBegin(FCEU_IQFCOUNT);
538 FrameSoundStuff(fcnt);
542 static uint32 ChannelBC[5];
544 static uint32 RectAmp[2][8];
546 static void FASTAPASS(1) CalcRectAmp(int P)
548 static int tal[4]={1,2,4,6};
551 uint32 *b=RectAmp[P];
556 //V=(PSG[P<<2]&15)<<4;
558 // V=decvolume[P]<<4;
559 m=tal[(PSG[P<<2]&0xC0)>>6];
566 static void RDoPCM(void)
571 uint32 out=PSG[0x11]<<3;
574 end=(SOUNDTS<<16)/soundtsinc;
575 if(end<=start) return;
583 for(V=start;V<end;V++)
585 PCMacc-=nesincsizeLL;
599 Wave[V>>4]+=PSG[0x11]<<3;
605 PCMBuffer=ARead[0x8000+PCMAddressIndex](0x8000+PCMAddressIndex);
606 PCMAddressIndex=(PCMAddressIndex+1)&0x7fff;
611 int t=(((PCMBuffer>>PCMBitIndex)&1)<<2)-2;
621 PCMBitIndex=(PCMBitIndex+1)&7;
623 Wave[V>>4]+=out; //(PSG[0x11]-64)<<3;
630 for(V=start;V<=(start|15);V++)
633 for(V=(start>>4)+1;V<(end>>4);V++)
636 for(V=end&(~15);V<end;V++)
640 for(V=start;V<end;V++)
646 static void RDoSQ1(void)
654 end=(SOUNDTS<<16)/soundtsinc;
655 if(end<=start) return;
658 if(curfreq[0]<8 || curfreq[0]>0x7ff)
660 if(!CheckFreq(curfreq[0],PSG[0x1]))
663 if(PSG[0x15]&1 && sqnon&1)
665 uint32 out=RectAmp[0][DutyCount[0]];
669 for(V=start;V<end;V++)
672 sqacc[0]-=nesincsizeLL;
678 if(sqacc[0]<=0) goto rea;
681 out=RectAmp[0][DutyCount[0]];
689 static void RDoSQ2(void)
697 end=(SOUNDTS<<16)/soundtsinc;
698 if(end<=start) return;
701 if(curfreq[1]<8 || curfreq[1]>0x7ff)
703 if(!CheckFreq(curfreq[1],PSG[0x5]))
706 if(PSG[0x15]&2 && sqnon&2)
708 uint32 out=RectAmp[1][DutyCount[1]];
713 for(V=start;V<end;V++)
716 sqacc[1]-=nesincsizeLL;
722 if(sqacc[1]<=0) goto rea;
725 out=RectAmp[1][DutyCount[1]];
734 static void RDoTriangle(void)
736 static uint32 tcout=0;
738 int32 start,end; //,freq;
739 int64 freq=(((PSG[0xa]|((PSG[0xb]&7)<<8))+1));
742 end=(SOUNDTS<<16)/soundtsinc;
743 if(end<=start) return;
746 if(! (PSG[0x15]&0x4 && sqnon&4 && tricoop) )
747 { // Counter is halted, but we still need to output.
748 for(V=start;V<end;V++)
751 else if(freq<=4) // 55.9Khz - Might be barely audible on a real NES, but
752 // it's too costly to generate audio at this high of a frequency
753 // (55.9Khz * 32 for the stepping).
754 // The same could probably be said for ~27.8Khz, so we'll
755 // take care of that too. We'll just output the average
756 // value(15/2 - scaled properly for our output format, of course).
757 // We'll also take care of ~18Khz and ~14Khz too, since they should be barely audible.
758 // (Some proof or anything to confirm/disprove this would be nice.).
760 for(V=start;V<end;V++)
761 Wave[V>>4]+=((0xF<<4)+(0xF<<2))>>1;
765 static int64 triacc=0;
769 for(V=start;V<end;V++)
771 triacc-=nesincsizeLL;
777 if(triacc<=0) goto rea;
780 if(tc&0x10) tcout^=0xF;
781 tcout=(tcout<<4)+(tcout<<2);
788 static void RDoNoise(void)
794 end=(SOUNDTS<<16)/soundtsinc;
795 if(end<=start) return;
798 if(PSG[0x15]&0x8 && sqnon&8)
804 amplitude=realvolume[2];
806 // amplitude=(PSG[0xC]&0xF);
808 // amplitude=decvolume[2]&0xF;
810 inc=NoiseFreqTable[PSG[0xE]&0xF];
811 amptab[0]=((amplitude<<2)+amplitude+amplitude)<<1;
817 if(PSG[0xE]&0x80) // "short" noise
818 for(V=start;V<end;V++)
825 feedback=((nreg>>8)&1)^((nreg>>14)&1);
826 nreg=(nreg<<1)+feedback;
834 for(V=start;V<end;V++)
841 feedback=((nreg>>13)&1)^((nreg>>14)&1);
842 nreg=(nreg<<1)+feedback;
854 void SetNESSoundMap(void)
856 SetWriteHandler(0x4000,0x4013,Write_PSG);
857 SetWriteHandler(0x4011,0x4011,Write0x11);
858 SetWriteHandler(0x4015,0x4015,Write_PSG);
859 SetWriteHandler(0x4017,0x4017,Write_PSG);
860 SetReadHandler(0x4015,0x4015,Read_PSG);
863 static int32 WaveNSF[256];
865 int32 highp; // 0 through 65536, 0 = no high pass, 65536 = max high pass
867 int32 lowp; // 0 through 65536, 65536 = max low pass(total attenuation)
868 // 65536 = no low pass
869 static void FilterSound(uint32 *in, int32 *out, int16 *outMono, int count)
871 static int64 acc=0, acc2=0;
874 //int16* outorig=out;
876 for(;count;count--,in++)//,out++)//,index++)
880 diff=((int64)*in<<24)-acc;
882 acc+=(diff*highp)>>16;
883 acc2+=((diff-acc2)*lowp)>>16;
886 // don't change the sound here
887 // *out=(acc2*(int64)FSettings.SoundVolume)>>(24+16);
888 // volume, 4 times louder by default??
889 // *out = acc2 >> 24;
890 // just a bit louder. Hope it's okay
893 if(*out<-32767) *out=-32767;
894 if(*out>32767) *out=32767;
898 tmp=(int16 *)(out-1);
899 // don't do this the first time
900 if (prev == -99999) continue;
901 // the middle one should be interpolated
902 tmp[1]=(int16)((*out + prev) >> 1);
905 //outMono[index] = (int16)*out;
906 *outMono = (int16)(acc2 >> 24);
907 //if(*outMono<-16384) *outMono=-16384;
908 //if(*outMono>16384) *outMono=16384;
911 // out=((int64)(acc2>>24)*(int64)FSettings.SoundVolume)>>16; //acc2>>24;
920 int FlushEmulateSound(void)
925 if(!timestamp) return(0);
927 if(!FSettings.SndRate || (soundvol == 0))
933 end=(SOUNDTS<<16)/soundtsinc;
940 if(GameExpSound.Fill)
941 GameExpSound.Fill(end&0xF);
943 // FilterSound(Wave,WaveFinal,end>>4);
944 FilterSound(Wave,WaveFinal,WaveFinalMono,end>>4);
945 // printf("count %d, num ints %d\n", end, (end >> 4));
946 if(FCEUGameInfo.type==GIT_NSF)
949 int x,s=0,si=end/1024; // Only want 1/4 of the output buffer to be displayed
952 WaveNSF[x]=WaveFinal[s>>4];
958 Wave[0]=Wave[(end>>4)];
963 ChannelBC[x]=end&0xF;
964 soundtsoffs=(soundtsinc*(end&0xF))>>16;
968 void GetSoundBuffer(int32 **W)
973 void PowerSound(void)
981 BWrite[0x4000+x](0x4000+x,0);
989 void ResetSound(void)
993 if(x!=1 && x!=5 && x!=0x14) BWrite[0x4000+x](0x4000+x,0);
1000 void SetSoundVariables(void)
1004 fhinc=PAL?16626:14915; // *2 CPU clock rate
1007 lengthtable[x]=Slengthtable[x]<<1;
1009 if(FSettings.SndRate)
1012 DoTriangle=RDoTriangle;
1019 DoNoise=DoTriangle=DoPCM=DoSQ1=DoSQ2=Dummyfunc;
1022 if(!FSettings.SndRate) return;
1023 if(GameExpSound.RChange)
1024 GameExpSound.RChange();
1026 nesincsizeLL=(int64)((int64)562949953421312*(long double)(PAL?PAL_CPU:NTSC_CPU)/(FSettings.SndRate OVERSAMPLE));
1027 PSG_base=(uint32)(PAL?(long double)PAL_CPU/16:(long double)NTSC_CPU/16);
1032 z=SNoiseFreqTable[x]<<1;
1033 z=(PAL?PAL_CPU:NTSC_CPU)/z;
1034 z=(long double)((uint32)((FSettings.SndRate OVERSAMPLE)<<12))/z;
1035 NoiseFreqTable[x]=z;
1037 soundtsinc=(uint32)((uint64)(PAL?(long double)PAL_CPU*65536:(long double)NTSC_CPU*65536)/(FSettings.SndRate OVERSAMPLE));
1038 memset(Wave,0,2048*4);
1041 highp=(250<<16)/FSettings.SndRate; // Arbitrary
1042 lowp=((int64)25000<<16)/FSettings.SndRate; // Arbitrary
1044 if(highp>(1<<16)) highp=1<<16;
1045 if(lowp>(1<<16)) lowp=1<<16;
1048 void FixOldSaveStateSFreq(void)
1053 curfreq[x]=PSG[0x2+(x<<2)]|((PSG[0x3+(x<<2)]&7)<<8);
1057 void FCEUI_Sound(int Rate)
1059 FSettings.SndRate=Rate;
1060 SetSoundVariables();
1063 void FCEUI_SetSoundVolume(uint32 volume)
1065 FSettings.SoundVolume=(volume<<16)/100;