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};
69 static const uint8 Slengthtable[0x20]=
71 0x5,0x7f,0xA,0x1,0x14,0x2,0x28,0x3,0x50,0x4,0x1E,0x5,0x7,0x6,0x0E,0x7,
72 0x6,0x08,0xC,0x9,0x18,0xa,0x30,0xb,0x60,0xc,0x24,0xd,0x8,0xe,0x10,0xf
75 static uint32 lengthtable[0x20];
77 static const uint32 SNoiseFreqTable[0x10]=
79 2,4,8,0x10,0x20,0x30,0x40,0x50,0x65,0x7f,0xbe,0xfe,0x17d,0x1fc,0x3f9,0x7f2
81 static uint32 NoiseFreqTable[0x10];
85 static const uint8 NTSCPCMTable[0x10]=
87 0xd6,0xbe,0xaa,0xa0,0x8f,0x7f,0x71,0x6b,
88 0x5f,0x50,0x47,0x40,0x35,0x2a,0x24,0x1b
91 static const uint8 PALPCMTable[0x10]= // These values are just guessed.
93 0xc6,0xb0,0x9d,0x94,0x84,0x75,0x68,0x63,
94 0x58,0x4a,0x41,0x3b,0x31,0x27,0x21,0x19
100 // $4011 - Actual data outputted
101 // $4012 - Address register: $c000 + V*64
102 // $4013 - Size register: Size in bytes = (V+1)*64
105 static int64 PCMacc=0;
109 uint32 PCMAddressIndex=0;
110 int32 PCMSizeIndex=0;
114 static void Dummyfunc(void) {};
116 static void (*DoNoise)(void)=Dummyfunc;
117 static void (*DoTriangle)(void)=Dummyfunc;
118 static void (*DoPCM)(void)=Dummyfunc;
119 static void (*DoSQ1)(void)=Dummyfunc;
120 static void (*DoSQ2)(void)=Dummyfunc;
122 static void CalcDPCMIRQ(void)
129 freq=(PALPCMTable[PSG[0x10]&0xF]<<4);
131 freq=(NTSCPCMTable[PSG[0x10]&0xF]<<4);
133 cycles=(((PSG[0x13]<<4)+1));
135 honk=((PSG[0x13]<<4)+1)*freq;
138 //else honk/=(double)113.66666666;
140 //PCMIRQCount=honk*3; //180;
141 //if(PAL) PCMIRQCount*=.93;
145 static void PrepDPCM()
147 PCMAddressIndex=0x4000+(PSG[0x12]<<6);
148 PCMSizeIndex=(PSG[0x13]<<4)+1;
150 //PCMBuffer=ARead[0x8000+PCMAddressIndex](0x8000+PCMAddressIndex);
152 PCMfreq=PALPCMTable[PSG[0x10]&0xF];
154 PCMfreq=NTSCPCMTable[PSG[0x10]&0xF];
155 PCMacc=(int64)PCMfreq<<50;
158 uint8 sweepon[2]={0,0};
159 int32 curfreq[2]={0,0};
165 uint8 DecCountTo1[3];
173 /* Instantaneous? Maybe the new freq value is being calculated all of the time... */
174 static int FASTAPASS(2) CheckFreq(uint32 cf, uint8 sr)
186 static DECLFW(Write0x11)
192 static uint8 DutyCount[2]={0,0};
194 static DECLFW(Write_PSG)
196 //if((A>=0x4004 && A<=0x4007) || A==0x4015)
197 //printf("$%04x:$%02x, %d\n",A,V,timestamp);
218 lengthcount[0]=lengthtable[(V>>3)&0x1f];
221 sweepon[0]=PSG[1]&0x80;
222 curfreq[0]=PSG[0x2]|((V&7)<<8);
224 DecCountTo1[0]=(PSG[0]&0xF)+1;
225 SweepCount[0]=((PSG[0x1]>>4)&7)+1;
227 sqacc[0]=((int64)curfreq[0]+1)<<50;
247 lengthcount[1]=lengthtable[(V>>3)&0x1f];
250 sweepon[1]=PSG[0x5]&0x80;
251 curfreq[1]=PSG[0x6]|((V&7)<<8);
253 DecCountTo1[1]=(PSG[0x4]&0xF)+1;
254 SweepCount[1]=((PSG[0x5]>>4)&7)+1;
256 sqacc[1]=((int64)curfreq[1]+1)<<50;
269 case 0xa:DoTriangle();
276 lengthcount[2]=lengthtable[(V>>3)&0x1f];
279 tricoop=PSG[0x8]&0x7f;
280 trimode=PSG[0x8]&0x80;
286 case 0xE:DoNoise();break;
292 lengthcount[3]=lengthtable[(V>>3)&0x1f];
295 DecCountTo1[2]=(PSG[0xC]&0xF)+1;
299 X6502_IRQEnd(FCEU_IQDPCM);
318 if(!(PSG[0x15]&0x10))
328 X6502_IRQEnd(FCEU_IQDPCM);
337 X6502_IRQEnd(FCEU_IQFCOUNT);
349 ret=(PSG[0x15]&(sqnon|0x10))|SIRQStat;
351 X6502_IRQEnd(/*FCEU_IQDPCM|*/FCEU_IQFCOUNT);
355 DECLFR(Read_PSGDummy)
359 ret=(PSG[0x15]&sqnon)|SIRQStat;
361 X6502_IRQEnd(/*FCEU_IQDPCM|*/FCEU_IQFCOUNT);
365 static void FASTAPASS(1) FrameSoundStuff(int V)
375 case 1: /* Envelope decay, linear counter, length counter, freq sweep */
376 if(PSG[0x15]&4 && sqnon&4)
382 if(lengthcount[2]<=0)
392 if(PSG[0x15]&(P+1) && sqnon&(P+1))
394 if(!(PSG[P<<2]&0x20))
399 if(lengthcount[P]<=0)
406 /* Frequency Sweep Code Here */
408 /* xxxx = hz. 120/(x+1)*/
413 if(SweepCount[P]>0) SweepCount[P]--;
416 SweepCount[P]=((PSG[(P<<2)+0x1]>>4)&7)+1; //+1;
418 if(PSG[(P<<2)+0x1]&0x8)
420 mod-=(P^1)+((curfreq[P])>>(PSG[(P<<2)+0x1]&7));
422 if(curfreq[P] && (PSG[(P<<2)+0x1]&7)/* && sweepon[P]&0x80*/)
429 mod=curfreq[P]>>(PSG[(P<<2)+0x1]&7);
430 if((mod+curfreq[P])&0x800)
437 if(curfreq[P] && (PSG[(P<<2)+0x1]&7)/* && sweepon[P]&0x80*/)
448 if(PSG[0x15]&0x8 && sqnon&8)
455 if(lengthcount[3]<=0)
463 case 0: /* Envelope decay + linear counter */
469 if(tricoop==1) DoTriangle();
476 if(DecCountTo1[P]>0) DecCountTo1[P]--;
477 if(DecCountTo1[P]<=0)
479 DecCountTo1[P]=(PSG[P<<2]&0xF)+1;
480 if(decvolume[P] || PSG[P<<2]&0x20)
483 /* Step from 0 to full volume seems to take twice as long
484 as the other steps. I don't know if this is the correct
485 way to double its length, though(or if it even matters).
487 if((PSG[P<<2]&0x20) && (decvolume[P]==0))
492 if(!(PSG[P<<2]&0x10))
493 realvolume[P]=decvolume[P];
496 if(DecCountTo1[2]>0) DecCountTo1[2]--;
497 if(DecCountTo1[2]<=0)
499 DecCountTo1[2]=(PSG[0xC]&0xF)+1;
500 if(decvolume[2] || PSG[0xC]&0x20)
503 /* Step from 0 to full volume seems to take twice as long
504 as the other steps. I don't know if this is the correct
505 way to double its length, though(or if it even matters).
507 if((PSG[0xC]&0x20) && (decvolume[2]==0))
513 realvolume[2]=decvolume[2];
520 void FrameSoundUpdate(void)
522 // Linear counter: Bit 0-6 of $4008
523 // Length counter: Bit 4-7 of $4003, $4007, $400b, $400f
529 if(!(PSG[0x17]&0xC0))
532 X6502_IRQBegin(FCEU_IQFCOUNT);
535 //if(SIRQStat&0x40) X6502_IRQBegin(FCEU_IQFCOUNT);
536 FrameSoundStuff(fcnt);
540 static uint32 ChannelBC[5];
542 static uint32 RectAmp[2][8];
544 static void FASTAPASS(1) CalcRectAmp(int P)
546 static int tal[4]={1,2,4,6};
549 uint32 *b=RectAmp[P];
554 //V=(PSG[P<<2]&15)<<4;
556 // V=decvolume[P]<<4;
557 m=tal[(PSG[P<<2]&0xC0)>>6];
564 static void RDoPCM(void)
569 uint32 out=PSG[0x11]<<3;
572 end=(timestamp<<16)/soundtsinc;
573 if(end<=start) return;
581 for(V=start;V<end;V++)
583 PCMacc-=nesincsizeLL;
597 Wave[V>>4]+=PSG[0x11]<<3;
603 PCMBuffer=ARead[0x8000+PCMAddressIndex](0x8000+PCMAddressIndex);
604 PCMAddressIndex=(PCMAddressIndex+1)&0x7fff;
609 int t=(((PCMBuffer>>PCMBitIndex)&1)<<2)-2;
619 PCMBitIndex=(PCMBitIndex+1)&7;
621 Wave[V>>4]+=out; //(PSG[0x11]-64)<<3;
628 for(V=start;V<=(start|15);V++)
631 for(V=(start>>4)+1;V<(end>>4);V++)
634 for(V=end&(~15);V<end;V++)
638 for(V=start;V<end;V++)
644 static void RDoSQ1(void)
652 end=(timestamp<<16)/soundtsinc;
653 if(end<=start) return;
656 if(curfreq[0]<8 || curfreq[0]>0x7ff)
658 if(!CheckFreq(curfreq[0],PSG[0x1]))
661 if(PSG[0x15]&1 && sqnon&1)
663 uint32 out=RectAmp[0][DutyCount[0]];
667 for(V=start;V<end;V++)
670 sqacc[0]-=nesincsizeLL;
676 if(sqacc[0]<=0) goto rea;
679 out=RectAmp[0][DutyCount[0]];
687 static void RDoSQ2(void)
695 end=(timestamp<<16)/soundtsinc;
696 if(end<=start) return;
699 if(curfreq[1]<8 || curfreq[1]>0x7ff)
701 if(!CheckFreq(curfreq[1],PSG[0x5]))
704 if(PSG[0x15]&2 && sqnon&2)
706 uint32 out=RectAmp[1][DutyCount[1]];
711 for(V=start;V<end;V++)
714 sqacc[1]-=nesincsizeLL;
720 if(sqacc[1]<=0) goto rea;
723 out=RectAmp[1][DutyCount[1]];
732 static void RDoTriangle(void)
734 static uint32 tcout=0;
736 int32 start,end; //,freq;
737 int64 freq=(((PSG[0xa]|((PSG[0xb]&7)<<8))+1));
740 end=(timestamp<<16)/soundtsinc;
741 if(end<=start) return;
744 if(! (PSG[0x15]&0x4 && sqnon&4 && tricoop) )
745 { // Counter is halted, but we still need to output.
746 for(V=start;V<end;V++)
749 else if(freq<=4) // 55.9Khz - Might be barely audible on a real NES, but
750 // it's too costly to generate audio at this high of a frequency
751 // (55.9Khz * 32 for the stepping).
752 // The same could probably be said for ~27.8Khz, so we'll
753 // take care of that too. We'll just output the average
754 // value(15/2 - scaled properly for our output format, of course).
755 // We'll also take care of ~18Khz and ~14Khz too, since they should be barely audible.
756 // (Some proof or anything to confirm/disprove this would be nice.).
758 for(V=start;V<end;V++)
759 Wave[V>>4]+=((0xF<<4)+(0xF<<2))>>1;
763 static int64 triacc=0;
767 for(V=start;V<end;V++)
769 triacc-=nesincsizeLL;
775 if(triacc<=0) goto rea;
778 if(tc&0x10) tcout^=0xF;
779 tcout=(tcout<<4)+(tcout<<2);
786 static void RDoNoise(void)
792 end=(timestamp<<16)/soundtsinc;
793 if(end<=start) return;
796 if(PSG[0x15]&0x8 && sqnon&8)
802 amplitude=realvolume[2];
804 // amplitude=(PSG[0xC]&0xF);
806 // amplitude=decvolume[2]&0xF;
808 inc=NoiseFreqTable[PSG[0xE]&0xF];
809 amptab[0]=((amplitude<<2)+amplitude+amplitude)<<1;
815 if(PSG[0xE]&0x80) // "short" noise
816 for(V=start;V<end;V++)
823 feedback=((nreg>>8)&1)^((nreg>>14)&1);
824 nreg=(nreg<<1)+feedback;
832 for(V=start;V<end;V++)
839 feedback=((nreg>>13)&1)^((nreg>>14)&1);
840 nreg=(nreg<<1)+feedback;
852 void SetNESSoundMap(void)
854 SetWriteHandler(0x4000,0x4013,Write_PSG);
855 SetWriteHandler(0x4011,0x4011,Write0x11);
856 SetWriteHandler(0x4015,0x4015,Write_PSG);
857 SetWriteHandler(0x4017,0x4017,Write_PSG);
858 SetReadHandler(0x4015,0x4015,Read_PSG);
861 static int32 WaveNSF[256];
863 int32 highp; // 0 through 65536, 0 = no high pass, 65536 = max high pass
865 int32 lowp; // 0 through 65536, 65536 = max low pass(total attenuation)
866 // 65536 = no low pass
867 static void FilterSound(uint32 *in, int32 *out, int16 *outMono, int count)
869 static int64 acc=0, acc2=0;
872 //int16* outorig=out;
874 for(;count;count--,in++)//,out++)//,index++)
878 diff=((int64)*in<<24)-acc;
880 acc+=(diff*highp)>>16;
881 acc2+=((diff-acc2)*lowp)>>16;
884 // don't change the sound here
885 // *out=(acc2*(int64)FSettings.SoundVolume)>>(24+16);
886 // volume, 4 times louder by default??
887 // *out = acc2 >> 24;
888 // just a bit louder. Hope it's okay
891 if(*out<-32767) *out=-32767;
892 if(*out>32767) *out=32767;
896 tmp=(int16 *)(out-1);
897 // don't do this the first time
898 if (prev == -99999) continue;
899 // the middle one should be interpolated
900 tmp[1]=(int16)((*out + prev) >> 1);
903 //outMono[index] = (int16)*out;
904 *outMono = (int16)(acc2 >> 24);
905 //if(*outMono<-16384) *outMono=-16384;
906 //if(*outMono>16384) *outMono=16384;
909 // out=((int64)(acc2>>24)*(int64)FSettings.SoundVolume)>>16; //acc2>>24;
918 int FlushEmulateSound(void)
924 if(!timestamp) return(0);
926 if(!FSettings.SndRate || (soundvol == 0))
932 end=(timestamp<<16)/soundtsinc;
939 if(GameExpSound.Fill)
940 GameExpSound.Fill(end&0xF);
942 // FilterSound(Wave,WaveFinal,end>>4);
943 FilterSound(Wave,WaveFinal,WaveFinalMono,end>>4);
944 // printf("count %d, num ints %d\n", end, (end >> 4));
945 if(FCEUGameInfo.type==GIT_NSF)
948 int x,s=0,si=end/1024; // Only want 1/4 of the output buffer to be displayed
951 WaveNSF[x]=WaveFinal[s>>4];
957 Wave[0]=Wave[(end>>4)];
962 ChannelBC[x]=end&0xF;
963 timestampbase+=timestamp;
964 timestamp=(soundtsinc*(end&0xF))>>16;
965 timestampbase-=timestamp;
969 void GetSoundBuffer(int32 **W)
974 void PowerSound(void)
982 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;