連立一次方程式の計算

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

Syntax

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

Parameters

a
Type: fvalgcli..::..FMATRIX_PTR
MxN 係数行列
x
Type: fvalgcli..::..FMATRIX_PTR
NxNRHS 解行列
b
Type: fvalgcli..::..FMATRIX_PTR
MxNRHS 右辺行列
method
Type: System..::..Int32
解法。下記のいずれかを指定。
  • 0 : QR または LQ 分解を使用した手法により計算
  • 1 : SVD を使用した手法により計算
residuals
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_linsolve2()
        {
            FMATRIX_PTR a = FMATRIX_PTR.Zero;
            FMATRIX_PTR b = FMATRIX_PTR.Zero;
            FMATRIX_PTR x = FMATRIX_PTR.Zero;
            double residual = new double();
            int rank = new int();
            try
            {
                int status = (int)f_err.F_ERR_NONE;
                a = api.fnFIE_mat_aalloc(3,3);
                b = api.fnFIE_mat_aalloc(3,2);
                x = api.fnFIE_mat_aalloc(3,2);

                a[0, 0] = 5; a[0, 1] = 6;
                a[1, 0] = 7; a[1, 1] = 8;

                x[0, 0] = 1; x[0, 1] = 2;
                x[1, 0] = 3; x[1, 1] = 4;

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

                status = api.fnFIE_mat_linsolve2(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
            {
                api.fnFIE_mat_afree(a);
                api.fnFIE_mat_afree(b);
                api.fnFIE_mat_afree(x);
            }
        }
    }
}

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_linsolve2()
        Dim a As FMATRIX_PTR = FMATRIX_PTR.Zero
        Dim b As FMATRIX_PTR = FMATRIX_PTR.Zero
        Dim x As FMATRIX_PTR = FMATRIX_PTR.Zero
        Dim residual As New Double()
        Dim rank As New Integer()
        Try
            Dim status As Integer = CInt(f_err.F_ERR_NONE)
            a = api.fnFIE_mat_aalloc(3, 3)
            b = api.fnFIE_mat_aalloc(3, 2)
            x = api.fnFIE_mat_aalloc(3, 2)

            a(0, 0) = 5
            a(0, 1) = 6
            a(1, 0) = 7
            a(1, 1) = 8

            x(0, 0) = 1
            x(0, 1) = 2
            x(1, 0) = 3
            x(1, 1) = 4

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

            'SVD
            'rcond
            status = api.fnFIE_mat_linsolve2(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
            api.fnFIE_mat_afree(a)
            api.fnFIE_mat_afree(b)
            api.fnFIE_mat_afree(x)
        End Try
    End Sub
End Class

See Also