// f01vb Example Program Text // C# version, NAG Copyright 2012 using System; using NagLibrary; namespace NagDotNetExamples { public class F01VBE { static string datafile = "ExampleData/f01vbe.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); // f01vb Example Program Text // Mark 25 Release. NAG Copyright 2013. int inc1 = 1, indent = 0, ncols = 80; string brac = "B", diag = "N", intlabel = "I", matrix = "G", nolabel = "N"; string form = "F5.2"; int i=0, ifail=0, info=0, lda=0, lenap=0, n=0; string title=""; string uplo=""; string[] clabs = new string[1]; string[] rlabs = new string[1]; Console.WriteLine(" {0}","F01VBF Example Program Results"); // Skip heading in data file sr.Reset(); Console.WriteLine(""); sr.Reset(); n = int.Parse(sr.Next()); uplo = sr.Next(); lda = (int)n; lenap = (n * (n + 1)) / 2; Complex[,] a = new Complex[lda, n]; Complex[] ap = new Complex[lenap]; // allocate sr.Reset(); // Read a triangular matrix of order n for (i = 1 ; i <= n ; i++) { for (int j = 1 ; j <= n ; j++) { a [ i - 1 , j - 1] = Complex.Parse(sr.Next()); } } // Print the unpacked matrix title = "Unpacked Matrix A:"; ifail = 0; X04.x04db(uplo, diag, n, n, a, brac, form, title, intlabel, rlabs, intlabel, clabs, ncols, indent, out ifail); Console.WriteLine(""); // Convert to packed vector form info = 0; // The NAG name equivalent of ztrttp is f01vbf F01.f01vb(uplo, n, a, ap, out info); // Print the packed vector title = "Packed Matrix AP: "; ifail = 0; X04.x04db(matrix, diag, lenap, inc1, ap, brac, form, title, intlabel, rlabs, nolabel, clabs, ncols, indent, out ifail); } catch (Exception e) { Console.WriteLine(e.Message); Console.WriteLine( "Exception Raised"); } } } }