BASIS_COMPARE
Compare two basis sets.

BASIS_COMPARE is an executable FORTRAN90 program, using double precision arithmetic, that can compare the amount of overlap between two basis sets.

In an M-dimensional space, a set of N orthonormal vectors U is provided, which span some subspace. We let Q be the M by N matrix whose columns are an orthogonal basis for the span of the vectors U(1:N). (We'll assume for now that the vectors U are linearly independent. Otherwise, Q could actually have fewer "interesting" columns than U.)

One or more M-vectors V are then considered. For each V, it is desired to know the proportion of the vector that lies within and outside the subspace spanned by U. V admits an orthogonal decomposition of the form:

        V = V1 + V2
      

• V1 is a vector in the subspace spanned by U;
• V2 is a vector orthogonal to the subspace spanned by U;
• sqrt ( ||V1||^2 + ||V2||^2 ) = ||V|| = Euclidean length of V.

The values

• ALPHA = ||V1|| / ||V||
• BETA = ||V2|| / ||V||

        TAN ( THETA ) = BETA / ALPHA.
      

To compute ALPHA, we realize that Q' * V produces the N by 1 vector of projection coefficients of V onto each U(I). Therefore, V1, the orthogonal projection of V into U is Q * ( Q' * V ):

        V1 =       Q * Q'   * V
      

        V2 = ( I - Q * Q' ) * V
      

Such sets are computed, for example, by CVT_BASIS_FLOW, POD_BASIS_FLOW, and SVD_BASIS, However, somewhat inconveniently, these programs are working with velocities, and store each vector in a pair of columns. For now, we have set up the program to expect the input data to have this format. Moreover, rather than comparing the two bases directly, we read in one basis to use as U, and then each vector of the other basis is checked against U. Also, U and V already consist of orthonormal columns, so the QR factorization of U is simply

        U = Q * R = U * I.
      

Related Data and Programs:

CVT_BASIS_FLOW is a FORTRAN90 program which extracts a basis from a set of velocity snapshots using CVT methods.

POD_BASIS_FLOW is a FORTRAN90 program which extracts a basis from a set of velocity snapshots using POD methods.

SVD_BASIS is a FORTRAN90 program which extracts a basis from any set of vectors, using POD methods.

Source Code:

• basis_compare.f90, the source code.
• basis_compare.csh, commands to compile and load the source code.

Examples and Tests:

• basis_compare_input.txt, sample input to the the program.
• basis_compare_output.txt, output from the program in response to the sample input.

List of Routines:

• MAIN is the main program for BASIS_COMPARE.
• CH_CAP capitalizes a single character.
• CH_EQI is a case insensitive comparison of two characters for equality.
• CH_IS_DIGIT returns .TRUE. if a character is a decimal digit.
• CH_TO_DIGIT returns the integer value of a base 10 digit.
• DATA_D2_READ reads a data set of pairs of real numbers stored in a file.
• DIGIT_INC increments a decimal digit.
• DIGIT_TO_CH returns the character representation of a decimal digit.
• FILE_NAME_INC generates the next filename in a series.
• FILE_ROW_COUNT counts the number of row records in a file.
• GET_UNIT returns a free FORTRAN unit number.
• S_INPUT prints a prompt string and reads a string from the user.
• S_TO_R8 reads an R8 from a string.
• TIMESTAMP prints the current YMDHMS date as a time stamp.
• TIMESTRING writes the current YMDHMS date into a string.

You can go up one level to the FORTRAN90 source codes.