/// /// This file is part of ILNumerics Community Edition. /// /// ILNumerics Community Edition - high performance computing for applications. /// Copyright (C) 2006 - 2012 Haymo Kutschbach, http://ilnumerics.net /// /// ILNumerics Community Edition is free software: you can redistribute it and/or modify /// it under the terms of the GNU General Public License version 3 as published by /// the Free Software Foundation. /// /// ILNumerics Community Edition is distributed in the hope that it will be useful, /// but WITHOUT ANY WARRANTY; without even the implied warranty of /// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the /// GNU General Public License for more details. /// /// You should have received a copy of the GNU General Public License /// along with ILNumerics Community Edition. See the file License.txt in the root /// of your distribution package. If not, see . /// /// In addition this software uses the following components and/or licenses: /// /// ================================================================================= /// The Open Toolkit Library License /// /// Copyright (c) 2006 - 2009 the Open Toolkit library. /// /// Permission is hereby granted, free of charge, to any person obtaining a copy /// of this software and associated documentation files (the "Software"), to deal /// in the Software without restriction, including without limitation the rights to /// use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of /// the Software, and to permit persons to whom the Software is furnished to do /// so, subject to the following conditions: /// /// The above copyright notice and this permission notice shall be included in all /// copies or substantial portions of the Software. /// /// ================================================================================= /// #pragma warning disable 1591 using System; using System.Collections.Generic; using System.Text; using ILNumerics.Native; using System.Runtime.InteropServices; using ILNumerics.Exceptions; using System.Security; using ILNumerics.Misc; namespace ILNumerics.Native { ///

/// Generic LAPACK implementation, unsupported processor types ///

public class ILLapackGeneric : IILLapack { #region DLL INCLUDES [DllImport("lapack_gen", EntryPoint = "ilaenv"),SuppressUnmanagedCodeSecurity] private static extern int lapack_ilaenv(ref int ispec, ref string name, ref string opts, ref int n1, ref int n2, ref int n3, ref int n4); /////////////////////////// DOUBLE LAPACK ///////////////////////////////// [DllImport("lapack_gen", EntryPoint = "dgemm"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dgemm(ref char TransA, ref char TransB, ref int M,ref int N, ref int K, ref double alpha, IntPtr A, ref int lda, IntPtr B, ref int ldb, ref double beta, double[] C, ref int ldc); [DllImport("lapack_gen", EntryPoint = "sgemm"),SuppressUnmanagedCodeSecurity] private static extern void lapack_sgemm(ref char TransA, ref char TransB, ref int M,ref int N, ref int K, ref float alpha, IntPtr A,ref int lda, IntPtr B, ref int ldb, ref float beta,float[] C, ref int ldc); [DllImport("lapack_gen", EntryPoint = "cgemm"),SuppressUnmanagedCodeSecurity] private static extern void lapack_cgemm(ref char TransA, ref char TransB, ref int M,ref int N, ref int K, ref fcomplex alpha, IntPtr A, ref int lda, IntPtr B, ref int ldb, ref fcomplex beta,[In,Out] fcomplex [] C, ref int ldc); [DllImport("lapack_gen", EntryPoint = "zgemm"),SuppressUnmanagedCodeSecurity] private static extern void lapack_zgemm(ref char TransA, ref char TransB, ref int M, ref int N, ref int K, ref complex alpha, IntPtr A, ref int lda, IntPtr B, ref int ldb, ref complex beta,[In,Out] complex [] C, ref int ldc); [DllImport("lapack_gen", EntryPoint = "dgesdd"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dgesdd(ref char jobz, ref int m, ref int n,double[] a, ref int lda, double[] s,double[] u, ref int ldu, double[] vt,ref int ldvt, double[] work, ref int lwork,int[] iwork, ref int info); [DllImport("lapack_gen", EntryPoint = "sgesdd"),SuppressUnmanagedCodeSecurity] private static extern void lapack_sgesdd(ref char jobz, ref int m, ref int n,float[] a, ref int lda, float[] s,float[] u, ref int ldu, float[] vt, ref int ldvt, float[] work, ref int lwork,int[] iwork, ref int info); [DllImport("lapack_gen", EntryPoint = "cgesdd"),SuppressUnmanagedCodeSecurity] private static extern void lapack_cgesdd(ref char jobz, ref int m, ref int n,[In, Out] fcomplex[] a, ref int lda, float[] s,[In, Out] fcomplex[] u, ref int ldu, [In, Out] fcomplex[] vt, ref int ldvt, [In, Out] fcomplex[] work, ref int lwork,[In,Out] float[] rwork, int[] iwork, ref int info); [DllImport("lapack_gen", EntryPoint = "zgesdd"),SuppressUnmanagedCodeSecurity] private static extern void lapack_zgesdd(ref char jobz, ref int m, ref int n,[In, Out] complex[] a, ref int lda, double[] s,[In, Out] complex[] u, ref int ldu, [In, Out] complex[] vt,ref int ldvt, [In, Out] complex[] work, ref int lwork,[In,Out] double[] rwork, int[] iwork, ref int info); [DllImport("lapack_gen", EntryPoint = "dgesvd"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dgesvd(ref char jobu, ref char jobvt, ref int m, ref int n,double[] a, ref int lda, double[] s,double[] u, ref int ldu, double[] vt,ref int ldvt, double[] work, ref int lwork,int[] iwork,ref int info); [DllImport("lapack_gen", EntryPoint = "sgesvd"),SuppressUnmanagedCodeSecurity] private static extern void lapack_sgesvd(ref char jobu, ref char jobvt, ref int m, ref int n, float[] a, ref int lda, float[] s,float[] u, ref int ldu, float[] vt,ref int ldvt, float[] work, ref int lwork,int[] iwork, ref int info); [DllImport("lapack_gen", EntryPoint = "cgesvd"),SuppressUnmanagedCodeSecurity] private static extern void lapack_cgesvd(ref char jobu, ref char jobvt, ref int m, ref int n,[In, Out] fcomplex[] a, ref int lda, float[] s,[In, Out] fcomplex[] u, ref int ldu,[In,Out] fcomplex[] vt,ref int ldvt, [In, Out] fcomplex[] work, ref int lwork,int[] iwork,ref int info); [DllImport("lapack_gen", EntryPoint = "zgesvd"),SuppressUnmanagedCodeSecurity] private static extern void lapack_zgesvd(ref char jobu, ref char jobvt, ref int m, ref int n,[In, Out] complex[] a, ref int lda, double[] s,[In, Out] complex[] u, ref int ldu, [In, Out] complex[] vt,ref int ldvt, [In, Out] complex[] work, ref int lwork,int[] iwork,ref int info); [DllImport("lapack_gen", EntryPoint = "dpotrf"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dpotrf (ref char uplo, ref int n, double [] A, ref int lda, ref int info); [DllImport("lapack_gen", EntryPoint = "spotrf"),SuppressUnmanagedCodeSecurity] private static extern void lapack_spotrf (ref char uplo, ref int n, float [] A, ref int lda, ref int info); [DllImport("lapack_gen", EntryPoint = "cpotrf"),SuppressUnmanagedCodeSecurity] private static extern void lapack_cpotrf (ref char uplo, ref int n,[In,Out] fcomplex [] A, ref int lda, ref int info); [DllImport("lapack_gen", EntryPoint = "zpotrf"),SuppressUnmanagedCodeSecurity] private static extern void lapack_zpotrf (ref char uplo, ref int n, [In,Out] complex [] A, ref int lda, ref int info); [DllImport("lapack_gen", EntryPoint = "dpotri"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dpotri (ref char uplo,ref int n, double [] A, ref int lda, ref int info); [DllImport("lapack_gen", EntryPoint = "spotri"),SuppressUnmanagedCodeSecurity] private static extern void lapack_spotri (ref char uplo, ref int n, float [] A, ref int lda, ref int info); [DllImport("lapack_gen", EntryPoint = "cpotri"),SuppressUnmanagedCodeSecurity] private static extern void lapack_cpotri (ref char uplo, ref int n,[In,Out] fcomplex [] A, ref int lda,ref int info); [DllImport("lapack_gen", EntryPoint = "zpotri"),SuppressUnmanagedCodeSecurity] private static extern void lapack_zpotri (ref char uplo, ref int n, [In,Out] complex [] A, ref int lda,ref int info); [DllImport("lapack_gen", EntryPoint = "dpotrs"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dpotrs (ref char uplo,ref int n, ref int NRHS, double [] A, ref int lda, double[] B, ref int ldb, ref int info); [DllImport("lapack_gen", EntryPoint = "spotrs"),SuppressUnmanagedCodeSecurity] private static extern void lapack_spotrs (ref char uplo,ref int n, ref int NRHS, float [] A, ref int lda, float[] B, ref int ldb, ref int info); [DllImport("lapack_gen", EntryPoint = "cpotrs"),SuppressUnmanagedCodeSecurity] private static extern void lapack_cpotrs (ref char uplo,ref int n, ref int NRHS, [In,Out] fcomplex [] A, ref int lda, [In,Out] fcomplex[] B, ref int ldb, ref int info); [DllImport("lapack_gen", EntryPoint = "zpotrs"),SuppressUnmanagedCodeSecurity] private static extern void lapack_zpotrs (ref char uplo,ref int n, ref int NRHS, [In,Out] complex [] A, ref int lda, [In,Out] complex[] B, ref int ldb, ref int info); [DllImport("lapack_gen", EntryPoint = "dgetrf"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dgetrf (ref int M, ref int N, double [] A, ref int LDA, int[] IPIV, ref int info); [DllImport("lapack_gen", EntryPoint = "sgetrf"),SuppressUnmanagedCodeSecurity] private static extern void lapack_sgetrf (ref int M, ref int N, float[] A, ref int LDA, int [] IPIV, ref int info); [DllImport("lapack_gen", EntryPoint = "cgetrf"),SuppressUnmanagedCodeSecurity] private static extern void lapack_cgetrf (ref int M, ref int N, [In,Out] fcomplex [] A, ref int LDA, int [] IPIV, ref int info); [DllImport("lapack_gen", EntryPoint = "zgetrf"),SuppressUnmanagedCodeSecurity] private static extern void lapack_zgetrf (ref int M, ref int N, [In,Out] complex [] A, ref int LDA, int [] IPIV, ref int info); [DllImport("lapack_gen", EntryPoint = "dgetri"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dgetri (ref int N, double[] A,ref int LDA, int[] IPIV, double[] work, ref int lwork, ref int info); [DllImport("lapack_gen", EntryPoint = "sgetri"),SuppressUnmanagedCodeSecurity] private static extern void lapack_sgetri (ref int N, float [] A, ref int LDA, int [] IPIV, float[] work, ref int lwork, ref int info); [DllImport("lapack_gen", EntryPoint = "cgetri"),SuppressUnmanagedCodeSecurity] private static extern void lapack_cgetri (ref int N, [In,Out] fcomplex [] A, ref int LDA,[In,Out] int [] IPIV, [In,Out] fcomplex[] work, ref int lwork, ref int info); [DllImport("lapack_gen", EntryPoint = "zgetri"),SuppressUnmanagedCodeSecurity] private static extern void lapack_zgetri (ref int N, [In,Out] complex [] A, ref int LDA, int [] IPIV, [In,Out] complex[] work, ref int lwork, ref int info); [DllImport("lapack_gen", EntryPoint = "dgeqrf"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dgeqrf (ref int M, ref int N, double [] A, ref int lda, double [] tau, double[] work, ref int lwork, ref int info); [DllImport("lapack_gen", EntryPoint = "sgeqrf"),SuppressUnmanagedCodeSecurity] private static extern void lapack_sgeqrf (ref int M, ref int N, float [] A, ref int lda, float [] tau, float[] work, ref int lwork, ref int info); [DllImport("lapack_gen", EntryPoint = "cgeqrf"),SuppressUnmanagedCodeSecurity] private static extern void lapack_cgeqrf (ref int M, ref int N, [In,Out] fcomplex [] A, ref int lda, [In,Out] fcomplex [] tau, [In,Out] fcomplex[] work, ref int lwork, ref int info); [DllImport("lapack_gen", EntryPoint = "zgeqrf"),SuppressUnmanagedCodeSecurity] private static extern void lapack_zgeqrf (ref int M, ref int N, [In,Out] complex [] A, ref int lda, [In,Out] complex [] tau, [In,Out] complex[] work, ref int lwork, ref int info); [DllImport("lapack_gen", EntryPoint = "dgeqp3"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dgeqp3 (ref int M, ref int N, double [] A, ref int LDA, int [] JPVT, double [] tau, double [] work, ref int lwork, ref int info ); [DllImport("lapack_gen", EntryPoint = "sgeqp3"),SuppressUnmanagedCodeSecurity] private static extern void lapack_sgeqp3 (ref int M, ref int N, float [] A, ref int LDA, int [] JPVT, float [] tau, float [] work, ref int lwork, ref int info ); [DllImport("lapack_gen", EntryPoint = "cgeqp3"),SuppressUnmanagedCodeSecurity] private static extern void lapack_cgeqp3 (ref int M, ref int N,[In,Out] fcomplex [] A, ref int LDA,[In,Out] int [] JPVT, [In,Out] fcomplex [] tau, [In,Out] fcomplex [] work, ref int lwork,[In,Out] float [] rwork, ref int info ); [DllImport("lapack_gen", EntryPoint = "zgeqp3"),SuppressUnmanagedCodeSecurity] private static extern void lapack_zgeqp3 (ref int M, ref int N, [In,Out] complex [] A, ref int LDA, [In,Out] int [] JPVT, [In,Out] complex [] tau, [In,Out] complex [] work, ref int lwork, [In,Out] double [] rwork, ref int info ); [DllImport("lapack_gen", EntryPoint = "dormqr"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dormqr (ref char side, ref char trans, ref int m, ref int n, ref int k, double[] A, int lda, double[] tau, double[] C, ref int ldc, ref int info); [DllImport("lapack_gen", EntryPoint = "sormqr"),SuppressUnmanagedCodeSecurity] private static extern void lapack_sormqr (ref char side, ref char trans, ref int m, ref int n, ref int k, float [] A, ref int lda, float [] tau , float [] C, ref int ldc, ref int info); [DllImport("lapack_gen", EntryPoint = "dorgqr"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dorgqr (ref int m, ref int n, ref int k, double[] A,ref int lda, double[] tau, double[] work, ref int lwork, ref int info); [DllImport("lapack_gen", EntryPoint = "sorgqr"),SuppressUnmanagedCodeSecurity] private static extern void lapack_sorgqr (ref int m, ref int n, ref int k, float [] A, ref int lda, float [] tau , float[] work, ref int lwork, ref int info); [DllImport("lapack_gen", EntryPoint = "cungqr"),SuppressUnmanagedCodeSecurity] private static extern void lapack_cungqr (ref int m, ref int n, ref int k, [In,Out] fcomplex[] A,ref int lda, [In,Out] fcomplex[] tau, [In,Out] fcomplex[] work, ref int lwork, ref int info); [DllImport("lapack_gen", EntryPoint = "zungqr"),SuppressUnmanagedCodeSecurity] private static extern void lapack_zungqr (ref int m, ref int n, ref int k, [In,Out] complex[] A,ref int lda, [In,Out] complex[] tau, [In,Out] complex[] work, ref int lwork, ref int info); [DllImport("lapack_gen", EntryPoint = "dtrtrs"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dtrtrs (ref char uplo, ref char transA, ref char diag, ref int N, ref int nrhs, IntPtr A, ref int LDA, IntPtr B, ref int LDB, ref int info); [DllImport("lapack_gen", EntryPoint = "strtrs"),SuppressUnmanagedCodeSecurity] private static extern void lapack_strtrs (ref char uplo, ref char transA, ref char diag, ref int N, ref int nrhs, IntPtr A, ref int LDA, IntPtr B, ref int LDB, ref int info); [DllImport("lapack_gen", EntryPoint = "ctrtrs"),SuppressUnmanagedCodeSecurity] private static extern void lapack_ctrtrs (ref char uplo, ref char transA, ref char diag, ref int N, ref int nrhs, IntPtr A, ref int LDA, IntPtr B, ref int LDB, ref int info); [DllImport("lapack_gen", EntryPoint = "ztrtrs"),SuppressUnmanagedCodeSecurity] private static extern void lapack_ztrtrs (ref char uplo, ref char transA, ref char diag, ref int N, ref int nrhs, IntPtr A, ref int LDA, IntPtr B, ref int LDB, ref int info); [DllImport("lapack_gen", EntryPoint = "dgetrs"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dgetrs(ref char trans, ref int N, ref int NRHS, double[] A, ref int LDA, int[] IPIV, double[] B, ref int LDB, ref int info); [DllImport("lapack_gen", EntryPoint = "sgetrs"),SuppressUnmanagedCodeSecurity] private static extern void lapack_sgetrs(ref char trans, ref int N, ref int NRHS, float[] A, ref int LDA, int[] IPIV, float[] B, ref int LDB, ref int info); [DllImport("lapack_gen", EntryPoint = "cgetrs"),SuppressUnmanagedCodeSecurity] private static extern void lapack_cgetrs(ref char trans, ref int N, ref int NRHS, [In,Out] fcomplex[] A, ref int LDA, int[] IPIV, [In,Out] fcomplex[] B, ref int LDB, ref int info); [DllImport("lapack_gen", EntryPoint = "zgetrs"),SuppressUnmanagedCodeSecurity] private static extern void lapack_zgetrs(ref char trans, ref int N, ref int NRHS, [In,Out] complex[] A, ref int LDA, int[] IPIV, [In,Out] complex[] B, ref int LDB, ref int info); [DllImport("lapack_gen", EntryPoint = "dgelsd"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dgelsd (ref int m,ref int n,ref int nrhs, double[] A,ref int lda, double[] B,ref int ldb, double[] S, ref double RCond, ref int rank, double[] work,ref int lwork, int[] iwork, ref int info); [DllImport("lapack_gen", EntryPoint = "sgelsd"),SuppressUnmanagedCodeSecurity] private static extern void lapack_sgelsd (ref int m,ref int n,ref int nrhs, float[] A,ref int lda, float[] B,ref int ldb, float[] S, ref float RCond, ref int rank, float[] work,ref int lwork, int[] iwork, ref int info); [DllImport("lapack_gen", EntryPoint = "cgelsd"),SuppressUnmanagedCodeSecurity] private static extern void lapack_cgelsd (ref int m,ref int n,ref int nrhs, [In,Out] fcomplex[] A,ref int lda, [In,Out] fcomplex[] B,ref int ldb, float[] S , ref float RCond, ref int rank, [In,Out] fcomplex[] work, ref int lwork, float[] rwork, int[] iwork, ref int info); [DllImport("lapack_gen", EntryPoint = "zgelsd"),SuppressUnmanagedCodeSecurity] private static extern void lapack_zgelsd (ref int m,ref int n,ref int nrhs, [In,Out] complex[] A,ref int lda, [In,Out] complex[] B,ref int ldb, double[] S, ref double RCond, ref int rank, [In,Out] complex[] work,ref int lwork, double[] rwork, int[] iwork, ref int info); [DllImport("lapack_gen", EntryPoint = "dgelsy"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dgelsy (ref int m,ref int n,ref int nrhs, double[] A,ref int lda, double[] B,ref int ldb, int[] JPVT0, ref double RCond, ref int rank, double[] work, ref int lwork, ref int info); [DllImport("lapack_gen", EntryPoint = "sgelsy"),SuppressUnmanagedCodeSecurity] private static extern void lapack_sgelsy (ref int m,ref int n,ref int nrhs, float[] A,ref int lda, float[] B,ref int ldb, int[] JPVT0, ref float RCond, ref int rank, float[] work,ref int lwork, ref int info); [DllImport("lapack_gen", EntryPoint = "cgelsy"),SuppressUnmanagedCodeSecurity] private static extern void lapack_cgelsy (ref int m,ref int n,ref int nrhs, [In,Out] fcomplex[] A,ref int lda, [In,Out] fcomplex[] B,ref int ldb, int[] JPVT0, ref float RCond, ref int rank, [In,Out] fcomplex[] work, ref int lwork, float[] rwork, ref int info); [DllImport("lapack_gen", EntryPoint = "zgelsy"),SuppressUnmanagedCodeSecurity] private static extern void lapack_zgelsy (ref int m,ref int n,ref int nrhs, [In,Out] complex[] A,ref int lda, [In,Out] complex[] B,ref int ldb, int[] JPVT0, ref double RCond, ref int rank, [In,Out] complex[] work, ref int lwork, double[] rwork, ref int info); [DllImport("lapack_gen", EntryPoint = "dgeevx"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dgeevx (ref char balance, ref char jobvl, ref char jobvr, ref char sense, ref int n, double[] A, ref int lda, double[] wr, double[] wi, double[] vl, ref int ldvl, double[] vr, ref int ldvr, ref int ilo, ref int ihi, double[] scale, ref double abnrm, double[] rconde, double[] rcondv, double[] work, ref int lwork, int [] iwork, ref int info); [DllImport("lapack_gen", EntryPoint = "sgeevx"),SuppressUnmanagedCodeSecurity] private static extern void lapack_sgeevx (ref char balance, ref char jobvl, ref char jobvr, ref char sense, ref int n, float[] A, ref int lda, float[] wr, float[] wi, float[] vl, ref int ldvl, float[] vr, ref int ldvr, ref int ilo, ref int ihi, float[] scale, ref float abnrm, float[] rconde, float[] rcondv, float [] work, ref int lwork, int [] iwork, ref int info); [DllImport("lapack_gen", EntryPoint = "cgeevx"),SuppressUnmanagedCodeSecurity] private static extern void lapack_cgeevx (ref char balance, ref char jobvl, ref char jobvr, ref char sense, ref int n, [In,Out] fcomplex[] A, ref int lda, [In,Out] fcomplex[] w, [In,Out] fcomplex[] vl, ref int ldvl, [In,Out] fcomplex[] vr, ref int ldvr, ref int ilo, ref int ihi, float[] scale, ref float abnrm, float[] rconde, float[] rcondv, [In,Out] fcomplex[] work, ref int lwork, float[] rwork, ref int info); [DllImport("lapack_gen", EntryPoint = "zgeevx"),SuppressUnmanagedCodeSecurity] private static extern void lapack_zgeevx (ref char balance, ref char jobvl, ref char jobvr, ref char sense, ref int n, [In,Out] complex[] A, ref int lda, [In,Out] complex[] w, [In,Out] complex[] vl, ref int ldvl, [In,Out] complex[] vr, ref int ldvr, ref int ilo, ref int ihi, double[] scale, ref double abnrm, double[] rconde, double[] rcondv, [In,Out] complex[] work, ref int lwork, double[] rwork, ref int info); [DllImport("lapack_gen", EntryPoint = "dsyevr"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dsyevr (ref char jobz, ref char range, ref char uplo, ref int n, double [] A, ref int lda, ref double vl, ref double vu, ref int il, ref int iu, ref double abstol, ref int m, double[] w, double [] z, ref int ldz, int[] isuppz, double[] work, ref int lwork, int[] iwork, ref int liwork, ref int info); [DllImport("lapack_gen", EntryPoint = "ssyevr"),SuppressUnmanagedCodeSecurity] private static extern void lapack_ssyevr (ref char jobz, ref char range, ref char uplo, ref int n, float [] A, ref int lda, ref float vl, ref float vu, ref int il, ref int iu, ref float abstol, ref int m, float [] w, float [] z, ref int ldz, int[] isuppz, float [] work, ref int lwork, int[] iwork, ref int liwork, ref int info); [DllImport("lapack_gen", EntryPoint = "cheevr"),SuppressUnmanagedCodeSecurity] private static extern void lapack_cheevr (ref char jobz, ref char range, ref char uplo, ref int n, [In,Out] fcomplex[] A, ref int lda, ref float vl, ref float vu, ref int il, ref int iu, ref float abstol, ref int m, float [] w, [In,Out] fcomplex[] z, ref int ldz, int[] isuppz, [In,Out] fcomplex[] work, ref int lwork, float[] rwork, ref int lrwork, int[] iwork, ref int liwork, ref int info); [DllImport("lapack_gen", EntryPoint = "zheevr"),SuppressUnmanagedCodeSecurity] private static extern void lapack_zheevr (ref char jobz, ref char range, ref char uplo, ref int n, [In,Out] complex [] A, ref int lda, ref double vl, ref double vu, ref int il, ref int iu, ref double abstol, ref int m, double[] w, [In,Out] complex [] z, ref int ldz, int[] isuppz, [In,Out] complex[] work, ref int lwork, double[] rwork, ref int lrwork, int[] iwork, ref int liwork, ref int info); [DllImport("lapack_gen", EntryPoint = "dsygv"),SuppressUnmanagedCodeSecurity] private static extern void lapack_dsygv (ref int itype, ref char jobz, ref char uplo, ref int n, double [] A, ref int lda, double [] B, ref int ldb, double [] w, double [] work, ref int lwork, ref int info); [DllImport("lapack_gen", EntryPoint = "ssygv"),SuppressUnmanagedCodeSecurity] private static extern void lapack_ssygv (ref int itype, ref char jobz, ref char uplo, ref int n, float [] A, ref int lda, float [] B, ref int ldb, float [] w, float [] work, ref int lwork, ref int info); [DllImport("lapack_gen", EntryPoint = "chegv"),SuppressUnmanagedCodeSecurity] private static extern void lapack_chegv (ref int itype, ref char jobz, ref char uplo, ref int n, [In,Out] fcomplex[] A, ref int lda, [In,Out] fcomplex [] B, ref int ldb, float [] w, [In,Out] fcomplex[] work , ref int lwork, float[] rwork, ref int info); [DllImport("lapack_gen", EntryPoint = "zhegv"),SuppressUnmanagedCodeSecurity] private static extern void lapack_zhegv (ref int itype, ref char jobz, ref char uplo, ref int n, [In,Out] complex [] A, ref int lda, [In,Out] complex [] B, ref int ldb, double [] w, [In,Out] complex[] work , ref int lwork, double[] rwork, ref int info); #endregion DLL INCLUDES #region ILLAPACK INTERFACE private int ILAENV (int ispec, string name, string opts, int n1, int n2, int n3, int n4) { return lapack_ilaenv(ref ispec, ref name, ref opts, ref n1, ref n2, ref n3, ref n4); } ///

/// Implement wrapper for ATLAS GeneralMatrixMultiply ///

/// Transposition state for matrix A: one of the constants in enum CBlas_Transpose /// Transposition state for matrix B: one of the constants in enum CBlas_Transpose /// Number of rows in A /// Number of columns in B /// Number of columns in A and number of rows in B /// multiplicationi factor for A /// pointer to double array A /// distance between first elements of each column for column based orientation or /// distance between first elements of each row for row based orientation for matrix B /// pointer to double array B /// distance between first elements of each column for column based orientation or /// distance between first elements of each row for row based orientation for matrix A /// multiplication faktor for matrix B /// pointer to predefined double array C of neccessary length /// distance between first elements of each column for column based orientation or /// distance between first elements of each row for row based orientation for matrix C /// All parameters except C are readonly. Only elements of matrix C will be altered. C must be a predefined /// continous double array of size MxN public void zgemm(char TransA, char TransB, int M, int N, int K, complex alpha, IntPtr A, int lda, IntPtr B, int ldb, complex beta, complex[] C, int ldc) { lapack_zgemm(ref TransA, ref TransB, ref M, ref N, ref K, ref alpha, A, ref lda, B, ref ldb, ref beta, C, ref ldc); } public void dgesdd (char jobz, int m, int n, double [] a, int lda, double [] s, double [] u, int ldu, double [] vt, int ldvt, ref int info) { try { double [] work = new double [1] { ( double )0.0 }; int lwork = -1; int[] iwork = new int[((m < n) ? m : n) * 8]; lapack_dgesdd (ref jobz, ref m, ref n, a, ref lda, s, u, ref ldu, vt, ref ldvt, work, ref lwork, iwork, ref info); if (work[0] != 0) { work = new double [(int)work[0]]; lwork = work.Length; lapack_dgesdd (ref jobz, ref m, ref n, a, ref lda, s, u, ref ldu, vt, ref ldvt, work, ref lwork, iwork, ref info); } } catch (Exception e) { if (e is OutOfMemoryException) { dgesvd (jobz, m, n, a, lda, s, u, ldu, vt, ldvt, ref info); } throw new ILException("Unable to do " + "dgesdd" + ".", e); } } #region HYCALPER AUTO GENERATED CODE public void sgesdd (char jobz, int m, int n, float [] a, int lda, float [] s, float [] u, int ldu, float [] vt, int ldvt, ref int info) { try { float [] work = new float [1] { ( float )0.0 }; int lwork = -1; int[] iwork = new int[((m < n) ? m : n) * 8]; lapack_sgesdd (ref jobz, ref m, ref n, a, ref lda, s, u, ref ldu, vt, ref ldvt, work, ref lwork, iwork, ref info); if (work[0] != 0) { work = new float [(int)work[0]]; lwork = work.Length; lapack_sgesdd (ref jobz, ref m, ref n, a, ref lda, s, u, ref ldu, vt, ref ldvt, work, ref lwork, iwork, ref info); } } catch (Exception e) { if (e is OutOfMemoryException) { sgesvd (jobz, m, n, a, lda, s, u, ldu, vt, ldvt, ref info); } throw new ILException("Unable to do " + "sgesdd" + ".", e); } } #endregion HYCALPER AUTO GENERATED CODE public void zgesdd (char jobz, int m, int n, complex [] a, int lda, double [] s, complex [] u, int ldu, complex [] vt, int ldvt, ref int info) { try { complex [] work = new complex [1] { ( complex )0.0 }; double [] rwork; int minMN = (m < n) ? m : n; if (jobz == 'N') { rwork = new double [minMN * 7]; } else { rwork = new double [5 * minMN * minMN + 5 * minMN]; } int lwork = -1; int[] iwork = new int[minMN * 8]; lapack_zgesdd (ref jobz, ref m, ref n, a, ref lda, s, u, ref ldu, vt, ref ldvt, work, ref lwork, rwork, iwork, ref info); if (work[0] != 0) { work = new complex [(int)work[0].real]; lwork = work.Length; lapack_zgesdd (ref jobz, ref m, ref n, a, ref lda, s, u, ref ldu, vt, ref ldvt, work, ref lwork,rwork, iwork, ref info); } } catch (Exception e) { if (e is OutOfMemoryException) { zgesvd (jobz, m, n, a, lda, s, u, ldu, vt, ldvt, ref info); } throw new ILException("Unable to do " + "zgesdd" + ".", e); } } #region HYCALPER AUTO GENERATED CODE public void cgesdd (char jobz, int m, int n, fcomplex [] a, int lda, float [] s, fcomplex [] u, int ldu, fcomplex [] vt, int ldvt, ref int info) { try { fcomplex [] work = new fcomplex [1] { ( fcomplex )0.0 }; float [] rwork; int minMN = (m < n) ? m : n; if (jobz == 'N') { rwork = new float [minMN * 7]; } else { rwork = new float [5 * minMN * minMN + 5 * minMN]; } int lwork = -1; int[] iwork = new int[minMN * 8]; lapack_cgesdd (ref jobz, ref m, ref n, a, ref lda, s, u, ref ldu, vt, ref ldvt, work, ref lwork, rwork, iwork, ref info); if (work[0] != 0) { work = new fcomplex [(int)work[0].real]; lwork = work.Length; lapack_cgesdd (ref jobz, ref m, ref n, a, ref lda, s, u, ref ldu, vt, ref ldvt, work, ref lwork,rwork, iwork, ref info); } } catch (Exception e) { if (e is OutOfMemoryException) { cgesvd (jobz, m, n, a, lda, s, u, ldu, vt, ldvt, ref info); } throw new ILException("Unable to do " + "zgesdd" + ".", e); } } #endregion HYCALPER AUTO GENERATED CODE ///

/// singular value decomposition ///

/// /// /// /// /// /// /// /// /// /// /// public void dgesvd (char jobz, int m, int n, double [] a, int lda, double [] s, double [] u, int ldu, double [] vt, int ldvt, ref int info) { if (jobz != 'A' && jobz != 'S' && jobz != 'N') throw new ILArgumentException("Argument jobz must be one of 'A','S' or 'N'"); try { double [] work = new double [1] { ( double )0.0 }; int lwork = -1; int[] iwork = new int[((m < n) ? m : n) * 8]; lapack_dgesvd (ref jobz, ref jobz, ref m, ref n, a, ref lda, s, u, ref ldu, vt, ref ldvt, work, ref lwork, iwork, ref info); if (work[0] != 0) { work = new double [(int)work[0]]; lwork = work.Length; lapack_dgesvd (ref jobz, ref jobz, ref m, ref n, a, ref lda, s, u, ref ldu, vt, ref ldvt, work, ref lwork, iwork, ref info); } } catch (Exception e) { if (e is OutOfMemoryException) { throw new ILMemoryException("Not enough memory for given arguments."); } throw new ILException("Unable to do gesvd.", e); } } #region HYCALPER AUTO GENERATED CODE ///

/// singular value decomposition ///

/// /// /// /// /// /// /// /// /// /// /// public void sgesvd (char jobz, int m, int n, float [] a, int lda, float [] s, float [] u, int ldu, float [] vt, int ldvt, ref int info) { if (jobz != 'A' && jobz != 'S' && jobz != 'N') throw new ILArgumentException("Argument jobz must be one of 'A','S' or 'N'"); try { float [] work = new float [1] { ( float )0.0 }; int lwork = -1; int[] iwork = new int[((m < n) ? m : n) * 8]; lapack_sgesvd (ref jobz, ref jobz, ref m, ref n, a, ref lda, s, u, ref ldu, vt, ref ldvt, work, ref lwork, iwork, ref info); if (work[0] != 0) { work = new float [(int)work[0]]; lwork = work.Length; lapack_sgesvd (ref jobz, ref jobz, ref m, ref n, a, ref lda, s, u, ref ldu, vt, ref ldvt, work, ref lwork, iwork, ref info); } } catch (Exception e) { if (e is OutOfMemoryException) { throw new ILMemoryException("Not enough memory for given arguments."); } throw new ILException("Unable to do gesvd.", e); } } ///

/// singular value decomposition ///

/// /// /// /// /// /// /// /// /// /// /// public void zgesvd (char jobz, int m, int n, complex [] a, int lda, double [] s, complex [] u, int ldu, complex [] vt, int ldvt, ref int info) { if (jobz != 'A' && jobz != 'S' && jobz != 'N') throw new ILArgumentException("Argument jobz must be one of 'A','S' or 'N'"); try { complex [] work = new complex [1] { ( complex )0.0 }; int lwork = -1; int[] iwork = new int[((m < n) ? m : n) * 8]; lapack_zgesvd (ref jobz, ref jobz, ref m, ref n, a, ref lda, s, u, ref ldu, vt, ref ldvt, work, ref lwork, iwork, ref info); if (work[0] != 0) { work = new complex [(int)work[0]]; lwork = work.Length; lapack_zgesvd (ref jobz, ref jobz, ref m, ref n, a, ref lda, s, u, ref ldu, vt, ref ldvt, work, ref lwork, iwork, ref info); } } catch (Exception e) { if (e is OutOfMemoryException) { throw new ILMemoryException("Not enough memory for given arguments."); } throw new ILException("Unable to do gesvd.", e); } } ///

/// singular value decomposition ///

/// /// /// /// /// /// /// /// /// /// /// public void cgesvd (char jobz, int m, int n, fcomplex [] a, int lda, float [] s, fcomplex [] u, int ldu, fcomplex [] vt, int ldvt, ref int info) { if (jobz != 'A' && jobz != 'S' && jobz != 'N') throw new ILArgumentException("Argument jobz must be one of 'A','S' or 'N'"); try { fcomplex [] work = new fcomplex [1] { ( fcomplex )0.0 }; int lwork = -1; int[] iwork = new int[((m < n) ? m : n) * 8]; lapack_cgesvd (ref jobz, ref jobz, ref m, ref n, a, ref lda, s, u, ref ldu, vt, ref ldvt, work, ref lwork, iwork, ref info); if (work[0] != 0) { work = new fcomplex [(int)work[0]]; lwork = work.Length; lapack_cgesvd (ref jobz, ref jobz, ref m, ref n, a, ref lda, s, u, ref ldu, vt, ref ldvt, work, ref lwork, iwork, ref info); } } catch (Exception e) { if (e is OutOfMemoryException) { throw new ILMemoryException("Not enough memory for given arguments."); } throw new ILException("Unable to do gesvd.", e); } } #endregion HYCALPER AUTO GENERATED CODE public void dpotrf(char uplo, int n, double[] A, int lda, ref int info) { lapack_dpotrf(ref uplo, ref n, A, ref lda, ref info); } public void spotrf(char uplo, int n, float[] A, int lda, ref int info) { lapack_spotrf(ref uplo, ref n, A, ref lda, ref info); } public void cpotrf(char uplo, int n, fcomplex[] A, int lda, ref int info) { lapack_cpotrf(ref uplo, ref n, A, ref lda, ref info); } public void zpotrf(char uplo, int n, complex[] A, int lda, ref int info) { lapack_zpotrf(ref uplo, ref n, A, ref lda, ref info); } public void dpotri(char uplo, int n, double[] A, int lda, ref int info) { lapack_dpotri(ref uplo, ref n, A, ref lda, ref info); } public void spotri(char uplo, int n, float[] A, int lda, ref int info) { lapack_spotri(ref uplo, ref n, A, ref lda, ref info); } public void cpotri(char uplo, int n, fcomplex[] A, int lda, ref int info) { lapack_cpotri(ref uplo, ref n, A, ref lda, ref info); } public void zpotri(char uplo, int n, complex[] A, int lda, ref int info) { lapack_zpotri(ref uplo, ref n, A, ref lda, ref info); } public void dgetrf(int M, int N, double[] A, int LDA, int[] IPIV, ref int info) { lapack_dgetrf(ref M, ref N, A, ref LDA, IPIV, ref info); } public void sgetrf(int M, int N, float[] A, int LDA, int[] IPIV, ref int info) { lapack_sgetrf(ref M, ref N, A, ref LDA, IPIV, ref info); } public void cgetrf(int M, int N, fcomplex[] A, int LDA, int[] IPIV, ref int info) { lapack_cgetrf(ref M, ref N, A, ref LDA, IPIV, ref info); } public void zgetrf(int M, int N, complex[] A, int LDA, int[] IPIV, ref int info) { lapack_zgetrf(ref M, ref N, A, ref LDA, IPIV, ref info); } public void dgetri (int N, double [] A, int LDA, int[] IPIV, ref int info) { double [] work = new double [1]; int lwork = -1; try { lapack_dgetri (ref N, A, ref LDA, IPIV, work, ref lwork, ref info); lwork = (int)work[0]; if (lwork > 0 && info == 0) { work = new double [lwork]; lapack_dgetri (ref N, A, ref LDA, IPIV, work, ref lwork, ref info); } else { throw new ILException("error in lapack_dgetri"); } } catch (OutOfMemoryException e) { throw new ILException("error on dgetri. Not enough memory! " + (lwork * Marshal.SizeOf( work[0] )).ToString() + " bytes has been requested.",e); } } #region HYCALPER AUTO GENERATED CODE public void sgetri (int N, float [] A, int LDA, int[] IPIV, ref int info) { float [] work = new float [1]; int lwork = -1; try { lapack_sgetri (ref N, A, ref LDA, IPIV, work, ref lwork, ref info); lwork = (int)work[0]; if (lwork > 0 && info == 0) { work = new float [lwork]; lapack_sgetri (ref N, A, ref LDA, IPIV, work, ref lwork, ref info); } else { throw new ILException("error in lapack_dgetri"); } } catch (OutOfMemoryException e) { throw new ILException("error on dgetri. Not enough memory! " + (lwork * Marshal.SizeOf( work[0] )).ToString() + " bytes has been requested.",e); } } public void zgetri (int N, complex [] A, int LDA, int[] IPIV, ref int info) { complex [] work = new complex [1]; int lwork = -1; try { lapack_zgetri (ref N, A, ref LDA, IPIV, work, ref lwork, ref info); lwork = (int)work[0]; if (lwork > 0 && info == 0) { work = new complex [lwork]; lapack_zgetri (ref N, A, ref LDA, IPIV, work, ref lwork, ref info); } else { throw new ILException("error in lapack_dgetri"); } } catch (OutOfMemoryException e) { throw new ILException("error on dgetri. Not enough memory! " + (lwork * Marshal.SizeOf( work[0] )).ToString() + " bytes has been requested.",e); } } public void cgetri (int N, fcomplex [] A, int LDA, int[] IPIV, ref int info) { fcomplex [] work = new fcomplex [1]; int lwork = -1; try { lapack_cgetri (ref N, A, ref LDA, IPIV, work, ref lwork, ref info); lwork = (int)work[0]; if (lwork > 0 && info == 0) { work = new fcomplex [lwork]; lapack_cgetri (ref N, A, ref LDA, IPIV, work, ref lwork, ref info); } else { throw new ILException("error in lapack_dgetri"); } } catch (OutOfMemoryException e) { throw new ILException("error on dgetri. Not enough memory! " + (lwork * Marshal.SizeOf( work[0] )).ToString() + " bytes has been requested.",e); } } #endregion HYCALPER AUTO GENERATED CODE public void dgeqrf (int M, int N, double [] A, int lda, double [] tau, ref int info) { double [] work = new double [1]; int lwork = -1; try { lapack_dgeqrf (ref M, ref N, A, ref lda, tau, work, ref lwork, ref info); lwork = (int)work[0]; if (lwork > 0 && info == 0) { work = new double [lwork]; lapack_dgeqrf (ref M, ref N, A, ref lda, tau, work, ref lwork, ref info); } else { throw new ILException("error in lapack_?geqrf"); } } catch (OutOfMemoryException e) { throw new ILException("error on ?geqrf. Not enough memory! " + (lwork * Marshal.SizeOf( work[0] )).ToString() + " bytes has been requested.",e); } } #region HYCALPER AUTO GENERATED CODE public void sgeqrf (int M, int N, float [] A, int lda, float [] tau, ref int info) { float [] work = new float [1]; int lwork = -1; try { lapack_sgeqrf (ref M, ref N, A, ref lda, tau, work, ref lwork, ref info); lwork = (int)work[0]; if (lwork > 0 && info == 0) { work = new float [lwork]; lapack_sgeqrf (ref M, ref N, A, ref lda, tau, work, ref lwork, ref info); } else { throw new ILException("error in lapack_?geqrf"); } } catch (OutOfMemoryException e) { throw new ILException("error on ?geqrf. Not enough memory! " + (lwork * Marshal.SizeOf( work[0] )).ToString() + " bytes has been requested.",e); } } public void zgeqrf (int M, int N, complex [] A, int lda, complex [] tau, ref int info) { complex [] work = new complex [1]; int lwork = -1; try { lapack_zgeqrf (ref M, ref N, A, ref lda, tau, work, ref lwork, ref info); lwork = (int)work[0]; if (lwork > 0 && info == 0) { work = new complex [lwork]; lapack_zgeqrf (ref M, ref N, A, ref lda, tau, work, ref lwork, ref info); } else { throw new ILException("error in lapack_?geqrf"); } } catch (OutOfMemoryException e) { throw new ILException("error on ?geqrf. Not enough memory! " + (lwork * Marshal.SizeOf( work[0] )).ToString() + " bytes has been requested.",e); } } public void cgeqrf (int M, int N, fcomplex [] A, int lda, fcomplex [] tau, ref int info) { fcomplex [] work = new fcomplex [1]; int lwork = -1; try { lapack_cgeqrf (ref M, ref N, A, ref lda, tau, work, ref lwork, ref info); lwork = (int)work[0]; if (lwork > 0 && info == 0) { work = new fcomplex [lwork]; lapack_cgeqrf (ref M, ref N, A, ref lda, tau, work, ref lwork, ref info); } else { throw new ILException("error in lapack_?geqrf"); } } catch (OutOfMemoryException e) { throw new ILException("error on ?geqrf. Not enough memory! " + (lwork * Marshal.SizeOf( work[0] )).ToString() + " bytes has been requested.",e); } } #endregion HYCALPER AUTO GENERATED CODE public void dormqr(char side, char trans, int m, int n, int k, double[] A, int lda, double[] tau, double[] C, int LDC, ref int info) { throw new Exception("The method or operation is not implemented."); } public void sormqr(char side, char trans, int m, int n, int k, float[] A, int lda, float[] tau, float[] C, int LDC, ref int info) { throw new Exception("The method or operation is not implemented."); } public void dorgqr (int M, int N, int K, double [] A, int lda, double [] tau, ref int info) { double [] work = new double [1]; int lwork = -1; try { /*!HC:lapack_***gqr*/ lapack_dorgqr (ref M, ref N, ref K, A, ref lda, tau, work, ref lwork, ref info); lwork = (int)work[0]; if (lwork > 0 && info == 0) { work = new double [lwork]; /*!HC:lapack_***gqr*/ lapack_dorgqr (ref M, ref N, ref K, A, ref lda, tau, work, ref lwork, ref info); } else { throw new ILException("error in lapack_?[un/or]gqr"); } } catch (OutOfMemoryException e) { throw new ILException("error on ?[un/or]gqr. Not enough memory! " + (lwork * Marshal.SizeOf( work[0] )).ToString() + " bytes has been requested.",e); } } #region HYCALPER AUTO GENERATED CODE public void sorgqr (int M, int N, int K, float [] A, int lda, float [] tau, ref int info) { float [] work = new float [1]; int lwork = -1; try { lapack_sorgqr (ref M, ref N, ref K, A, ref lda, tau, work, ref lwork, ref info); lwork = (int)work[0]; if (lwork > 0 && info == 0) { work = new float [lwork]; lapack_sorgqr (ref M, ref N, ref K, A, ref lda, tau, work, ref lwork, ref info); } else { throw new ILException("error in lapack_?[un/or]gqr"); } } catch (OutOfMemoryException e) { throw new ILException("error on ?[un/or]gqr. Not enough memory! " + (lwork * Marshal.SizeOf( work[0] )).ToString() + " bytes has been requested.",e); } } public void zungqr (int M, int N, int K, complex [] A, int lda, complex [] tau, ref int info) { complex [] work = new complex [1]; int lwork = -1; try { lapack_zungqr (ref M, ref N, ref K, A, ref lda, tau, work, ref lwork, ref info); lwork = (int)work[0]; if (lwork > 0 && info == 0) { work = new complex [lwork]; lapack_zungqr (ref M, ref N, ref K, A, ref lda, tau, work, ref lwork, ref info); } else { throw new ILException("error in lapack_?[un/or]gqr"); } } catch (OutOfMemoryException e) { throw new ILException("error on ?[un/or]gqr. Not enough memory! " + (lwork * Marshal.SizeOf( work[0] )).ToString() + " bytes has been requested.",e); } } public void cungqr (int M, int N, int K, fcomplex [] A, int lda, fcomplex [] tau, ref int info) { fcomplex [] work = new fcomplex [1]; int lwork = -1; try { lapack_cungqr (ref M, ref N, ref K, A, ref lda, tau, work, ref lwork, ref info); lwork = (int)work[0]; if (lwork > 0 && info == 0) { work = new fcomplex [lwork]; lapack_cungqr (ref M, ref N, ref K, A, ref lda, tau, work, ref lwork, ref info); } else { throw new ILException("error in lapack_?[un/or]gqr"); } } catch (OutOfMemoryException e) { throw new ILException("error on ?[un/or]gqr. Not enough memory! " + (lwork * Marshal.SizeOf( work[0] )).ToString() + " bytes has been requested.",e); } } #endregion HYCALPER AUTO GENERATED CODE public void dgeqp3 ( int M,int N, double [] A,int LDA,int [] JPVT, double [] tau, ref int info ) { double [] work = new double [1]; int lwork = -1; try { /*dummy*/ lapack_dgeqp3 (ref M, ref N, A, ref LDA, JPVT, tau, work, ref lwork, ref info); lwork = (int)work[0]; if (lwork > 0 && info == 0) { work = new double [lwork]; lapack_dgeqp3 (ref M, ref N, A, ref LDA, JPVT, tau, work, ref lwork, ref info); } else { throw new ILException("error in lapack_?geqp3"); } } catch (OutOfMemoryException e) { throw new ILException("error on ?geqp3. Not enough memory! " + (lwork * Marshal.SizeOf( work[0] )).ToString() + " bytes has been requested.",e); } } #region HYCALPER AUTO GENERATED CODE public void sgeqp3 ( int M,int N, float [] A,int LDA,int [] JPVT, float [] tau, ref int info ) { float [] work = new float [1]; int lwork = -1; try { lapack_sgeqp3 (ref M, ref N, A, ref LDA, JPVT, tau, work, ref lwork, ref info); lwork = (int)work[0]; if (lwork > 0 && info == 0) { work = new float [lwork]; lapack_sgeqp3 (ref M, ref N, A, ref LDA, JPVT, tau, work, ref lwork, ref info); } else { throw new ILException("error in lapack_?geqp3"); } } catch (OutOfMemoryException e) { throw new ILException("error on ?geqp3. Not enough memory! " + (lwork * Marshal.SizeOf( work[0] )).ToString() + " bytes has been requested.",e); } } public void zgeqp3 ( int M,int N, complex [] A,int LDA,int [] JPVT, complex [] tau, ref int info ) { complex [] work = new complex [1]; int lwork = -1; try { double[] rwork = new double[2 * N]; lapack_zgeqp3 (ref M, ref N, A, ref LDA, JPVT, tau, work, ref lwork, rwork, ref info); lwork = (int)work[0]; if (lwork > 0 && info == 0) { work = new complex [lwork]; lapack_zgeqp3 (ref M, ref N, A, ref LDA, JPVT, tau, work, ref lwork, rwork, ref info); } else { throw new ILException("error in lapack_?geqp3"); } } catch (OutOfMemoryException e) { throw new ILException("error on ?geqp3. Not enough memory! " + (lwork * Marshal.SizeOf( work[0] )).ToString() + " bytes has been requested.",e); } } public void cgeqp3 ( int M,int N, fcomplex [] A,int LDA,int [] JPVT, fcomplex [] tau, ref int info ) { fcomplex [] work = new fcomplex [1]; int lwork = -1; try { float[] rwork = new float[2 * N]; lapack_cgeqp3 (ref M, ref N, A, ref LDA, JPVT, tau, work, ref lwork, rwork, ref info); lwork = (int)work[0]; if (lwork > 0 && info == 0) { work = new fcomplex [lwork]; lapack_cgeqp3 (ref M, ref N, A, ref LDA, JPVT, tau, work, ref lwork, rwork, ref info); } else { throw new ILException("error in lapack_?geqp3"); } } catch (OutOfMemoryException e) { throw new ILException("error on ?geqp3. Not enough memory! " + (lwork * Marshal.SizeOf( work[0] )).ToString() + " bytes has been requested.",e); } } #endregion HYCALPER AUTO GENERATED CODE public void dtrtrs(char uplo, char transA, char diag, int N, int nrhs, IntPtr A, int LDA, IntPtr B, int LDB, ref int info) { lapack_dtrtrs(ref uplo, ref transA, ref diag, ref N, ref nrhs, A, ref LDA, B, ref LDB, ref info); } public void strtrs(char uplo, char transA, char diag, int N, int nrhs, IntPtr A, int LDA, IntPtr B, int LDB, ref int info) { lapack_strtrs(ref uplo, ref transA, ref diag, ref N, ref nrhs, A, ref LDA, B, ref LDB, ref info); } public void ctrtrs(char uplo, char transA, char diag, int N, int nrhs, IntPtr A, int LDA, IntPtr B, int LDB, ref int info) { lapack_ctrtrs(ref uplo, ref transA, ref diag, ref N, ref nrhs, A, ref LDA, B, ref LDB, ref info); } public void ztrtrs(char uplo, char transA, char diag, int N, int nrhs, IntPtr A, int LDA, IntPtr B, int LDB, ref int info) { lapack_ztrtrs(ref uplo, ref transA, ref diag, ref N, ref nrhs, A, ref LDA, B, ref LDB, ref info); } public void dgetrs(char trans, int N, int NRHS, double[] A, int LDA, int[] IPIV, double[] B, int LDB, ref int info) { lapack_dgetrs(ref trans,ref N,ref NRHS, A, ref LDA, IPIV, B, ref LDB, ref info); } public void sgetrs(char trans, int N, int NRHS, float[] A, int LDA, int[] IPIV, float[] B, int LDB, ref int info) { lapack_sgetrs(ref trans,ref N,ref NRHS, A, ref LDA, IPIV, B, ref LDB, ref info); } public void cgetrs(char trans, int N, int NRHS, fcomplex[] A, int LDA, int[] IPIV, fcomplex[] B, int LDB, ref int info) { lapack_cgetrs(ref trans,ref N,ref NRHS, A, ref LDA, IPIV, B, ref LDB, ref info); } public void zgetrs(char trans, int N, int NRHS, complex[] A, int LDA, int[] IPIV, complex[] B, int LDB, ref int info) { lapack_zgetrs(ref trans,ref N,ref NRHS, A, ref LDA, IPIV, B, ref LDB, ref info); } public void dpotrs(char uplo, int n, int nrhs, double[] A, int lda, double[] B, int ldb, ref int info) { lapack_dpotrs(ref uplo, ref n, ref nrhs, A, ref lda, B, ref ldb, ref info); } public void spotrs(char uplo, int n, int nrhs, float[] A, int lda, float[] B, int ldb, ref int info) { lapack_spotrs(ref uplo, ref n, ref nrhs, A, ref lda, B, ref ldb, ref info); } public void cpotrs(char uplo, int n, int nrhs, fcomplex[] A, int lda, fcomplex[] B, int ldb, ref int info) { lapack_cpotrs(ref uplo, ref n, ref nrhs, A, ref lda, B, ref ldb, ref info); } public void zpotrs(char uplo, int n, int nrhs, complex[] A, int lda, complex[] B, int ldb, ref int info) { lapack_zpotrs(ref uplo, ref n, ref nrhs, A, ref lda, B, ref ldb, ref info); } public void dgelsd (int m, int n, int nrhs, double [] A, int lda, double [] B, int ldb, double [] S, double RCond, ref int rank, ref int info) { double [] work = new double [1]; int [] iwork = new int[1]; int lwork = -1; /*HC:HycalpTag1*/ lapack_dgelsd (ref m, ref n, ref nrhs, A, ref lda, B, ref ldb, S,ref RCond, ref rank, work, ref lwork, iwork, ref info); if (info != 0) throw new ILArgumentException("dgelsd: invalid parameter: #" + (-info).ToString()); lwork = ILAENV(9, "dgelsd", " ",0,0,0,0); if (lwork <= 0) throw new ILArgumentException("dgelsd: unknown error determining working size liwork"); iwork = new int[lwork]; lwork = (int)work[0]; work = new double [lwork]; /*HC:HycalpTag2*/ lapack_dgelsd (ref m, ref n, ref nrhs, A, ref lda, B, ref ldb, S,ref RCond, ref rank, work, ref lwork, iwork, ref info); } #region HYCALPER AUTO GENERATED CODE public void sgelsd (int m, int n, int nrhs, float [] A, int lda, float [] B, int ldb, float [] S, float RCond, ref int rank, ref int info) { float [] work = new float [1]; int [] iwork = new int[1]; int lwork = -1; lapack_sgelsd (ref m, ref n, ref nrhs, A, ref lda, B, ref ldb, S,ref RCond, ref rank, work, ref lwork, iwork, ref info); if (info != 0) throw new ILArgumentException("dgelsd: invalid parameter: #" + (-info).ToString()); lwork = ILAENV(9, "dgelsd", " ",0,0,0,0); if (lwork <= 0) throw new ILArgumentException("dgelsd: unknown error determining working size liwork"); iwork = new int[lwork]; lwork = (int)work[0]; work = new float [lwork]; lapack_sgelsd (ref m, ref n, ref nrhs, A, ref lda, B, ref ldb, S,ref RCond, ref rank, work, ref lwork, iwork, ref info); } public void zgelsd (int m, int n, int nrhs, complex [] A, int lda, complex [] B, int ldb, double [] S, double RCond, ref int rank, ref int info) { complex [] work = new complex [1]; int [] iwork = new int[1]; int lwork = -1; double [] rwork = new double [1]; lapack_zgelsd (ref m, ref n, ref nrhs, A, ref lda, B, ref ldb, S,ref RCond, ref rank, work, ref lwork, rwork, iwork, ref info); if (info != 0) throw new ILArgumentException("dgelsd: invalid parameter: #" + (-info).ToString()); lwork = ILAENV(9, "dgelsd", " ",0,0,0,0); if (lwork <= 0) throw new ILArgumentException("dgelsd: unknown error determining working size liwork"); iwork = new int[lwork]; lwork = (int)work[0]; work = new complex [lwork]; rwork = new double [(int)rwork[0]]; lapack_zgelsd (ref m, ref n, ref nrhs, A, ref lda, B, ref ldb, S,ref RCond, ref rank, work, ref lwork, rwork, iwork, ref info); } public void cgelsd (int m, int n, int nrhs, fcomplex [] A, int lda, fcomplex [] B, int ldb, float [] S, float RCond, ref int rank, ref int info) { fcomplex [] work = new fcomplex [1]; int [] iwork = new int[1]; int lwork = -1; float [] rwork = new float [1]; lapack_cgelsd (ref m, ref n, ref nrhs, A, ref lda, B, ref ldb, S,ref RCond, ref rank, work, ref lwork, rwork, iwork, ref info); if (info != 0) throw new ILArgumentException("dgelsd: invalid parameter: #" + (-info).ToString()); lwork = ILAENV(9, "dgelsd", " ",0,0,0,0); if (lwork <= 0) throw new ILArgumentException("dgelsd: unknown error determining working size liwork"); iwork = new int[lwork]; lwork = (int)work[0]; work = new fcomplex [lwork]; rwork = new float [(int)rwork[0]]; lapack_cgelsd (ref m, ref n, ref nrhs, A, ref lda, B, ref ldb, S,ref RCond, ref rank, work, ref lwork, rwork, iwork, ref info); } #endregion HYCALPER AUTO GENERATED CODE public void dgelsy (int m, int n, int nrhs, double [] A, int lda, double [] B, int ldb, int[] JPVT0, double RCond, ref int rank, ref int info) { int lwork = -1; double [] work = new double [1]; // /*!HC:lapack_?gelsy*/ lapack_dgelsy (ref m,ref n,ref nrhs,A,ref lda,B,ref ldb,JPVT0,ref RCond,ref rank,work,ref lwork,ref info); if (info != 0) throw new ILArgumentException("?gelsy: unable to determine optimal block size. cancelling..."); lwork = (int) work[0]; work = new double [lwork]; // /*!HC:lapack_?gelsy*/ lapack_dgelsy (ref m,ref n,ref nrhs,A,ref lda,B,ref ldb,JPVT0,ref RCond,ref rank,work,ref lwork,ref info); } #region HYCALPER AUTO GENERATED CODE public void sgelsy (int m, int n, int nrhs, float [] A, int lda, float [] B, int ldb, int[] JPVT0, float RCond, ref int rank, ref int info) { int lwork = -1; float [] work = new float [1]; lapack_sgelsy (ref m,ref n,ref nrhs,A,ref lda,B,ref ldb,JPVT0,ref RCond,ref rank,work,ref lwork,ref info); if (info != 0) throw new ILArgumentException("?gelsy: unable to determine optimal block size. cancelling..."); lwork = (int) work[0]; work = new float [lwork]; lapack_sgelsy (ref m,ref n,ref nrhs,A,ref lda,B,ref ldb,JPVT0,ref RCond,ref rank,work,ref lwork,ref info); } public void zgelsy (int m, int n, int nrhs, complex [] A, int lda, complex [] B, int ldb, int[] JPVT0, double RCond, ref int rank, ref int info) { int lwork = -1; complex [] work = new complex [1]; double[] rwork = new double[1]; lapack_zgelsy (ref m,ref n,ref nrhs,A,ref lda,B,ref ldb,JPVT0,ref RCond,ref rank,work,ref lwork, rwork, ref info); if (info != 0) throw new ILArgumentException("?gelsy: unable to determine optimal block size. cancelling..."); lwork = (int) work[0]; work = new complex [lwork]; rwork = new double[lwork]; lapack_zgelsy (ref m,ref n,ref nrhs,A,ref lda,B,ref ldb,JPVT0,ref RCond,ref rank,work,ref lwork, rwork, ref info); } public void cgelsy (int m, int n, int nrhs, fcomplex [] A, int lda, fcomplex [] B, int ldb, int[] JPVT0, float RCond, ref int rank, ref int info) { int lwork = -1; fcomplex [] work = new fcomplex [1]; float[] rwork = new float[1]; lapack_cgelsy (ref m,ref n,ref nrhs,A,ref lda,B,ref ldb,JPVT0,ref RCond,ref rank,work,ref lwork, rwork, ref info); if (info != 0) throw new ILArgumentException("?gelsy: unable to determine optimal block size. cancelling..."); lwork = (int) work[0]; work = new fcomplex [lwork]; rwork = new float[lwork]; lapack_cgelsy (ref m,ref n,ref nrhs,A,ref lda,B,ref ldb,JPVT0,ref RCond,ref rank,work,ref lwork, rwork, ref info); } #endregion HYCALPER AUTO GENERATED CODE #region ?GEEVX public void dgeevx(char balance, char jobvl, char jobvr, char sense, int n, double[] A, int lda, double[] wr, double[] wi, double[] vl, int ldvl, double[] vr, int ldvr, ref int ilo, ref int ihi, double[] scale, ref double abnrm, double[] rconde, double[] rcondv, ref int info) { double [] work = new double[1]; int lwork = -1; int [] iwork = new int[2 * n - 2]; lapack_dgeevx(ref balance, ref jobvl, ref jobvr, ref sense, ref n, A, ref lda, wr, wi, vl, ref ldvl, vr, ref ldvr, ref ilo, ref ihi, scale, ref abnrm, rconde, rcondv, work, ref lwork, iwork, ref info); if (info != 0) throw new ILArgumentException("error in lapack call: ?geevx. (" + info + ")"); lwork = (int)work[0]; work = new double[lwork]; lapack_dgeevx(ref balance, ref jobvl, ref jobvr, ref sense, ref n, A, ref lda, wr, wi, vl, ref ldvl, vr, ref ldvr, ref ilo, ref ihi, scale, ref abnrm, rconde, rcondv, work, ref lwork, iwork, ref info); } public void sgeevx(char balance, char jobvl, char jobvr, char sense, int n, float[] A, int lda, float[] wr, float[] wi, float[] vl, int ldvl, float[] vr, int ldvr, ref int ilo, ref int ihi, float[] scale, ref float abnrm, float[] rconde, float[] rcondv, ref int info) { float [] work = new float[1]; int lwork = -1; int [] iwork = new int[2 * n - 2]; lapack_sgeevx(ref balance, ref jobvl, ref jobvr, ref sense, ref n, A, ref lda, wr, wi, vl, ref ldvl, vr, ref ldvr, ref ilo, ref ihi, scale, ref abnrm, rconde, rcondv, work, ref lwork, iwork, ref info); if (info != 0) throw new ILArgumentException("error in lapack call: ?geevx. (" + info + ")"); lwork = (int)work[0]; work = new float[lwork]; lapack_sgeevx(ref balance, ref jobvl, ref jobvr, ref sense, ref n, A, ref lda, wr, wi, vl, ref ldvl, vr, ref ldvr, ref ilo, ref ihi, scale, ref abnrm, rconde, rcondv, work, ref lwork, iwork, ref info); } public void cgeevx(char balance, char jobvl, char jobvr, char sense, int n, fcomplex[] A, int lda, fcomplex[] w, fcomplex[] vl, int ldvl, fcomplex[] vr, int ldvr, ref int ilo, ref int ihi, float[] scale, ref float abnrm, float[] rconde, float[] rcondv, ref int info) { fcomplex [] work = new fcomplex[1]; int lwork = -1; float[] rwork = new float[2 * n]; lapack_cgeevx(ref balance, ref jobvl, ref jobvr, ref sense, ref n, A, ref lda, w, vl, ref ldvl, vr, ref ldvr, ref ilo, ref ihi, scale, ref abnrm, rconde, rcondv, work, ref lwork, rwork, ref info); if (info != 0) throw new ILArgumentException("error in lapack call: ?geevx. (" + info + ")"); lwork = (int)work[0]; work = new fcomplex[lwork]; lapack_cgeevx(ref balance, ref jobvl, ref jobvr, ref sense, ref n, A, ref lda, w, vl, ref ldvl, vr, ref ldvr, ref ilo, ref ihi, scale, ref abnrm, rconde, rcondv, work, ref lwork, rwork, ref info); } public void zgeevx(char balance, char jobvl, char jobvr, char sense, int n, complex[] A, int lda, complex[] w, complex[] vl, int ldvl, complex[] vr, int ldvr, ref int ilo, ref int ihi, double[] scale, ref double abnrm, double[] rconde, double[] rcondv, ref int info) { complex [] work = new complex[1]; int lwork = -1; double[] rwork = new double[2 * n]; lapack_zgeevx(ref balance, ref jobvl, ref jobvr, ref sense, ref n, A, ref lda, w, vl, ref ldvl, vr, ref ldvr, ref ilo, ref ihi, scale, ref abnrm, rconde, rcondv, work, ref lwork, rwork, ref info); if (info != 0) throw new ILArgumentException("error in lapack call: ?geevx. (" + info + ")"); lwork = (int)work[0]; work = new complex[lwork]; lapack_zgeevx(ref balance, ref jobvl, ref jobvr, ref sense, ref n, A, ref lda, w, vl, ref ldvl, vr, ref ldvr, ref ilo, ref ihi, scale, ref abnrm, rconde, rcondv, work, ref lwork, rwork, ref info); } #endregion #region ?SYEVR public void dsyevr(char jobz, char range, char uplo, int n, double[] A, int lda, double vl, double vu, int il, int iu, double abstol, ref int m, double[] w, double[] z, int ldz, int[] isuppz, ref int info) { double [] work = new double[1]; int lwork = -1, liwork = -1; int [] iwork = new int[1]; lapack_dsyevr(ref jobz,ref range,ref uplo,ref n, A,ref lda, ref vl, ref vu, ref il, ref iu, ref abstol,ref m, w,z,ref ldz,isuppz, work,ref lwork,iwork,ref liwork,ref info); if (info != 0) { throw new ILArgumentException("?syevr: error returned from lapack: " + info); } lwork = (int)work[0]; liwork = (int) iwork[0]; //work = ILMemoryPool.Pool.New(lwork); //iwork = ILMemoryPool.Pool.New(liwork); work = new double[lwork]; iwork = new int[liwork]; lapack_dsyevr(ref jobz,ref range,ref uplo,ref n, A,ref lda, ref vl, ref vu, ref il, ref iu, ref abstol,ref m, w,z,ref ldz,isuppz, work,ref lwork,iwork,ref liwork,ref info); //ILMemoryPool.Pool.RegisterObject(work); //ILMemoryPool.Pool.RegisterObject(iwork); if (info != 0) { throw new ILArgumentException("?syevr: error returned from lapack: " + info); } } public void ssyevr(char jobz, char range, char uplo, int n, float[] A, int lda, float vl, float vu, int il, int iu, float abstol, ref int m, float[] w, float[] z, int ldz, int[] isuppz, ref int info) { float [] work = new float[1]; int lwork = -1, liwork = -1; int [] iwork = new int[1]; lapack_ssyevr(ref jobz,ref range,ref uplo,ref n, A,ref lda, ref vl, ref vu, ref il, ref iu, ref abstol,ref m, w,z,ref ldz,isuppz, work,ref lwork,iwork,ref liwork,ref info); if (info != 0) { throw new ILArgumentException("?syevr: error returned from lapack: " + info); } lwork = (int)work[0]; work = ILMemoryPool.Pool.New(lwork); liwork = (int) iwork[0]; iwork = ILMemoryPool.Pool.New(liwork); lapack_ssyevr(ref jobz,ref range,ref uplo,ref n, A,ref lda, ref vl, ref vu, ref il, ref iu, ref abstol,ref m, w,z,ref ldz,isuppz, work,ref lwork,iwork,ref liwork,ref info); ILMemoryPool.Pool.Free(work); ILMemoryPool.Pool.Free(iwork); if (info != 0) { throw new ILArgumentException("?syevr: error returned from lapack: " + info); } } public void cheevr(char jobz, char range, char uplo, int n, fcomplex[] A, int lda, float vl, float vu, int il, int iu, float abstol, ref int m, float[] w, fcomplex[] z, int ldz, int[] isuppz, ref int info) { fcomplex[] work = new fcomplex[1]; float [] rwork = new float[1]; int [] iwork = new int[1]; int lrwork = -1, liwork = -1, lwork = -1; lapack_cheevr(ref jobz,ref range,ref uplo,ref n, A,ref lda, ref vl, ref vu, ref il, ref iu, ref abstol,ref m, w,z,ref ldz,isuppz, work,ref lwork,rwork,ref lrwork, iwork,ref liwork,ref info); if (info != 0) { throw new ILArgumentException("?chevr: error returned from lapack: " + info); } lrwork = (int)rwork[0]; rwork = ILMemoryPool.Pool.New(lrwork); lwork = (int) work[0]; work = ILMemoryPool.Pool.New(lwork); liwork = (int) iwork[0]; iwork = ILMemoryPool.Pool.New(liwork); lapack_cheevr(ref jobz,ref range,ref uplo,ref n, A,ref lda, ref vl, ref vu, ref il, ref iu, ref abstol,ref m, w,z,ref ldz,isuppz, work,ref lwork,rwork,ref lrwork,iwork,ref liwork,ref info); ILMemoryPool.Pool.Free(rwork); ILMemoryPool.Pool.Free(work); ILMemoryPool.Pool.Free(iwork); if (info != 0) { throw new ILArgumentException("?syevr: error returned from lapack: " + info); } } public void zheevr(char jobz, char range, char uplo, int n, complex[] A, int lda, double vl, double vu, int il, int iu, double abstol, ref int m, double[] w, complex[] z, int ldz, int[] isuppz, ref int info) { complex[] work = new complex[1]; double [] rwork = new double[1]; int [] iwork = new int[1]; int lrwork = -1, liwork = -1, lwork = -1; lapack_zheevr(ref jobz,ref range,ref uplo,ref n, A,ref lda, ref vl, ref vu, ref il, ref iu, ref abstol,ref m, w,z,ref ldz,isuppz, work,ref lwork,rwork,ref lrwork, iwork,ref liwork,ref info); if (info != 0) { throw new ILArgumentException("?heevr: error returned from lapack: " + info); } lrwork = (int)rwork[0]; rwork = ILMemoryPool.Pool.New(lrwork); lwork = (int) work[0]; work = ILMemoryPool.Pool.New(lwork); liwork = (int) iwork[0]; iwork = ILMemoryPool.Pool.New(liwork); lapack_zheevr(ref jobz,ref range,ref uplo,ref n, A,ref lda, ref vl, ref vu, ref il, ref iu, ref abstol,ref m, w,z,ref ldz,isuppz, work,ref lwork,rwork,ref lrwork,iwork,ref liwork,ref info); ILMemoryPool.Pool.Free(rwork); ILMemoryPool.Pool.Free(work); ILMemoryPool.Pool.Free(iwork); if (info != 0) { throw new ILArgumentException("?syevr: error returned from lapack: " + info); } } #endregion #region ?[he/sy]gv - generalized eigenproblem public void dsygv (int itype, char jobz, char uplo, int n, double [] A, int lda, double [] B, int ldb, double [] w, ref int info) { // query workspace int lwork = -1; double [] work = new double[1] {0.0}; lapack_dsygv (ref itype,ref jobz,ref uplo,ref n, A,ref lda, B,ref ldb, w, work, ref lwork, ref info); if (info != 0 || work[0] <= 0.0) return; // create temporary array(s) lwork = (int) work[0]; work = ILMemoryPool.Pool.New(lwork); lapack_dsygv (ref itype,ref jobz,ref uplo,ref n, A,ref lda, B,ref ldb, w, work, ref lwork, ref info); ILMemoryPool.Pool.Free(work); } public void ssygv (int itype, char jobz, char uplo, int n, float [] A, int lda, float [] B, int ldb, float [] w, ref int info) { // query workspace int lwork = -1; float [] work = new float[1] {0.0f}; lapack_ssygv (ref itype,ref jobz,ref uplo,ref n, A,ref lda, B,ref ldb, w, work, ref lwork, ref info); if (info != 0 || work[0] <= 0.0) return; // create temporary array(s) lwork = (int) work[0]; work = ILMemoryPool.Pool.New(lwork); lapack_ssygv (ref itype,ref jobz,ref uplo,ref n, A,ref lda, B,ref ldb, w, work, ref lwork, ref info); ILMemoryPool.Pool.Free(work); } public void chegv (int itype, char jobz, char uplo, int n, fcomplex[] A, int lda, fcomplex[] B, int ldb, float [] w, ref int info) { int lwork = (int)Math.Max(1,2*n-1); bool dummy; float [] rwork = ILMemoryPool.Pool.New(Math.Max(1,3*n-2),false, out dummy); fcomplex [] work = ILMemoryPool.Pool.New(lwork,false, out dummy); // NOTE: this implementation of LAPACK (netlib.org, not optimized) seems to overwrite B even if only a workspace query is assumed. // therefore this hack uses the (may not optimal efficiently) working size 2*n-1 for computation. lapack_chegv(ref itype,ref jobz,ref uplo,ref n, A,ref lda, B,ref ldb, w, work, ref lwork, rwork, ref info); ILMemoryPool.Pool.Free(rwork); ILMemoryPool.Pool.Free(work); } public void zhegv (int itype, char jobz, char uplo, int n, complex [] A, int lda, complex [] B, int ldb, double [] w, ref int info) { // NOTE: this implementation of LAPACK (netlib.org, not optimized) seems to overwrite B even if only a workspace query is assumed. // therefore this hack uses the (may not optimal efficiently) working size 2*n-1 for computation. int lwork = (int)Math.Max(1,2*n-1); bool dummy; double [] rwork = ILMemoryPool.Pool.New(Math.Max(1,3*n-2)); complex [] work = ILMemoryPool.Pool.New(lwork,false, out dummy); lapack_zhegv(ref itype,ref jobz,ref uplo,ref n, A,ref lda, B,ref ldb, w, work, ref lwork, rwork, ref info); ILMemoryPool.Pool.Free(rwork); ILMemoryPool.Pool.Free(work); } #endregion #endregion } }