Free cookie consent management tool by TermsFeed Policy Generator

source: branches/HeuristicLab.Problems.GaussianProcessTuning/ILNumerics.2.14.4735.573/Functions/builtin/eye.cs @ 11316

Last change on this file since 11316 was 9102, checked in by gkronber, 12 years ago

#1967: ILNumerics source for experimentation

File size: 4.1 KB
RevLine 
[9102]1///
2///    This file is part of ILNumerics Community Edition.
3///
4///    ILNumerics Community Edition - high performance computing for applications.
5///    Copyright (C) 2006 - 2012 Haymo Kutschbach, http://ilnumerics.net
6///
7///    ILNumerics Community Edition is free software: you can redistribute it and/or modify
8///    it under the terms of the GNU General Public License version 3 as published by
9///    the Free Software Foundation.
10///
11///    ILNumerics Community Edition is distributed in the hope that it will be useful,
12///    but WITHOUT ANY WARRANTY; without even the implied warranty of
13///    MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14///    GNU General Public License for more details.
15///
16///    You should have received a copy of the GNU General Public License
17///    along with ILNumerics Community Edition. See the file License.txt in the root
18///    of your distribution package. If not, see <http://www.gnu.org/licenses/>.
19///
20///    In addition this software uses the following components and/or licenses:
21///
22///    =================================================================================
23///    The Open Toolkit Library License
24///   
25///    Copyright (c) 2006 - 2009 the Open Toolkit library.
26///   
27///    Permission is hereby granted, free of charge, to any person obtaining a copy
28///    of this software and associated documentation files (the "Software"), to deal
29///    in the Software without restriction, including without limitation the rights to
30///    use, copy, modify, merge, publish, distribute, sublicense, and/or sell copies of
31///    the Software, and to permit persons to whom the Software is furnished to do
32///    so, subject to the following conditions:
33///
34///    The above copyright notice and this permission notice shall be included in all
35///    copies or substantial portions of the Software.
36///
37///    =================================================================================
38///   
39
40using System;
41using System.Collections.Generic;
42using System.Text;
43using System.Runtime.InteropServices;
44using ILNumerics.Storage;
45using ILNumerics.Misc;
46using ILNumerics.Native;
47using ILNumerics.Exceptions;
48
49namespace ILNumerics {
50
51    public partial class ILMath {
52
53        /// <summary>
54        /// Create double matrix having unity diagonal values 
55        /// </summary>
56        /// <param name="rows">Number of rows</param>
57        /// <param name="columns">Number of columns</param>
58        /// <returns>Unity matrix (diagonal matrix) of type double</returns>
59        public static ILRetArray<double> eye(int rows, int columns) {
60            using (ILScope.Enter()) {
61                ILArray<double> ret = zeros(rows, columns);
62                int diagLen = Math.Min(rows, columns);
63                for (int i = 0; i < diagLen; i++)
64                    ret.SetValue(1.0, i, i);
65                return ret;
66            }
67        }
68        /// <summary>
69        /// Create unity matrix, arbitrary numeric type
70        /// </summary>
71        /// <param name="rows">Number of rows</param>
72        /// <param name="columns">Number of columns</param>
73        /// <returns>Unity matrix (diagonal matrix) of element type T</returns>
74        /// <typeparam name="T">Element type</typeparam>
75        /// <exception cref="ILNumerics.Exceptions.ILArgumentException">If the type specified is not supported. Supported types are: double, float, complex, fcomplex, int, long, short, byte</exception>
76        public static ILRetArray<T> eye<T>(int rows, int columns) {
77            using (ILScope.Enter()) {
78                int diagLen = Math.Min(rows, columns);
79                ILArray<T> o = ones<T>(new ILSize(1,diagLen));
80                if (isnull(o))
81                    throw new ILArgumentException("eye does not support elements of type '" + typeof(T).Name + "'");
82                ILArray<T> ret = zeros<T>(rows, columns);
83                for (int i = 0; i < diagLen; i++)
84                    ret.SetValue(o.GetValue(i), i, i);
85                return ret;
86            }
87        }
88
89    }
90}
Note: See TracBrowser for help on using the repository browser.