コード例 #1
0
int correctDigit(const char *orifile, const char *recofile)
{
	FILE *fp1 = fopen(orifile, "r");
	FILE *fp2 = fopen(recofile, "r");

	int ret = 0;
	char str1[50], str2[50];
	while (fscanf(fp1, "%s", str1) == 1) {
		fscanf(fp2, "%s", str2);
		ret += correctChar(str1, str2);
	}

	fclose(fp1);
	fclose(fp2);
	return ret;
}
コード例 #2
0
bool RMatrix::operator () (const LongInt& i, RVector& w, const LongInt& j, const RVector& v)
 {
    bool value = true;

    integer incx = 1;
    integer m    = numberOfRows();
    integer n    = numberOfColumns();
    integer lda  = m;

    // dgemv (  correctChar(), &m, &n, (LongReal*) &(Complex::complexUnit), (LongReal*) (&_pelm[0]), 
    //         &lda, (LongReal*) &v(j), &incx, (LongReal*) &(Complex::complexUnit), 
    //         (LongReal*) &w(i), &incx
    //       );
    TensorCalculus::Blas<double>::gemv(*correctChar(), m, n, 1.0, &_pelm[0], 
            lda, &v(j), incx, 1.0, &w(i), incx);

    clearTranspose();
   return value;
 }
コード例 #3
0
RVector RMatrix :: operator * (const RVector& v)
 {
    integer incx = 1;
    integer m    = attr_numberOfRows;
    integer n    = attr_numberOfColumns;
    integer lda  = attr_numberOfRows;
    
	RVector result (attr_numberOfRows);

	// dgemv ( correctChar(), &m, &n, (LongReal*) &(Complex::complexUnit), (LongReal*) (&_pelm[0]), 
        //     &lda, (LongReal*) &v(0), &incx, (LongReal*) &(Complex::complexNull), 
        //     (LongReal*)&result(0), &incx
        //   );
	TensorCalculus::Blas<double>::gemv ( *correctChar(), m, n, 1.0, (&_pelm[0]), 
            lda, &v(0), incx, 0.0, &result(0), incx);

        clearTranspose();
   return result;
 }