Theses and Dissertations
Permanent URI for this collection
Browse
Recent Submissions
Item Open Access A local characterization of domino evacuation-shuffling(Colorado State University. Libraries, 2024) McCann, Jacob, author; Gillespie, Maria, advisor; Peterson, Christopher, committee member; Huang, Dongzhou, committee memberWe consider linear intersection problems in the Grassmanian (the space of k-dimensional subspaces of Cn), where the dimension of the intersection is 2. These spaces are called Schubert surfaces. We build of the previous work of Speyer [1] and Gillespie and Levinson [2]. Speyer showed there is a combinatorial interpretation for what happens to fibers of Schubert intersections above a "wall crossing", where marked points corresponding to the coordinates of partitions coincide. Building off Speyer's work, Levinson showed there is a combinatorial operation associated with the monodromy operator on Schubert curves, involving rectification, promotion, and shuffling of Littlewood-Richardson Young Tableaux, which overall is christened evacuation-shuffling. Gillespie and Levinson [2] further developed a localization of the evacuation-shuffling algorithm for Schubert curves. We fully develop a local description of the monodromy operator on certain classes of curves embedded inside Schubert surfaces [3].Item Open Access Data-driven methods for compact modeling of stochastic processes(Colorado State University. Libraries, 2024) Johnson, Mats S., author; Aristoff, David, advisor; Cheney, Margaret, committee member; Pinaud, Olivier, committee member; Krapf, Diego, committee memberStochastic dynamics are prevalent throughout many scientific disciplines where finding useful compact models is an ongoing pursuit. However, the simulations involved are often high-dimensional, complex problems necessitating vast amounts of data. This thesis addresses two approaches for handling such complications, coarse graining and neural networks. First, by combining Markov renewal processes with Mori-Zwanzig theory, coarse graining error can be eliminated when modeling the transition probabilities of the system. Second, instead of explicitly defining the low-dimensional approximation, using kernel approximations and a scaling matrix the appropriate subspace is uncovered through iteration. The algorithm, named the Fast Committor Machine, applies the recent Recursive Feature Machine of Radhakrishnan et al. to the committor problem using randomized numerical linear algebra. Both projects outline practical data-driven methods for estimating quantities of interest in stochastic processes that are tunable with only a few hyperparameters. The success of these methods is demonstrated numerically against standard methods on the biomolecule alanine dipeptide.Item Open Access Pseudostable Hodge classes(Colorado State University. Libraries, 2024) Williams, Matthew M., author; Cavalieri, Renzo, advisor; Gillespie, Maria, committee member; Peterson, Chris, committee member; Tucker, Dustin, committee memberWe study the relationship between Hodge classes on moduli spaces of pseudostable and stable curves given by the contraction morphism T. While Mumford's relation does not hold in the pseudostable case, we show that one can express the (pullback via T of the) Chern classes of E ⊕ E^∨ solely in terms of strata and ψ classes. We organize the combinatorial structure of the pullback of products of two pseudostable λ classes and obtain an explicit comparison of arbitrary pseudostable and stable quadratic Hodge integrals, as well as certain families of cubic and higher degree pseudostable Hodge integrals.Item Open Access A quantum H*(T)-module via quasimap invariants(Colorado State University. Libraries, 2024) Lee, Jae Hwang, author; Shoemaker, Mark, advisor; Cavalieri, Renzo, advisor; Gillespie, Maria, committee member; Peterson, Christopher, committee member; Hulpke, Alexander, committee member; Chen, Hua, committee memberFor X a smooth projective variety, the quantum cohomology ring QH*(X) is a deformation of the usual cohomology ring H*(X), where the product structure is modified to incorporate quantum corrections. These correction terms are defined using Gromov-Witten invariants. When X is toric with geometric quotient description V//T, the cohomology ring H*(V//T) also has the structure of a H*(T)-module. In this paper, we introduce a new deformation of the cohomology of X using quasimap invariants with a light point. This defines a quantum H*(T)-module structure on H*(X) through a modified version of the WDVV equations. We explicitly compute this structure for the Hirzebruch surface of type 2. We conjecture that this new quantum module structure is isomorphic to the natural module structure of the Batyrev ring for a semipositive toric variety.Item Open Access Sparse matrix varieties, Daubechies spaces, and good compression regions of Grassmann manifolds(Colorado State University. Libraries, 2024) Collery, Brian, author; Peterson, Chris, advisor; Shonkwiler, Clayton, advisor; Cavalieri, Renzo, committee member; Kirby, Michael, committee member; Pouchet, Louis-Nöel, committee memberThe Grassmann manifold Gr(k, n) is a geometric object whose points parameterize k dimensional subspaces of Rn. The flag manifold is a generalization in that its points parameterize flags of vector spaces in Rn. This thesis concerns applications of the geometry of the Grassmann and flag manifolds, with an emphasis on image compression. As a motivating example, the discrete versions of Daubechies wavelets generate distinguished n-dimensional subspaces of R2n that can be considered as distinguished points on Gr(n, 2n). We show that geodesic paths between "Daubechies points" parameterize families of "good" image compression matrices. Furthermore, we show that these paths lie on a distinguished Schubert cell in the Grassmannian. Inspired by the structure of Daubechies wavelets, we define and explore sparse matrix varieties as a generalization. Keeping in that theme, we are interested in understanding geometric considerations that constrain the "good" compression region of a Grassmann manifold.Item Open Access Investigating individually expressed motives and collectively generated goals for equity-oriented reform in undergraduate mathematics education(Colorado State University. Libraries, 2024) Tremaine, Rachel, author; Hagman, Jess Ellis, advisor; Arnold, Elizabeth, committee member; Miranda, Rick, committee member; Basile, Vincent, committee member; DiGregorio, Gaye, committee memberSupporting diversity, equity, and inclusion (DEI) is an explicitly stated goal of many mathematics departments across the country, and addressing ongoing disparities in outcomes and experiences within undergraduate mathematics is a shared responsibility among undergraduate mathematics community members. Despite the prevalence of ideological, political, and contextual barriers to equity-oriented action within undergraduate mathematics spaces, many community members can and do take a responsive stance toward enhancing DEI within their department and at their institution. Understanding how mathematics faculty members, administrators, and students are personally motivated to take up work toward these aims within their own mathematics departments is paramount in ensuring that such work continues. In this dissertation I present two investigations which draw on cultural historical activity theory (CHAT) as a conceptual and theoretical lens. In the first investigation, I analyze the motives of 30 undergraduate mathematics community members (five administrators, 17 faculty members, and eight students) across three institutions to understand their reasoning for participation in an intradepartmental community focused on creating transformative, equity-oriented change within introductory mathematics courses. A reflexive thematic analysis of journal entries and individual interviews with participants resulted in five themes which motivated participation in collaborative equity reform within their mathematics department: a relational motive, a self-improvement motive, a student experience motive, an influence motive, and a values to action motive. With these themes in mind, I then consider how a Networked Improvement Community (NIC) at one institution developed a shared object for their work through a CHAT lens, highlighting what rules, communities, subjects, artifacts, and divisions of labor proved salient to this development. The prevalence and pervasiveness of self-interests, identity neutrality, and paternalism are critically discussed within the context of these investigations, and I build on existing literature to produce recommendations for disrupting such ideologies to produce transformative change in undergraduate mathematics environments. Among these recommendations are the need for critical engagement to see beyond self-interest in the context of one's own reform work, and the need for collaborative reform groups to not only position students as experts on their own experiences, but to also conceptualize instructors as novices on student experiences. I conclude with a discussion of future work supporting continued theorizing of the link between individually expressed motives and collectively generated goals in undergraduate mathematics reform efforts.Item Open Access Cracking open the black box: a geometric and topological analysis of neural networks(Colorado State University. Libraries, 2024) Cole, Christina, author; Kirby, Michael, advisor; Peterson, Chris, advisor; Cheney, Margaret, committee member; Draper, Bruce, committee memberDeep learning is a subfield of machine learning that has exploded in recent years in terms of publications and commercial consumption. Despite their increasing prevalence in performing high-risk tasks, deep learning algorithms have outpaced our understanding of them. In this work, we hone in on neural networks, the backbone of deep learning, and reduce them to their scaffolding defined by polyhedral decompositions. With these decompositions explicitly defined for low-dimensional examples, we utilize novel visualization techniques to build a geometric and topological understanding of them. From there, we develop methods of implicitly accessing neural networks' polyhedral skeletons, which provide substantial computational and memory savings compared to those requiring explicit access. While much of the related work using neural network polyhedral decompositions is limited to toy models and datasets, the savings provided by our method allow us to use state-of-the-art neural networks and datasets in our analyses. Our experiments alone demonstrate the viability of a polyhedral view of neural networks and our results show its usefulness. More specifically, we show that the geometry that a polyhedral decomposition imposes on its neural network's domain contains signals that distinguish between original and adversarial images. We conclude our work with suggested future directions. Therefore, we (1) contribute toward closing the gap between our use of neural networks and our understanding of them through geometric and topological analyses and (2) outline avenues for extensions upon this work.Item Open Access Number of 4-cycles of the genus 2 superspecial isogeny graph(Colorado State University. Libraries, 2024) Sworski, Vladimir P., author; Pries, Rachel, advisor; Hulpke, Alexander, committee member; Rajopadhye, Sanjay, committee member; Shoemaker, Mark, committee memberThe genus 2 superspecial degree-2 isogeny graph over a finite field of size p2 is a network graph whose vertices are constructed from genus 2 superspecial curves and whose edges are the degree 2 isogenies between them. Flynn and Ti discovered 4-cycles in the graph, which pose problems for applications in cryptography. Florit and Smith constructed an atlas which describes what the neighborhood of each vertex looks like. We wrote a program in SageMath that can calculate neighborhoods of these graphs for small primes. Much of our work is motivated by these computations. We examine the prevalence of 4-cycles in the graph and, motivated by work of Arpin, et al. in the genus 1 situation, in the subgraph called the spine. We calculate the number of 4-cycles that pass through vertices of 12 of the 14 kinds possible. This also resulted in constructing the neighborhood of all vertices two steps or fewer away for three special types of curves. We also establish conjectures about the number of vertices and cycles in small neighborhoods of the spine.Item Open Access Persistence and simplicial metric thickenings(Colorado State University. Libraries, 2024) Moy, Michael, author; Adams, Henry, advisor; Patel, Amit, committee member; Peterson, Christopher, committee member; Ben-Hur, Asa, committee memberThis dissertation examines the theory of one-dimensional persistence with an emphasis on simplicial metric thickenings and studies two particular filtrations of simplicial metric thickenings in detail. It gives self-contained proofs of foundational results on one-parameter persistence modules of vector spaces, including interval decomposability, existence of persistence diagrams and barcodes, and the isometry theorem. These results are applied to prove the stability of persistent homology for sublevel set filtrations, simplicial complexes, and simplicial metric thickenings. The filtrations of simplicial metric thickenings studied in detail are the Vietoris–Rips and anti-Vietoris–Rips metric thickenings of the circle. The study of the Vietoris–Rips metric thickenings is motivated by persistent homology and its use in applied topology, and it builds on previous work on their simplicial complex counterparts. On the other hand, the study of the anti-Vietoris–Rips metric thickenings is motivated by their connections to graph colorings. In both cases, the homotopy types of these spaces are shown to be odd-dimensional spheres, with dimensions depending on the scale parameters.Item Open Access Relative oriented class groups of quadratic extensions(Colorado State University. Libraries, 2024) O'Connor, Kelly A., author; Pries, Rachel, advisor; Achter, Jeffrey, committee member; Shoemaker, Mark, committee member; Rugenstein, Maria, committee memberIn 2018 Zemková defined relative oriented class groups associated to quadratic extensions of number fields L/K, extending work of Bhargava concerning composition laws for binary quadratic forms over number fields of higher degree. This work generalized the classical correspondence between ideal classes of quadratic orders and classes of integral binary quadratic forms to any base number field of narrow class number 1. Zemková explicitly computed these relative oriented class groups for quadratic extensions of the rationals. We consider extended versions of this work and develop general strategies to compute relative oriented class groups for quadratic extensions of higher degree number fields by way of the action of Gal(K/Q) on the set of real embeddings of K. We also investigate the binary quadratic forms side of Zemková's bijection and determine conditions for representability of elements of K. Another project comprising work done jointly with Lian Duan, Ning Ma, and Xiyuan Wang is included in this thesis. Our project investigates a principal version of the Chebotarev density theorem, a famous theorem in algebraic number theory which describes the splitting of primes in number field extensions. We provide an overview of the formulation of the principal density and describe its connection to the splitting behavior of the Hilbert exact sequence.Item Open Access Counting isogeny classes of Drinfeld modules over finite fields via Frobenius distributions(Colorado State University. Libraries, 2024) Bray, Amie M., author; Achter, Jeffrey, advisor; Gillespie, Maria, committee member; Hulpke, Alexander, committee member; Pallickara, Shrideep, committee member; Pries, Rachel, committee memberClassically, the size of an isogeny class of an elliptic curve -- or more generally, a principally polarized abelian variety -- over a finite field is given by a suitable class number. Gekeler expressed the size of an isogeny class of an elliptic curve over a prime field in terms of a product over all primes of local density functions. These local density functions are what one might expect given a random matrix heuristic. In his proof, Gekeler shows that the product of these factors gives the size of an isogeny class by appealing to class numbers of imaginary quadratic orders. Achter, Altug, Garcia, and Gordon generalized Gekeler's product formula to higher dimensional abelian varieties over prime power fields without the calculation of class numbers. Their proof uses the formula of Langlands and Kottwitz that expresses the size of an isogeny class in terms of adelic orbital integrals. This dissertation focuses on the function field analog of the same problem. Due to Laumon, one can express the size of an isogeny class of Drinfeld modules over finite fields via adelic orbital integrals. Meanwhile, Gekeler proved a product formula for rank two Drinfeld modules using a similar argument to that for elliptic curves. We generalize Gekeler's formula to higher rank Drinfeld modules by the direct comparison of Gekeler-style density functions with orbital integralsItem Open Access Compound-Gaussian-regularized inverse problems: theory, algorithms, and neural networks(Colorado State University. Libraries, 2024) Lyons, Carter, author; Cheney, Margaret, advisor; Raj, Raghu G., advisor; Azimi, Mahmood, committee member; King, Emily, committee member; Mueller, Jennifer, committee memberLinear inverse problems are frequently encountered in a variety of applications including compressive sensing, radar, sonar, medical, and tomographic imaging. Model-based and data-driven methods are two prevalent classes of approaches used to solve linear inverse problems. Model-based methods incorporate certain assumptions, such as the image prior distribution, into an iterative estimation algorithm, often, as an example, solving a regularized least squares problem. Instead, data-driven methods learn the inverse reconstruction mapping directly by training a neural network structure on actual signal and signal measurement pairs. Alternatively, algorithm unrolling, a recent approach to inverse problems, combines model-based and data-driven methods through the implementation of an iterative estimation algorithm as a deep neural network (DNN). This approach offers a vehicle to embed domain-level and algorithmic insights into the design of neural networks such that the network layers are interpretable. The performance, in reconstructed signal quality, of unrolled DNNs often exceeds that of corresponding iterative algorithms and standard DNNs while doing so in a computationally efficient fashion. In this work, we leverage algorithm unrolling to combine a powerful statistical prior, the compound Gaussian (CG) prior, with the powerful representational ability of machine learning and DNN approaches. Specifically, first we construct a novel iterative CG-regularized least squares algorithm for signal reconstruction and provide a computational theory for this algorithm. Second, using algorithm unrolling, the newly developed CG-based least squares iterative algorithm is transformed into an original DNN in a manner to facilitate the learning of the optimization landscape geometry. Third, a generalization on the newly constructed CG regularized least squares iterative algorithm is developed, theoretically analyzed, and unrolled to yield a novel state-of-the-art DNN that provides a partial learning of the prior distribution constrained to the CG class of distributions. Fourth, techniques in statistical learning theory are employed for deriving original generalization error bounds on both unrolled DNNs to substantiate theoretical guarantees of each neural network when estimating signals from linear measurements after training. Finally, ample numerical experimentation is conducted for every new CG-based iterative and DNN approach proposed in this paper. Simulation results show our methods outperform previous state-of-the-art iterative signal estimation algorithms and deep-learning-based methods, especially with limited training datasets.Item Open Access Properties of the reconstruction algorithm and associated scattering transform for admittivities in the plane(Colorado State University. Libraries, 2009) Von Herrmann, Alan, author; Mueller, Jennifer, advisorWe consider the inverse admittivity problem in dimension two. The focus of this dissertation is to develop some properties of the scattering transform Sγ(k) with γ ϵ W1,p(Ω) and to develop properties of the exponentially growing solutions to the admittivity equation. We consider the case when the potential matrix is Hermitian and the definition of the potential matrix used by Francini [Inverse Problems, 16, 2000]. These exponentially growing solutions play a role in developing a reconstruction algorithm from the Dirichlet-to-Neumann map of γ. A boundary integral equation is derived relating the Dirichlet-to-Neumann map of γ to the exponentially growing solutions to the admittivity equation.Item Open Access Radial basis functions for color conversion(Colorado State University. Libraries, 2008) Qiao, Yue, author; Kirby, Michael, advisorThe most difficult and challenging task in color printing is to reduce costs while maintaining superior quality. This dissertation proposes significant enhancements to printer color conversion techniques including accurate nonlinear models that incorporate perceptual color difference metrics, lossless gray component replacement (GCR) transformations, optimized toner saving algorithms and numerical/perceptual based gamut mapping methods. Radial Basis Functions (RBFs) combined with the Lp norm approximation with emphasis on L1, L2, and L∞ were developed for color conversion. The exchange algorithm was employed in the L∞ and L1 approximations with RBFs that satisfy the Haar condition. Both the Barrodale and Phillips (BP) algorithm for solving the dual problem and the Bartels and Conn's (BC) algorithm for solving the primal were extended to multidimensional color conversion. A new approach for lossless GCR was achieved by finding one dimensional color manifolds in the CMIYK color space using multidimensional optimization techniques. We proposed objective functions for toner savings, cost savings, etc., with no quality degradation. The color conversion with the toner/ink limitation problem was solved via both L1 and L∞ approximation algorithms in the neutral and saturated color regions respectively. The L1 algorithm was a modified Barrodale and Roberts (BR) primal algorithm with an added constraint, while the L∞ algorithm was developed based on the BP dual algorithm which extended the three-stage algorithm to a four-stage algorithm. A novel gamut mapping algorithm was developed based on the numerical model guided by a perceptual color difference model. The direction of the gamut mapping is not fixed as in other methods. The algorithm transformed a connected out-of-gamut colors to connected colors around the boundary of the device gamut. The out-of-gamut colors in a small neighborhood vary continuously and smoothly. Our results indicated that the color conversion quality was significantly improved. The lossless GCR algorithm is accurate and efficient. Both the BP and BC algorithms for solving the toner/ink limitation are able to convert colors from CIELab to CMY with any given toner/ink limitation. We foresee this research will have significant impact on the color reproduction industry.Item Open Access Mathematical methods for fluid-solid interfaces: meandering streams and sand ripples(Colorado State University. Libraries, 2008) Mertens, Keith, author; Putkaradze, Vakhtang, advisorThis thesis presents several mathematical methods for modeling free surfaces, interfaces, and fluid-solid interactions. This work is done in the context of two physical systems. In the first two sections, the focus will be to understand the the physics of streams flowing down inclined substrates. Models will be derived to investigate both steady state and dynamic meandering profiles. It will be shown that, through the right approximation techniques, many physical insights can be drawn about this system. These results include: a complete understanding of the steady states, transitions between steady states, mechanism of meandering, forces involved in meandering, and spectral scaling laws of long-time ensemble averaged meandering stream profiles. In the third section, the focus will shift to how one can model underlying physics when it becomes too complicated to address from first principles. Here, the power of symmetries and conservation laws are explored to derive an amplitude equation describing the interface between sand and water when the water is subjected to oscillatory flow. The thesis will then close by posing a novel way to study scaling laws with respect to parameters using Lie's prolongation algorithm. Through this work various tools will be combined from the fields of physics, engineering, applied and pure mathematics to develop approaches for reducing complex systems into tractable pieces which can be studied carefully.Item Open Access An adaptive algorithm for an elliptic optimization problem, and stochastic-deterministic coupling: a mathematical framework(Colorado State University. Libraries, 2008) Lee, Sheldon, author; Estep, Donald, advisor; Tavener, Simon, advisorThis dissertation consists of two parts. In the first part, we study optimization of a quantity of interest of a solution of an elliptic problem, with respect to parameters in the data using a gradient search algorithm. We use the generalized Green's function as an efficient way to compute the gradient. We analyze the effect of numerical error on a gradient search, and develop an efficient way to control these errors using a posteriori error analysis. Specifically, we devise an adaptive algorithm to refine and unrefine the finite element mesh at each step in the descent search algorithm. We give basic examples and apply this technique to a model of a healing wound. In the second part, we construct a mathematical framework for coupling atomistic models with continuum models. We first study the case of coupling two deterministic diffusive regions with a common interface. We construct a fixed point map by repeatedly solving the problems, while passing the flux in one direction and the concentration in the other direction. We examine criteria for the fixed point iteration to converge, and offer remedies such as reversing the direction of the coupling, or relaxation, for the case it does not. We then study the one dimensional case where the particles undergo a random walk on a lattice, next to a continuum region. As the atomistic region is random, this technique yields a fixed point iteration of distributions. We run numerical tests to study the long term behavior of such an iteration, and compare the results with the deterministic case. We also discuss a probability transition matrix approach, in which we assume that the boundary conditions at each iterations follow a Markov chain.Item Open Access Large-scale computational analysis of National Animal Identification System mock data, including traceback and trace forward(Colorado State University. Libraries, 2008) Ladd, Joshua, author; Burns, Patrick J., advisorCattle production is the single largest segment of U.S. agriculture. Animal disease, whether a single incident or a full-scale outbreak, can result in significantly restricted access to both foreign and domestic markets. Regaining consumer confidence is difficult. If a disease cannot be traced back to a common source, then only time can tell whether or not eradication and containment efforts have been successful. Simply "waiting it out" can result in long-term economic losses on a National scale especially when diseases which are prone to epizootic outbreaks or those with long incubation periods are involved. The United States Department of Agriculture (USDA) maintains that traceability is the key to protecting animal health and marketability: The National Animal Identification System (NAIS) is a voluntary disease traceability framework released by the USDA. Many of the efforts surrounding the development of the NAIS have encompassed the identification of livestock production and handling premises as well as individuals or herds of animals, whereas little effort has been directed toward the ultimate goal of animal traceback in 48 hours. In this dissertation, computational science is applied to the problem of animal disease traceability. In particular, a computational model is developed for the purpose of conducting large-scale traceability simulations. The model consists of two components; the first being a parallel, Monte Carlo discrete events simulator capable of generating large, NAIS-compliant, mock datasets representative of the processing requirements of actual NAIS data. The second component is a large-scale, parallel disease tracing algorithm that is mapped onto an SMP supercomputer where high-performance is achieved by adopting a hybrid parallel programming model that mixes a shared memory multi-threading model (OpenMP) with a distributed memory message passing model (MPI). The objectives of this dissertation are to characterize the computational requirements of the NAIS, identify computational platforms and programming paradigms well suited to this effort, and to identify and address computational performance bottlenecks associated with large-scale tracing algorithms.Item Open Access Modeling spatio-temporal systems with skew radial basis functions: theory, algorithms and applications(Colorado State University. Libraries, 2008) Jamshidi, Arthur (Arta) Amir, author; Kirby, Michael, advisorThe discovery of knowledge in large data sets can often be formulated as a problem in nonlinear function approximation. The inherent challenge in such an approach is that the data is often high dimensional, scattered and sparse. Given a limited number of exemplars one would like to construct models that can generalize to new regions or events. Additionally, underlying physical processes may not be stationary and the nature of the nonlinear relationships may evolve. Ideally, a good model would be adaptive and remain valid over extended regions in space and time. In this work we propose a new Radial Basis Function (RBF) algorithm for constructing nonlinear models from high-dimensional scattered data. The algorithm progresses iteratively adding a new function at each step to refine the model. The placement of the functions is driven by one or more statistical hypotheses tests that reveal geometric structure in the data when it fails. At each step the added function is fit to data contained in a spatio-temporally defined local region to determine the parameters, in particular, the scale of the local model. Unlike prior techniques for nonlinear function fitting over scattered data, the proposed method requires no ad hoc parameters and it behaves effectively like a black box. Thus, the number of basis functions required for an accurate fit is determined automatically by the algorithm. An extension of the algorithms to multivariate case, i.e., the dimension of the range of the mapping is greater or equal to two, is also carried out. This approach produces more parsimonious models by exploiting the correlation among the various range dimensions. The convergence properties of the algorithms are shown from different prospectives. To further enhance the order and conditioning of the models we introduce several new compactly supported RBFs for approximating functions in LP(Rd) via over-determined least squares. We also propose a skew-radial basis function expansion for the empirical model fitting problem to achieve more accuracy and lower model orders. This is accomplished by modulating or skewing, each RBF by an asymmetric shape function which increases the number of degrees of freedom available to fit the data. We show that if the original RBF interpolation problem is positive definite, then so is the skew-radial basis function when it is viewed as a bounded perturbation of the RBF. We illustrate the utility of the theoretic and algorithmic innovations via several applications including modeling data on manifolds, prediction of financial and chaotic time-series and prediction of the maximum wind intensity of a hurricane. In addition, the skew-radial basis functions are shown to provide good approximations to data with jumps. While the algorithms presented here are in the context of RBFs, in principle they can be employed with other methods for function approximation such as multi-layer perceptrons.Item Open Access Automorphism towers of general linear groups(Colorado State University. Libraries, 2008) Jónsdóttir, Margrét Sóley, author; Hulpke, Alexander, advisorLet G0, be a group, G 1 be the automorphism group of G0, G2 the automorphism group of G1 etc. The sequence of these groups together with the natural homomorphisms πi,i+1 : Gi → Gi+1, which take each element to the inner automorphism it induces, is called the automorphism tower of G 0. If πi,i+1 is an isomorphism for some i then the automorphism tower of G is said to terminate. For a given group it is in general not easy to say whether its automorphism tower terminates. Wielandt showed in 1939 that if G is finite with a trivial center then the automorphism tower of G will terminate in a finite number of steps. Since then, some sporadic examples of automorphism towers of finite groups have been described but no general results have been proven. In this thesis we study automorphism towers of finite groups with a non-trivial center. We look at the two extremes: (1) Groups which are center-rich. (2) Groups which have a small but non-trivial center. We show that when looking for an infinite family of groups with terminating automorphism towers the first case is unfeasible. We then turn our attention to the latter case, specifically general linear groups of dimension at least two. In odd characteristic GL(2, q) is not a split extension of the center. The first thing we do is to calculate the automorphism group of GL(2, q) for odd prime powers q. We provide explicit generators and describe the structure of Aut(GL(2, q)) in terms of well-known groups. In this case, the first automorphism group in the tower is a subdirect product of two characteristic factors. This structure is propagated through the tower and we use it to reduce the problem to studying subgroups of automorphism groups of smaller groups. We then use this structure to compute examples of automorphism towers of GL(2, q).Item Open Access A ratio ergodic theorem on Borel actions of Zd and Rd(Colorado State University. Libraries, 2009) Holt, Eric Norman, author; Rudolph, Daniel, advisorWe prove a ratio ergodic theorem for free Borel actions of Zd and Rd on a standard Borel probability space. The proof employs an extension of the Besicovitch Covering Lemma, as well as a notion of coarse dimension that originates in an upcoming paper of Hochman. Due to possible singularity of the measure, we cannot use functional analytic arguments and therefore diffuse the measure onto the orbits of the action. This diffused measure is denoted μx, and our averages are of the form 1/μx(Bn) ∫ Bn f o T-v(x)dμx. A Følner condition on the orbits of the action is shown, which is the main tool used in the proof of the ergodic theorem. Also, an extension of a known example of divergence of a ratio average is presented for which the action is both conservative and free.