Previous Program Introduction Next Program

Program Formtovec

Usage

Formtovec 'file1' 'file2' [-m] [-d]
file1: matrix_TYP
file2: matrix_TYP or bravais_TYP.

Description

For each matrix A in file1 a vector V is calculated with the following property:

A = 1/V[NO+1] * (V[1] * F1 + V[2] * F2 +...+ V[NO] * FNO)

where Fi are the matrices in 'file2' if 'file2' is a matrix_TYP, otherwise are the matrices describing the form space of the bravais group in 'file2'.

Options

 -m: Use a modular (but exact) method to calculate
     the result. The result is calculated for a couple
     of big primes and fitted together with the chinese
     remainder theorem.
     This method is much faster and avoids trouble with
     overflow, BUT IS ABLE TO HANDLE THE CASE V[NO+1] == 1
     ONLY (AND WILL RUN INTO AN INFINITE LOOP OTHERWISE).
 -d: give the denominator in the NO+1-th colunm.

Remarks

CAUTION: if not used with the option -d, the denominator is not printed, so you will get a vector with only NO columns.

See also for Vectoform.

Examples

Previous Program Introduction Next Program