test.c

上传个模糊聚类代码

//parse DNA DATAS
#include <stdio.h>
#define SIZE 40

char* inputFileName="inputData.txt";
char* outputFileName="outputData.txt";

int main()
{	
	//DNA SEQUENCE:A,T,G,C
	int dataArr[SIZE][4];
	int count;
	char myChar;
	FILE* inputFile,* outputFile;
	int i,j;
	
	if((inputFile=fopen(inputFileName,"rb"))==NULL)
	{
		printf("open input file failure,exit!\n");
		getchar();
		return 1;
	}
	
	if((outputFile=fopen(outputFileName,"wb"))==NULL)
	{
		printf("open out file failure,exit!\n");
		getchar();
        return 1;
	}
	
	//data initialize
	for (i=0 ;i<SIZE;i++)
		for (j=0;j<4;j++)
			dataArr[i][j]=0;
	
	count=0;
	
	do
	{  	
        do
		{	
 	          do
 	          {
		   	          fscanf(inputFile,"%c",&myChar);
   	                  if(feof(inputFile))//should add this,or you will get no end of parse
  		                    break;
               }while(myChar=='\n'||myChar=='\t'&&myChar==' '||myChar==13||myChar==10);
			
   		    if(myChar=='a'||myChar=='A')
			{
				dataArr[count][0]++;		
			}
			else if(myChar=='t'||myChar=='T')
			{
				dataArr[count][1]++;
			}
			else if(myChar=='g'||myChar=='G')
			{
				dataArr[count][2]++;
			}
			else if(myChar=='c'||myChar=='C')
			{
				dataArr[count][3]++;
			}
			else if(myChar!=-1&&myChar!='\0'&&myChar!=48&&!feof(inputFile))
			{
				printf("\nabnormal char scaned:");
				printf("%d\n",myChar);
                return -1;
			}
			printf("%c",myChar);
		}
		while(myChar!=48);//end of a dna list line	
		
		       do
 	          {
		   	          fscanf(inputFile,"%c",&myChar);
   	                  if(feof(inputFile))//should add this,or you will get no end of parse
  		                    break;
               }while(myChar=='\n'||myChar=='\t'&&myChar==' '||myChar==13||myChar==10);
        
        count++;
	}while(!feof(inputFile));
		
	//output
	for(i=0;i<SIZE;i++)
	{
		for(j=0;j<3;j++)
			fprintf(outputFile,"%d ",dataArr[i][j]);
		fprintf(outputFile,"%d \r\n",dataArr[i][j]);
	}
	printf("end of test program\n");
	getchar();
	return 0;
}