// f01ve Example Program Text // C# version, NAG Copyright 2012 using System; using NagLibrary; using System.Globalization; namespace NagDotNetExamples { public class F01VEE { static string datafile = "ExampleData/f01vee.d"; // Specify data file as a command line argument. // It defaults to the named file specified below otherwise. static void Main(String[] args) { if (args.Length == 1) { datafile = args[0]; } StartExample(); } public static void StartExample() { try { DataReader sr = new DataReader(datafile); // f01ve Example Program Text // Mark 25 Release. NAG Copyright 2013. int inc1 = 1, indent = 0, ncols = 80; string diag = "N", intlabel = "I", matrix = "G", nolabel = "N"; string form = "F5.2"; int i=0, ifail=0, info=0, lda=0, lenarf=0, n=0; string title=""; string transr="", uplo=""; double[,] a; double[] arf; string[] clabs = new string[1]; string[] rlabs = new string[1]; Console.WriteLine(" {0}","F01VEF Example Program Results"); // Skip heading in data file sr.Reset(); Console.WriteLine(""); sr.Reset(); n = int.Parse(sr.Next()); uplo = sr.Next(); transr = sr.Next(); lda = (int)n; a = new double[lda, n]; lenarf = n*(n+1)/2; arf = new double[lenarf]; // allocate sr.Reset(); // Read a triangular matrix of order n into array A for (i = 1 ; i <= n ; i++) { for (int j = 1 ; j <= n ; j++) { a [ i - 1 , j - 1] = double.Parse(sr.Next(), CultureInfo.InvariantCulture); } } // Print the unpacked array A title = "Unpacked Matrix A: "; ifail = 0; X04.x04cb(uplo, diag, n, n, a, form, title, intlabel, rlabs, intlabel, clabs, ncols, indent, out ifail); Console.WriteLine(""); // Convert to Rectangular Full Packed form info = 0; // The NAG name equivalent of dtrttf is f01vef F01.f01ve(transr, uplo, n, a, arf, out info); // Print the Rectangular Full Packed array title = "RFP Packed Array ARF:"; ifail = 0; X04.x04cb(matrix, diag, lenarf, inc1, arf, form, title, intlabel, rlabs, nolabel, clabs, ncols, indent, out ifail); } catch (Exception e) { Console.WriteLine(e.Message); Console.WriteLine( "Exception Raised"); } } } }