pcg_example.m File Reference

Functions
function	pcg_example (dim) if nargin< 1 dim
	b (1:n)
end	A (idx, idx)
	disp ('L2 norm|x_matlab-x_paralution|is') disp(norm(err))

Variables
	nnz = n*5-dim*4
	A = sparse(n, n)
for	i
end if	j
end end end	x_ref = pcg(A, b', 1e-6, 500)
	x_paralution = paralution_pcg(A, b', 1e-6, 500)
	err = x_ref - x_paralution

Function Documentation

end A	(	idx	,
		idx
	)

b

(

1:n

)

Referenced by paralution::SIRA< OperatorTypeH, VectorTypeH, ValueTypeH, OperatorTypeL, VectorTypeL, ValueTypeL >::Givens(), main(), and mexFunction().

disp

(

'L2 norm|x_matlab-x_paralution|is'

)

function pcg_example

(

dim

)

Variable Documentation

A = sparse(n, n)

Referenced by paralution::HostMatrixDENSE< ValueType >::MatMatMult(), paralution::HostMatrixCSR< ValueType >::MatMatMult(), mexFunction(), paralution::HostMatrixCSR< ValueType >::NumericMatMatMult(), and paralution::HostMatrixCSR< ValueType >::SymbolicMatMatMult().

err = x_ref - x_paralution

Referenced by paralution::OCLAcceleratorMatrixCSR< ValueType >::AddScalar(), paralution::LocalMatrix< ValueType >::AddScalar(), paralution::OCLAcceleratorMatrixCSR< ValueType >::AddScalarDiagonal(), paralution::LocalMatrix< ValueType >::AddScalarDiagonal(), paralution::OCLAcceleratorMatrixCSR< ValueType >::AddScalarOffDiagonal(), paralution::LocalMatrix< ValueType >::AddScalarOffDiagonal(), paralution::OCLAcceleratorVector< ValueType >::AddScale(), paralution::allocate_ocl(), paralution::OCLAcceleratorVector< ValueType >::Amax(), paralution::LocalMatrix< ValueType >::AMGAggregate(), paralution::LocalMatrix< ValueType >::AMGAggregation(), paralution::LocalMatrix< ValueType >::AMGConnect(), paralution::LocalMatrix< ValueType >::AMGSmoothedAggregation(), paralution::OCLAcceleratorMatrixMCSR< ValueType >::Apply(), paralution::OCLAcceleratorMatrixBCSR< ValueType >::Apply(), paralution::OCLAcceleratorMatrixDENSE< ValueType >::Apply(), paralution::OCLAcceleratorMatrixCOO< ValueType >::Apply(), paralution::OCLAcceleratorMatrixDIA< ValueType >::Apply(), paralution::OCLAcceleratorMatrixELL< ValueType >::Apply(), paralution::OCLAcceleratorMatrixHYB< ValueType >::Apply(), paralution::OCLAcceleratorMatrixCSR< ValueType >::Apply(), paralution::OCLAcceleratorMatrixMCSR< ValueType >::ApplyAdd(), paralution::OCLAcceleratorMatrixCOO< ValueType >::ApplyAdd(), paralution::OCLAcceleratorMatrixDIA< ValueType >::ApplyAdd(), paralution::OCLAcceleratorMatrixELL< ValueType >::ApplyAdd(), paralution::OCLAcceleratorMatrixHYB< ValueType >::ApplyAdd(), paralution::OCLAcceleratorMatrixCSR< ValueType >::ApplyAdd(), paralution::LocalMatrix< ValueType >::Assemble(), paralution::LocalMatrix< ValueType >::AssembleUpdate(), paralution::OCLAcceleratorVector< ValueType >::Asum(), paralution::LocalMatrix< ValueType >::CMK(), paralution::OCLAcceleratorMatrixCSR< ValueType >::Compress(), paralution::LocalMatrix< ValueType >::Compress(), paralution::OCLAcceleratorMatrixELL< ValueType >::ConvertFrom(), paralution::OCLAcceleratorVector< ValueType >::CopyFrom(), paralution::OCLAcceleratorVector< ValueType >::CopyFromDouble(), paralution::OCLAcceleratorVector< ValueType >::CopyFromFloat(), paralution::OCLAcceleratorVector< ValueType >::CopyFromPermute(), paralution::OCLAcceleratorVector< ValueType >::CopyFromPermuteBackward(), paralution::LocalMatrix< ValueType >::CreateFromMap(), paralution::OCLAcceleratorMatrixCSR< ValueType >::DiagonalMatrixMultL(), paralution::LocalMatrix< ValueType >::DiagonalMatrixMultL(), paralution::OCLAcceleratorMatrixCSR< ValueType >::DiagonalMatrixMultR(), paralution::LocalMatrix< ValueType >::DiagonalMatrixMultR(), paralution::OCLAcceleratorVector< ValueType >::Dot(), paralution::OCLAcceleratorVector< ValueType >::DotNonConj(), paralution::OCLAcceleratorMatrixDENSE< ValueType >::ExtractColumnVector(), paralution::OCLAcceleratorMatrixCSR< ValueType >::ExtractColumnVector(), paralution::LocalMatrix< ValueType >::ExtractColumnVector(), paralution::OCLAcceleratorMatrixCSR< ValueType >::ExtractDiagonal(), paralution::LocalMatrix< ValueType >::ExtractDiagonal(), paralution::OCLAcceleratorMatrixCSR< ValueType >::ExtractInverseDiagonal(), paralution::LocalMatrix< ValueType >::ExtractInverseDiagonal(), paralution::OCLAcceleratorMatrixCSR< ValueType >::ExtractL(), paralution::LocalMatrix< ValueType >::ExtractL(), paralution::OCLAcceleratorMatrixCSR< ValueType >::ExtractLDiagonal(), paralution::OCLAcceleratorMatrixDENSE< ValueType >::ExtractRowVector(), paralution::LocalMatrix< ValueType >::ExtractRowVector(), paralution::OCLAcceleratorMatrixCSR< ValueType >::ExtractSubMatrix(), paralution::LocalMatrix< ValueType >::ExtractSubMatrix(), paralution::OCLAcceleratorMatrixCSR< ValueType >::ExtractU(), paralution::LocalMatrix< ValueType >::ExtractU(), paralution::OCLAcceleratorMatrixCSR< ValueType >::ExtractUDiagonal(), paralution::free_ocl(), paralution::LocalMatrix< ValueType >::FSAI(), paralution::LocalMatrix< ValueType >::Gershgorin(), paralution::LocalMatrix< ValueType >::Householder(), paralution::LocalMatrix< ValueType >::ICFactorize(), paralution::LocalMatrix< ValueType >::ILU0Factorize(), paralution::LocalMatrix< ValueType >::ILUpFactorize(), paralution::LocalMatrix< ValueType >::ILUTFactorize(), paralution::LocalMatrix< ValueType >::Invert(), paralution::LocalMatrix< ValueType >::LLSolve(), paralution::LocalMatrix< ValueType >::LSolve(), paralution::LocalMatrix< ValueType >::LUFactorize(), paralution::LocalMatrix< ValueType >::LUSolve(), paralution::OCLAcceleratorMatrixCSR< ValueType >::MatrixAdd(), paralution::LocalMatrix< ValueType >::MatrixAdd(), paralution::LocalMatrix< ValueType >::MatrixMult(), paralution::LocalMatrix< ValueType >::MaximalIndependentSet(), paralution::LocalMatrix< ValueType >::MultiColoring(), paralution::OCLAcceleratorVector< ValueType >::Norm(), paralution::ocl_dev2dev(), paralution::ocl_dev2host(), paralution::ocl_host2dev(), paralution::ocl_set_to(), paralution::paralution_get_device_ocl(), paralution::paralution_get_platform_ocl(), paralution::paralution_info_ocl(), paralution::paralution_init_ocl(), paralution::paralution_set_kernels_ocl(), paralution::paralution_stop_ocl(), paralution::OCLAcceleratorMatrixCSR< ValueType >::Permute(), paralution::OCLAcceleratorVector< ValueType >::Permute(), paralution::LocalMatrix< ValueType >::Permute(), paralution::OCLAcceleratorVector< ValueType >::PermuteBackward(), paralution::LocalMatrix< ValueType >::PermuteBackward(), paralution::OCLAcceleratorVector< ValueType >::PointWiseMult(), paralution::OCLAcceleratorVector< ValueType >::Power(), paralution::LocalVector< ValueType >::Prolongation(), paralution::LocalMatrix< ValueType >::QRDecompose(), paralution::LocalMatrix< ValueType >::QRSolve(), paralution::LocalMatrix< ValueType >::RCMK(), paralution::LocalMatrix< ValueType >::ReadFileCSR(), paralution::LocalMatrix< ValueType >::ReadFileMTX(), paralution::OCLAcceleratorVector< ValueType >::Reduce(), paralution::OCLAcceleratorMatrixDENSE< ValueType >::ReplaceColumnVector(), paralution::OCLAcceleratorMatrixCSR< ValueType >::ReplaceColumnVector(), paralution::LocalMatrix< ValueType >::ReplaceColumnVector(), paralution::OCLAcceleratorMatrixDENSE< ValueType >::ReplaceRowVector(), paralution::LocalMatrix< ValueType >::ReplaceRowVector(), paralution::LocalVector< ValueType >::Restriction(), paralution::OCLAcceleratorMatrixCSR< ValueType >::Scale(), paralution::OCLAcceleratorVector< ValueType >::Scale(), paralution::LocalMatrix< ValueType >::Scale(), paralution::OCLAcceleratorVector< ValueType >::ScaleAdd(), paralution::OCLAcceleratorVector< ValueType >::ScaleAdd2(), paralution::OCLAcceleratorVector< ValueType >::ScaleAddScale(), paralution::OCLAcceleratorMatrixCSR< ValueType >::ScaleDiagonal(), paralution::LocalMatrix< ValueType >::ScaleDiagonal(), paralution::OCLAcceleratorMatrixCSR< ValueType >::ScaleOffDiagonal(), paralution::LocalMatrix< ValueType >::ScaleOffDiagonal(), paralution::LocalMatrix< ValueType >::Sort(), paralution::LocalMatrix< ValueType >::SPAI(), paralution::LocalMatrix< ValueType >::SymbolicPower(), paralution::LocalMatrix< ValueType >::Transpose(), paralution::LocalMatrix< ValueType >::USolve(), paralution::LocalMatrix< ValueType >::WriteFileCSR(), paralution::LocalMatrix< ValueType >::WriteFileMTX(), paralution::LocalMatrix< ValueType >::ZeroBlockPermutation(), and paralution::LocalMatrix< ValueType >::Zeros().

end if i

Initial value:

=0:dim-1
    for j=0:dim-1

      idx = i*dim+j+1

Referenced by paralution::HostVector< ValueType >::AddScale(), paralution::HostVector< ValueType >::Amax(), paralution::HostVector< ValueType >::Assemble(), paralution::HostVector< ValueType >::Asum(), paralution::HostVector< ValueType >::CopyFrom(), paralution::HostVector< ValueType >::CopyFromPermute(), paralution::HostVector< ValueType >::CopyFromPermuteBackward(), paralution::HostVector< ValueType >::Dot(), paralution::HostVector< ValueType >::DotNonConj(), export_dealii_vector(), getwd(), import_dealii_matrix(), import_dealii_vector(), lowerc(), main(), mexFunction(), paralution::HostVector< ValueType >::Norm(), paralution_fortran_solve_coo(), paralution_fortran_solve_csr(), paralution_solve(), paralution::HostVector< ValueType >::PartialSum(), paralution::HostVector< ValueType >::Permute(), paralution::HostVector< ValueType >::PermuteBackward(), paralution::HostVector< ValueType >::PointWiseMult(), paralution::HostVector< ValueType >::Power(), paralution::HostVector< ValueType >::Prolongation(), paralution::HostVector< ValueType >::Reduce(), paralution::HostVector< ValueType >::Restriction(), paralution::HostVector< ValueType >::Scale(), paralution::HostVector< ValueType >::ScaleAdd(), paralution::HostVector< ValueType >::ScaleAdd2(), and paralution::HostVector< ValueType >::ScaleAddScale().

if j

Initial value:

= 0
        A(idx,idx-1) = -1.0

Referenced by paralution::HostMatrixCSR< ValueType >::AMGAggregate(), paralution::HostMatrixCSR< ValueType >::AMGAggregation(), paralution::HostMatrixCSR< ValueType >::AMGConnect(), paralution::HostMatrixCSR< ValueType >::AMGSmoothedAggregation(), paralution::HostStencilLaplace2D< ValueType >::Apply(), paralution::HostMatrixDIA< ValueType >::Apply(), paralution::HostStencilLaplace2D< ValueType >::ApplyAdd(), paralution::HostMatrixDIA< ValueType >::ApplyAdd(), paralution::FGMRES< OperatorType, VectorType, ValueType >::BackSubstitute_(), paralution::GMRES< OperatorType, VectorType, ValueType >::BackSubstitute_(), paralution::IDR< OperatorType, VectorType, ValueType >::Build(), paralution::AS< OperatorType, VectorType, ValueType >::Build(), paralution::BlockPreconditioner< OperatorType, VectorType, ValueType >::Build(), paralution::MultiColored< OperatorType, VectorType, ValueType >::Clear(), paralution::BlockPreconditioner< OperatorType, VectorType, ValueType >::Clear(), paralution::HostMatrixCSR< ValueType >::CMK(), paralution::HostMatrixCSR< ValueType >::Compress(), paralution::coo_to_csr(), paralution::HostMatrixMCSR< ValueType >::CopyFrom(), paralution::HostMatrixDENSE< ValueType >::CopyFrom(), paralution::HostMatrixDIA< ValueType >::CopyFrom(), paralution::HostMatrixCOO< ValueType >::CopyFrom(), paralution::HostMatrixCSR< ValueType >::CopyFrom(), paralution::HostMatrixCOO< ValueType >::CopyFromCOO(), paralution::HostMatrixCSR< ValueType >::CopyFromCSR(), paralution::HostMatrixCOO< ValueType >::CopyToCOO(), paralution::HostMatrixCSR< ValueType >::CopyToCSR(), paralution::csr_to_coo(), paralution::csr_to_dense(), paralution::csr_to_dia(), paralution::csr_to_ell(), paralution::csr_to_hyb(), paralution::csr_to_mcsr(), paralution::MultiColored< OperatorType, VectorType, ValueType >::Decompose_(), paralution::dense_to_csr(), paralution::dia_to_csr(), paralution::SIRA< OperatorTypeH, VectorTypeH, ValueTypeH, OperatorTypeL, VectorTypeL, ValueTypeL >::Eigpair(), paralution::LocalMatrix< ValueType >::ExtractSubMatrices(), paralution::HostMatrixCSR< ValueType >::FSAI(), paralution::HostMatrixCSR< ValueType >::ICFactorize(), paralution::HostMatrixCSR< ValueType >::ILU0Factorize(), paralution::HostMatrixCSR< ValueType >::ILUpFactorizeNumeric(), paralution::HostMatrixCSR< ValueType >::ILUTFactorize(), paralution::HostMatrixDENSE< ValueType >::Invert(), paralution::kernel_dia_convert(), paralution::kernel_dia_diag_map(), paralution::kernel_permute_cols(), paralution::SIRA< OperatorTypeH, VectorTypeH, ValueTypeH, OperatorTypeL, VectorTypeL, ValueTypeL >::Locking(), paralution::HostMatrixDENSE< ValueType >::LUFactorize(), paralution::HostMatrixDENSE< ValueType >::LUSolve(), main(), paralution::DPCG< OperatorType, VectorType, ValueType >::MakeZ_COO(), paralution::DPCG< OperatorType, VectorType, ValueType >::MakeZ_CSR(), paralution::HostMatrixDENSE< ValueType >::MatMatMult(), paralution::HostMatrixCSR< ValueType >::MatMatMult(), paralution::HostMatrixCSR< ValueType >::MatrixAdd(), paralution::mcsr_to_csr(), paralution::BlockPreconditioner< OperatorType, VectorType, ValueType >::MoveToAcceleratorLocalData_(), paralution::MultiColored< OperatorType, VectorType, ValueType >::MoveToAcceleratorLocalData_(), paralution::BlockPreconditioner< OperatorType, VectorType, ValueType >::MoveToHostLocalData_(), paralution::MultiColored< OperatorType, VectorType, ValueType >::MoveToHostLocalData_(), paralution::HostMatrixCSR< ValueType >::NumericMatMatMult(), paralution::HostMatrixCSR< ValueType >::Permute(), paralution::HostVector< ValueType >::Power(), paralution::HostMatrixDENSE< ValueType >::QRSolve(), paralution::MultiColoredSGS< OperatorType, VectorType, ValueType >::ReBuildNumeric(), paralution::MultiColoredILU< OperatorType, VectorType, ValueType >::ReBuildNumeric(), paralution::HostMatrixCSR< ValueType >::ReplaceColumnVector(), paralution::HostMatrixCSR< ValueType >::ReplaceRowVector(), paralution::BlockPreconditioner< OperatorType, VectorType, ValueType >::Solve(), paralution::MultiColoredSGS< OperatorType, VectorType, ValueType >::SolveL_(), paralution::MultiColoredILU< OperatorType, VectorType, ValueType >::SolveL_(), paralution::GMRES< OperatorType, VectorType, ValueType >::SolveNonPrecond_(), paralution::FGMRES< OperatorType, VectorType, ValueType >::SolveNonPrecond_(), paralution::IDR< OperatorType, VectorType, ValueType >::SolveNonPrecond_(), paralution::FGMRES< OperatorType, VectorType, ValueType >::SolvePrecond_(), paralution::GMRES< OperatorType, VectorType, ValueType >::SolvePrecond_(), paralution::IDR< OperatorType, VectorType, ValueType >::SolvePrecond_(), paralution::MultiColoredSGS< OperatorType, VectorType, ValueType >::SolveR_(), paralution::MultiColoredILU< OperatorType, VectorType, ValueType >::SolveR_(), paralution::MultiColoredGS< OperatorType, VectorType, ValueType >::SolveR_(), paralution::HostMatrixCSR< ValueType >::Sort(), paralution::HostMatrixCSR< ValueType >::SPAI(), paralution::spmv_add_dia(), paralution::spmv_dia(), and paralution::HostMatrixCSR< ValueType >::SymbolicMatMatMult().

nnz = n*5-dim*4

Referenced by paralution::OCLAcceleratorMatrixCSR< ValueType >::AddScalar(), paralution::HostMatrixBCSR< ValueType >::AllocateBCSR(), paralution::HostMatrixCOO< ValueType >::AllocateCOO(), paralution::MICAcceleratorMatrixCOO< ValueType >::AllocateCOO(), paralution::OCLAcceleratorMatrixCOO< ValueType >::AllocateCOO(), paralution::HostMatrixCSR< ValueType >::AllocateCSR(), paralution::MICAcceleratorMatrixCSR< ValueType >::AllocateCSR(), paralution::OCLAcceleratorMatrixCSR< ValueType >::AllocateCSR(), paralution::MICAcceleratorMatrixDIA< ValueType >::AllocateDIA(), paralution::HostMatrixDIA< ValueType >::AllocateDIA(), paralution::OCLAcceleratorMatrixDIA< ValueType >::AllocateDIA(), paralution::HostMatrixELL< ValueType >::AllocateELL(), paralution::MICAcceleratorMatrixELL< ValueType >::AllocateELL(), paralution::OCLAcceleratorMatrixELL< ValueType >::AllocateELL(), paralution::HostMatrixMCSR< ValueType >::AllocateMCSR(), paralution::MICAcceleratorMatrixMCSR< ValueType >::AllocateMCSR(), paralution::OCLAcceleratorMatrixMCSR< ValueType >::AllocateMCSR(), paralution::HostMatrixCSR< ValueType >::AMGSmoothedAggregation(), paralution::OCLAcceleratorMatrixCOO< ValueType >::Apply(), paralution::OCLAcceleratorMatrixHYB< ValueType >::Apply(), paralution::OCLAcceleratorMatrixCOO< ValueType >::ApplyAdd(), paralution::OCLAcceleratorMatrixHYB< ValueType >::ApplyAdd(), paralution::HostMatrixCSR< ValueType >::Assemble(), paralution::HostMatrixBCSR< ValueType >::ConvertFrom(), paralution::HostMatrixDIA< ValueType >::ConvertFrom(), paralution::HostMatrixHYB< ValueType >::ConvertFrom(), paralution::HostMatrixELL< ValueType >::ConvertFrom(), paralution::HostMatrixCSR< ValueType >::ConvertFrom(), paralution::coo_to_csr(), paralution::HostMatrixELL< ValueType >::CopyFrom(), paralution::HostMatrixCSR< ValueType >::CreateFromMap(), paralution::csr_to_coo(), paralution::dense_to_csr(), paralution::HostMatrixCSR< ValueType >::ExtractL(), paralution::HostMatrixCSR< ValueType >::ExtractLDiagonal(), paralution::HostMatrixCSR< ValueType >::ExtractU(), paralution::HostMatrixCSR< ValueType >::ExtractUDiagonal(), paralution::HostMatrixCSR< ValueType >::FSAI(), paralution::HostMatrixCSR< ValueType >::ILUpFactorizeNumeric(), paralution::HostMatrixCSR< ValueType >::ILUTFactorize(), import_dealii_matrix(), paralution::kernel_coo_permute(), paralution::kernel_csr_extract_submatrix_row_nnz(), main(), mexFunction(), paralution_fortran_solve_coo(), paralution_fortran_solve_csr(), paralution_solve(), paralution::read_coordinate_stream(), paralution::HostMatrixCSR< ValueType >::ReadFileCSR(), paralution::HostMatrixCOO< ValueType >::ReadFileMTX(), paralution::OCLAcceleratorMatrixCSR< ValueType >::ReplaceColumnVector(), paralution::HostMatrixCSR< ValueType >::ReplaceColumnVector(), paralution::HostMatrixCSR< ValueType >::ReplaceRowVector(), paralution::OCLAcceleratorMatrixCSR< ValueType >::Scale(), paralution::HostMatrixCOO< ValueType >::SetDataPtrCOO(), paralution::HostMatrixCSR< ValueType >::SetDataPtrCSR(), paralution::HostMatrixDIA< ValueType >::SetDataPtrDIA(), paralution::HostMatrixELL< ValueType >::SetDataPtrELL(), paralution::HostMatrixCSR< ValueType >::SPAI(), paralution::spmv_add_coo(), paralution::spmv_coo(), paralution::LocalMatrix< ValueType >::UpdateValuesCSR(), and paralution::write_matrix_mtx().

x_paralution = paralution_pcg(A, b', 1e-6, 500)

end end end x_ref = pcg(A, b', 1e-6, 500)