連立一次方程式の計算

Namespace: fvalgcli
Assembly: fvalgcli (in fvalgcli.dll) Version: 3.1.0.0 (3.1.0.11)

Syntax

C#
public static int fnFIE_mat_linsolve(
	FMATRIX_PTR a,
	FVECTOR_PTR x,
	FVECTOR_PTR b,
	int method,
	ref double residual,
	double rcond,
	ref int rank
)
Visual Basic
Public Shared Function fnFIE_mat_linsolve ( 
	a As FMATRIX_PTR,
	x As FVECTOR_PTR,
	b As FVECTOR_PTR,
	method As Integer,
	ByRef residual As Double,
	rcond As Double,
	ByRef rank As Integer
) As Integer

Parameters

a
Type: fvalgcli..::..FMATRIX_PTR
MxN 係数行列
x
Type: fvalgcli..::..FVECTOR_PTR
解ベクトル。x の次元はNでなければいけません。
b
Type: fvalgcli..::..FVECTOR_PTR
右辺ベクトル。b の次元はMでなければいけません。
method
Type: System..::..Int32
解法。下記のいずれかを指定。
  • 0 : QR または LQ 分解を使用した手法により計算
  • 1 : SVD を使用した手法により計算
residual
Type: System..::..Double%
誤差二乗和。 a.row ≦ a.col の場合、本パラメータは使用されません。 不要な場合は IntPtr.Zero を指定してください。
rcond
Type: System..::..Double
ゼロ判定閾値 method != 1 の場合本パラメータは無視されます。 ランク計算時に 特異値 sigma[i] ≦ rcond*sigma[0] の値をゼロと見なします。 rcond＜0.0 の値が指定されたときは、max(sigma) * eps の値を使用します。
rank
Type: System..::..Int32%
行列 a のランク。

Return Value

Type: Int32
以下のエラーコードを返します。

エラーコード：
f_err内容
F_ERR_NONE正常終了
F_ERR_NOMEMORYメモリ不足エラー
F_ERR_CALC_IMPOSSIBLE計算不能エラー
F_ERR_INVALID_PARAMパラメータ異常
F_ERR_NO_LICENCEライセンスエラー、または未初期化エラー

Examples

C# Copy imageCopy
//    $Revision: 1.1 $

using System;
using System.Collections.Generic;
using System.Text;
using fvalgcli;

namespace TC
{
    public partial class FIE
    {
        [FvPluginExecute]
        public void fnFIE_mat_linsolve()
        {
            FMATRIX_PTR a = FMATRIX_PTR.Zero;
            FVECTOR_PTR b = FVECTOR_PTR.Zero;
            FVECTOR_PTR x = FVECTOR_PTR.Zero;
            double residual = 0;
            int rank = 0;
            try
            {
                int status = (int)f_err.F_ERR_NONE;
                a = api.fnFIE_mat_aalloc(5, 3);
                b = api.fnFIE_mat_valloc(5);
                x = api.fnFIE_mat_valloc(3);

                // 値を格納.
                a[0, 0] = 4; a[0, 1] = 5; a[0, 2] = 8;
                a[1, 0] = 3; a[1, 1] = 2; a[1, 2] = 1;
                a[2, 0] = 4; a[2, 1] = 5; a[2, 2] = 6;
                a[3, 0] = 1; a[3, 1] = 3; a[3, 2] = -1;
                a[4, 0] = 9; a[4, 1] = 6; a[4, 2] = 2;

                // 値を格納.
                b[0] = 38; b[1] = 10; b[2] = 32; b[3] = 4; b[4] = 27;

                Assert.IsTrue(a != FMATRIX_PTR.Zero, "エラーが発生しました。");

                status = api.fnFIE_mat_linsolve(a, x, b, 1/*SVD*/,  ref residual, 1e-10/*rcond*/, ref rank);
                Assert.IsTrue(status == (int)f_err.F_ERR_NONE, "エラーが発生しました。({0})", (f_err)status);

            }
            catch (System.Exception ex)
            {
                Assert.Fail(ex.Message);
            }
            finally
            {
                a.Dispose();
                b.Dispose();
                x.Dispose();
            }
        }
    }
}

Visual Basic Copy imageCopy
'    $Revision: 1.1 $

Imports System.Collections.Generic
Imports System.Text
Imports fvalgcli

Public Partial Class FIE
    <FvPluginExecute> _
    Public Sub fnFIE_mat_linsolve()
        Dim a As FMATRIX_PTR = FMATRIX_PTR.Zero
        Dim b As FVECTOR_PTR = FVECTOR_PTR.Zero
        Dim x As FVECTOR_PTR = FVECTOR_PTR.Zero
        Dim residual As Double = 0
        Dim rank As Integer = 0
        Try
            Dim status As Integer = CInt(f_err.F_ERR_NONE)
            a = api.fnFIE_mat_aalloc(5, 3)
            b = api.fnFIE_mat_valloc(5)
            x = api.fnFIE_mat_valloc(3)

            ' 値を格納.
            a(0, 0) = 4
            a(0, 1) = 5
            a(0, 2) = 8
            a(1, 0) = 3
            a(1, 1) = 2
            a(1, 2) = 1
            a(2, 0) = 4
            a(2, 1) = 5
            a(2, 2) = 6
            a(3, 0) = 1
            a(3, 1) = 3
            a(3, 2) = -1
            a(4, 0) = 9
            a(4, 1) = 6
            a(4, 2) = 2

            ' 値を格納.
            b(0) = 38
            b(1) = 10
            b(2) = 32
            b(3) = 4
            b(4) = 27

            Assert.IsTrue(a <> FMATRIX_PTR.Zero, "エラーが発生しました。")

            'SVD
            'rcond
            status = api.fnFIE_mat_linsolve(a, x, b, 1, residual, 1E-10, _
                rank)

            Assert.IsTrue(status = CInt(f_err.F_ERR_NONE), "エラーが発生しました。({0})", CType(status, f_err))
        Catch ex As System.Exception
            Assert.Fail(ex.Message)
        Finally
            a.Dispose()
            b.Dispose()
            x.Dispose()
        End Try
    End Sub
End Class

See Also