9 static char name[M][B];
15 scoremtx = NOTSPECIFIED;
18 PreRead( stdin, &njob, &nlenmax );
24 seq = AllocateCharMtx( njob, nlenmax+1 );
25 mtx = AllocateDoubleMtx( njob, njob );
28 FRead( stdin, name, nlen, seq );
30 readData( stdin, name, nlen, seq );
33 for( i=0; i<njob-1; i++ )
35 fprintf( stderr, "%4d/%4d\r", i+1, njob );
36 for( j=i+1; j<njob; j++ )
37 mtx[i][j] = (double)substitution_score( seq[i], seq[j] );
41 for( i=0; i<njob-1; i++ ) for( j=i+1; j<njob; j++ )
42 fprintf( stdout, "i=%d, j=%d, mtx[][] = %f\n", i, j, mtx[i][j] );
45 fp = fopen( "hat2", "w" );
46 WriteHat2( fp, njob, name, mtx );
50 res = system( ALNDIR "/spgsdl < hat2" );