Bayesian closed-skew Gaussian inversion is defined as a generalization of traditional Bayesian Gaussian inversion. Bayesian inversion is often used in seismic inversion, and the closed-skew model is able to capture the skewness in the variable of interest. Different stationary prior models are presented, but the generalization comes at a cost, simulation from high-dimensional pdfs and parameter inference from data is more complicated. An efficient algorithm to generate realizations from the high-dimensional closed-skew Gaussian distribution is presented. A full-likelihood is used for parameter estimation of stationary prior models under exponential dependence structure. The simulation algorithms and estimators are evaluated on synthetic examples. Also a closed-skew T-distribution is presented to include heavy tails in the pdf and the model is presented with some examples. In the last part the simulation algorithm, the different prior models and parameter estimators are demonstrated on real data from a well in the Sleipner Øst field. The full-likelihood estimator seems to be the best estimator for data with exponential dependence structure

Spectral discretization in space and time of the weak formulation of a partial differential equations (PDE) is studied. The exact solution to the PDE, with either Dirichlet or Neumann boundary conditions imposed, is approximated using high order polynomials. This is known as a spectral Galerkin method. The main focus of this work is the solution algorithm for the arising algebraic system of equations. A direct fast tensor-product solver is presented for the Poisson problem in a rectangular domain. We also explore the possibility of using a similar method in deformed domains, where the geometry of the domain is approximated using high order polynomials. Furthermore, time-dependent PDE's are studied. For the linear convection-diffusion equation in $mathbb{R}$ we present a tensor-product solver allowing for parallel implementation, solving $mathcal{O}(N)$ independent systems of equations. Lastly, an iterative tensor-product solver is considered for a nonlinear time-dependent PDE. For most algorithms implemented, the computational cost is $mathcal O (N^{p+1})$ floating point operations and a memory required of $mathcal O (N^{p})$ floating point numbers for $mathcal O (N^{p})$ unknowns. In this work we only consider $p=2$, but the theory is easily extended to apply in higher dimensions. Numerical results verify the expected convergence for both the iterative method and the spectral discretization. Exponential convergence is obtained when the solution and domain geometry are infinitely smooth.

Heterogeneous and parallel computing systems are increasingly appealing to high-performance computing. Among heterogeneous systems, the GPUs have become an attractive device for compute-intensive problems. Their many-core architecture, primarily customized for graphics processing, is now widely available through programming architectures that exploit parallelism in GPUs. We follow this new trend and attempt an implementation of a classical mathematical model describing incompressible single-phase fluid flow through a porous medium. The porous medium is an oil reservoir represented by means of corner-point grids. Important geological and mathematical properties of corner-point grids will be discussed. The model will also incorporate pressure- and rate-controlled wells to be used for some realistic simulations. Among the test models are the 10th SPE Comparative Solution Project Model 2. After deriving the underlying mathematical model, it will be discretised using the numerical technique of Mimetic Finite Difference methods. The heterogeneous system utilised is a desktop computer with an NVIDIA GPU, and the programming architecture to be used is CUDA, which will be described. Two different versions of the final discretised system have been implemented; a traditional way of using an assembled global stiffness sparse matrix, and a Matrix-free version, in which only the element stiffness matrices are used. The former version evaluates two GPU libraries; CUSP and THRUST. These libraries will be briefly decribed. The linear system is solved using the iterative Jacobi-preconditioned conjugate gradient method. Numerical tests on realistic and complex reservoir models shows significant performance benefits compared to corresponding CPU implementations.

This thesis evaluates risk measures for interest rate portfolios. First a model for interest rates is established: the LIBOR market model. The model is applied to Norwegian and international interest rate data and used to calculate the value of the portfolio by using Monte Carlo simulation. Estimation of volatility and correlation is discussed as well as the two risk measures value at risk and expected tail loss. The data used is analysed before the results of the backtesting evaluating the two risk measures are presented.

A block-structured approach for solving 2-dimensional finite element approximations of the Poisson equation on graphics processing units(GPUs) is developed. Linear triangular elements are used, and a matrix-free version of the conjugate gradient method is utilized for solving test problems with over 30 million elements. A speedup of 24 is achieved on a NVIDIA Tesla C1060 GPU when compared to a serial CPU version of the same solution approach, and a comparison is made with previous GPU implementations of the same problem.

I denne rapporten sammenligner vi tre varianter av Markov Chain Monte Carlo (MCMC) - simulering av Isingmodellen. Vi sammenligner enkeltnode-oppdatering, naiv blokkoppdatering og rekursiv blokkoppdatering. Vi begynner med å gi en generell introduksjon til markovfelt og Isingmodellen. Deretter viser vi det teoretiske fundamentet som MCMC-metoder hviler på. Etter det gir vi en teoretisk introduksjon til enkeltnode-oppdatering. Så gir vi en innføring i naiv blokkoppdatering som er den tradisjonelle metoden å utføre blokkoppdatering på. Deretter gir vi en tilsvarende innføring i en nylig foreslått metode for å gjennomføre blokkoppdatering, nemlig rekursiv blokkoppdatering. Blokkoppdatering er en metode som har vist seg nyttig med hensyn på miksing når vi simulerer. Med det menes at blokkoppdatering har vist seg nyttig i forhold til å utforske utfallsrommet til fordelingen vi er interessert i med færre iterasjoner enn enkeltnode-oppdatering. Problemet med naiv blokkoppdatering er imidlertid at vi raskt får en høy beregningsmengde ved at hver iterasjon tar veldig lang tid. Vi prøver også ut rekursiv blokkoppdatering. Denne tar sikte på å redusere beregningsmengden for hver iterasjon når vi utfører blokkoppdatering på et markovfelt. Vi viser så simuleringsalgoritmer og resultater. Vi har simulert Isingmodellen med enkeltnode-oppdatering, naiv blokkoppdatering og rekursiv blokkoppdatering. Det vi sammenligner er antall iterasjoner før markovfeltet konvergerer og spesielt beregningstiden pr iterasjon. Vi viser at beregningsmengden pr iterasjon øker med 91000 ganger med naiv blokkoppdatering dersom vi går fra en 3 × 3 blokk til en 5 × 5 blokk. Tilsvarende tall for rekursiv blokkoppdatering er en økning på 83 ganger fra en 3 × 3 blokk til en 5 × 5 blokk. Vi sammenligner også tiden det tar før Isingmodellen konvergerer. Når vi benytter naiv blokkoppdatering finner vi at Isingmodellen bruker 15 sekunder på å konvergere med en 3 × 3 blokk, 910 sekunder på å konvergere med en 4×4 blokk og 182000 sekunder med en 5×5 blokk. Tilsvarende tall for rekursiv blokkoppdatering er 3.74 sekunder for en 3 × 3 blokk, 72 sekunder for en 4 × 4 blokk og 141.2 sekunder for en 5×5 blokk. Når vi benytter enkeltnode-oppdatering bruker feltet 6.6 sekunder på å konvergere.

<p>In this paper we present the first efficient point counting algorithm due to Schoof, before giving a significant improvement due to Elkies. In the final section I give Satoh's algorithm which is even faster and has paved the way for the field of $p$-adic point counting.</p>

<p>Sammendrag I kapittel 1 gjøres klart hvilke forkunnskaper som er nødvendige for å kunne forstå og ha utbytte av å lese denne oppgaven. Det henvises også til litteratur som kan brukes som oppslagsverk hvis noen begreper under lesningen skulle være ukjente. I kapittel 2 klargjøres sentrale begreper som fikspunkt, kontraksjon og kompletthet. Dernest presenteres Banachs fikspunktteorem med bevis. Kapitlet peker også mot ulike anvendelser av teoremet. Kapittel 3 handler i sin helhet om den velkjente Newtons metode. I sin enkelthet demonstrerer den hvordan vi kan finne tilnærmede løsninger ved iterasjoner. Arbeid med denne metoden er derfor nyttig for forståelse av Banachs fikspunktteorem og anvendelser av teoremet i mer abstrakte sammenhenger. I kapittel 4 viser vi hvordan vi kan finne en tilnærmet løsning av en partiell differensialligning ved hjelp Jacobi-iterasjon. For å demonstrere metoden tar vi utgangspunkt i et velkjent varmeledningsproblem. Vi bruker Banachs fikspunktteorem for å vise at metoden kan brukes i vårt tilfelle. Dette leder oss inn i flere disipliner av matematikken, og vi får derved demonstrert ulike teknikker for å løse denne type problemer. I kapittel 5 dreier det seg om differensialligninger og integralligninger. Som vi skal se har Picards entydighets- og eksistensteorem nær tilknytning til Banachs fikspunktteorem. Vi vil få demonstrert hvordan metoden kan brukes til å finne tilnærmede løsninger og i noen tilfeller eksakte løsninger av ordinære differensialligninger. Vi gir også et enkelt eksempel på bruk av Banachs fikspunktteorem i forhold til en spesiell type integralligning, Fredholm-ligninger. I kapittel 6 bruker vi Banachs fikspunktteorem til å bevise det inverse funksjonsteoremet. I kapittel 7 er siktemålet å sette The Scottish Café inn i en kulturhistorisk sammenheng. Vi begynner med å skissere et historisk bakteppe hvor vi søker å gi en viss oversikt over Polens historie. Behandlingen av jøder i Polen fram til 1940 blir et hovedtema. Dernest beskriver vi historien til The Scottish Café. Vi gir også en vurdering av hvilke faktorer som hadde betydning for utviklingen av miljøet. Noe av hovedkonklusjonen er at utviklinv vi gen av de sterke matematikkmiljøene i Polen i mellomkrigstiden langt på vei kan sees på som en del av nasjonsbyggingen i Polen. Utviklingen av det spesielle og særegne miljøet i The Scottish Café kan derimot sees som en reaksjon på antisemittismen i det polske samfunnet. I kapittel 8 finner vi noen refleksjoner over det vi har arbeidet med.</p>

<p>Vi innfører klassen av kvasikonforme avbildninger i det komplekse planet, og viser en rekke fundamentale egenskaper ved disse. Herunder den målbare avbildningssatsen. Heretter anvendes disse resultatene til å vise klassifikasjonssetningen for fikserte fatoukomponenter i sin optimale form. Det vil si, vi viser riktigheten av Fatous formodning om vandrende komponenter, og eksistens av siegeldisker og hermanringer. Fokus i bevisførelsen er på anvendelser av kvasikonforme avbildninger.</p>

<p>Oppgaven handler om fouriertransformasjon av generaliserte funksjoner, med spesiell vekt på fouriertransformasjon av hyperfunksjoner. Transformasjonen på hyperfunksjoner er deretter sammenlignet med Carlemans fouriertransform, som er en av de tidlige generaliseringene av den klassiske fouriertransformen. Det er vist at begge transformene er symmetriske, det vil si at invers fouriertransform er også definert på samme rom. Videre vises det at begge transformene generaliserer både den klassiske fouriertransformen og Schwartz' fouriertransform av distribusjoner med kompakt support.</p>