SuperLU is a general purpose library for the direct solution of large, sparse, nonsymmetric systems of linear equations on high performance machines. The library is written in C and is callable from either C or Fortran. The library routines will perform an LU decomposition with partial pivoting and triangular system solves through forward and back substitution. The LU factorization routines can handle non-square matrices but the triangular solves are performed only for square matrices. The matrix columns may be preordered (before factorization) either through library or user supplied routines. This preordering for sparsity is completely separate from the factorization. Working precision iterative refinement subroutines are provided for improved backward stability. Routines are also provided to equilibrate the system, estimate the condition number, calculate the relative backward error, and estimate error bounds for the refined solutions.
cordic methods describe essentially the same algorithm that with suitably chosen inputs can be used to calculate a whole range of scientific functions including sin, cos, tan, arctan, arcsin, arccos, sinh, cosh, tanh, arctanh, log, exp, square root and even multiply and divide.
the method dates back to volder [1959], and due to its versatility and compactness, it made possible the microcoding of the hp35 pocket scientific calculator in 1972.
here is some code to illustrate the techniques. ive split the methods into three parts linear, circular and hyperbolic. in the hp35 microcode these would be unified into one function (for space reasons). because the linear mode can perform multiply and divide, you only need add/subtract and shift to complete the implementation.
you can select in the code whether to do the multiples and divides also by cordic means. other multiplies and divides are all powers of 2 (these dont count). to eliminate these too, would involve ieee hackery.
The radially averaged power spectrum (RAPS) is the direction-independent mean spectrum, i.e. the average of all possible directional power spectra. This function computes and plots the RAPS of an input image. The function was written to produce the RAPS for square weather radar reflectivity images so you will have to relabel the plot axes and title to suit your needs. The radially averaged power spectrum provides a convenient means to view and compare information contained in 2D spectra in 1D
Main program running when workpiece is ready on deferent belt(deferent_ready=ture).
* Call Square_Wave subroutine to generate 0.5ms square wave on P1.2 to drive
* electromotor,then drive deferent belt step forward. When it steps to the measure
* zone, it stops to be measured. Then call A_D subroutine to transform analog
* signals to digital signals , after then call serial subroutine to transfer
* digital signals to PC. Call square wave subroutine to drive deferent belt step to
* original position waitting for defere ready flag to run the next circle.
