problem in DSS_FACTOR_complex

haijunwu — Mon, 10 May 2010 14:35:22 GMT

I just modified the test program dss_sym_f90 in examples\\solver\\source

by converting the real data to complex type with real part are zero. the program complies well without error,

but exit at calling DSS_FACTOR_complex. I do not know why?

!*******************************************************************************

! INTEL CONFIDENTIAL

! The source code contained or described herein and all documents related to

! the source code ("Material") are owned by Intel Corporation or its suppliers

! or licensors. Title to the Material remains with Intel Corporation or its

! suppliers and licensors. The Material contains trade secrets and proprietary

! and confidential information of Intel or its suppliers and licensors. The

! Material is protected by worldwide copyright and trade secret laws and

! treaty provisions. No part of the Material may be used, copied, reproduced,

! modified, published, uploaded, posted, transmitted, distributed or disclosed

! in any way without Intel's prior express written permission.

! No license under any patent, copyright, trade secret or other intellectual

! property right is granted to or conferred upon you by disclosure or delivery

! of the Materials, either expressly, by implication, inducement, estoppel or

! otherwise. Any license under such intellectual property rights must be

! express and approved by Intel in writing.

!*******************************************************************************

! Content : MKL DSS Fortran-90 example

!*******************************************************************************

!--------------------------------------------------------------------------

! Example program for solving a symmetric positive definite system of

! equations.

!--------------------------------------------------------------------------

INCLUDE 'mkl_dss.f90' ! Include the standard DSS "header file."

PROGRAM solver_f90_test

use mkl_dss

IMPLICIT NONE

INTEGER, PARAMETER :: dp = KIND(1.0D0)

INTEGER :: error

INTEGER :: i

INTEGER, PARAMETER :: bufLen = 20

! Define the data arrays and the solution and rhs vectors.

INTEGER, ALLOCATABLE :: columns( : )

INTEGER :: nCols

INTEGER :: nNonZeros

INTEGER :: nRhs

INTEGER :: nRows

REAL(KIND=DP), ALLOCATABLE :: rhst( : )

complex(kind=dp),allocatable :: rhs(:)

INTEGER, ALLOCATABLE :: rowIndex( : )

!REAL(KIND=DP), ALLOCATABLE :: solution( : )

!REAL(KIND=DP), ALLOCATABLE :: values( : )

complex(KIND=DP), ALLOCATABLE :: solution( : )

REAL(KIND=DP), ALLOCATABLE :: valuest( : )

complex(KIND=DP), ALLOCATABLE :: values(:)

TYPE(MKL_DSS_HANDLE) :: handle ! Allocate storage for the solver handle.

REAL(KIND=DP),ALLOCATABLE::statOUt( : )

CHARACTER*15 statIn

INTEGER perm(1)

INTEGER buff(bufLen)

complex(kind=dp) ci

! Set the problem to be solved.

ci = (0_DP,1_DP)

nRows = 5

nCols = 5

nNonZeros = 9

nRhs = 1

perm(1) = 0

ALLOCATE( rowIndex( nRows + 1 ) )

rowIndex = (/ 1, 6, 7, 8, 9, 10 /)

ALLOCATE( columns( nNonZeros ) )

columns = (/ 1, 2, 3, 4, 5, 2, 3, 4, 5 /)

ALLOCATE( valuest( nNonZeros ), values(nNonZeros ))

valuest = (/ 9.0_DP, 1.5_DP, 6.0_DP, 0.75_DP, 3.0_DP, 0.5_DP, 12.0_DP, &

& 0.625_DP, 16.0_DP /)

values = valuest*ci

ALLOCATE( rhs( nRows ),rhst(nrows) )

rhst = (/ 1.0_DP, 2.0_DP, 3.0_DP, 4.0_DP, 5.0_DP /)

rhs = rhst*ci

! Initialize the solver.

error = DSS_CREATE( handle, MKL_DSS_DEFAULTS )

IF (error /= MKL_DSS_SUCCESS) GOTO 999

! Define the non-zero structure of the matrix.

error = DSS_DEFINE_STRUCTURE( handle, MKL_DSS_SYMMETRIC_COMPLEX, rowIndex, nRows, &

& nCols, columns, nNonZeros )

IF (error /= MKL_DSS_SUCCESS) GOTO 999

! Reorder the matrix.

error = DSS_REORDER( handle, MKL_DSS_DEFAULTS, perm )

IF (error /= MKL_DSS_SUCCESS) GOTO 999

! Factor the matrix.

error = DSS_FACTOR_complex( handle, MKL_DSS_DEFAULTS, values )

IF (error /= MKL_DSS_SUCCESS) GOTO 999

! Allocate the solution vector and solve the problem.

ALLOCATE( solution( nRows ) )

error = DSS_SOLVE_complex(handle, MKL_DSS_DEFAULTS, rhs, nRhs, solution )

IF (error /= MKL_DSS_SUCCESS) GOTO 999

! Print Out the determinant of the matrix (no statistics for a diagonal matrix)

IF( nRows .LT. nNonZeros ) THEN

ALLOCATE(statOut( 5 ) )

statIn = 'determinant'

call mkl_cvt_to_null_terminated_str(buff,bufLen,statIn)

error = DSS_STATISTICS(handle, MKL_DSS_DEFAULTS, buff, statOut )

IF (error /= MKL_DSS_SUCCESS) GOTO 999

WRITE(*,"('pow of determinant is '(5F10.3))") ( statOut(1) )

WRITE(*,"('base of determinant is '(5F10.3))") ( statOut(2) )

WRITE(*,"('Determinant is '(5F10.3))") ( (10**statOut(1))*statOut(2) )

END IF

! Deallocate solver storage and various local arrays.

error = DSS_DELETE( handle, MKL_DSS_DEFAULTS )

IF (error /= MKL_DSS_SUCCESS ) GOTO 999

IF ( ALLOCATED( rowIndex) ) DEALLOCATE( rowIndex )

IF ( ALLOCATED( columns ) ) DEALLOCATE( columns )

IF ( ALLOCATED( values ) ) DEALLOCATE( values )

IF ( ALLOCATED( rhs ) ) DEALLOCATE( rhs )

IF ( ALLOCATED( statOut ) ) DEALLOCATE( statOut )

! Print the solution vector, deallocate it and exit

WRITE(*,"('Solution Array: '(5F10.3))") ( solution(i), i = 1, nCols )

IF ( ALLOCATED( solution ) ) DEALLOCATE( solution )

GOTO 1000

! Print an error message and exit

999 WRITE(*,*) "Solver returned error code ", error

STOP 1

1000 CONTINUE

END PROGRAM solver_f90_test

problem in DSS_FACTOR_complex

Sergey_Solovev__Inte — Tue, 11 May 2010 05:11:20 GMT

Example dss_sym_f90handles withpositive definite matrices. If you convert it into complex test with zero real part, the new matrix becomes indefinite. So, you should factorize matrix as indefinite:
error = DSS_FACTOR_complex( handle, MKL_DSS_INDEFINITE, values )

problem in DSS_FACTOR_complex

haijunwu — Wed, 12 May 2010 02:03:27 GMT

Thanks it works.

For general complex matrix, if I just can make sure itsdeterminant is not zero, which option should I set toDSS_FACTOR_complex.

topic problem in DSS_FACTOR_complex in Intel® oneAPI Math Kernel Library

problem in DSS_FACTOR_complex

problem in DSS_FACTOR_complex

problem in DSS_FACTOR_complex