Dear all, we hope all is well with final exams and project 5. We have been made aware of the fact that many of you this year have exams close to our deadline for project 5, December 15.

We decided thus to move the final deadline to December 18 at the devilrish time 2359 (1159pm), hoping this can be of help. That is the final deadline. 

Dear All, this is our last week and our last lecture is Thursday (no lecture on Friday, only lab). 

We will end our discussion of Partial differential equations (chapter 10 of lecture notes) with hints and practicalities related to the various variants of project 5, before we summarize what we have done this semester. We will also present some future research directions in computational science, with potential master of science thesis topics and/or courses of interest for many of you. 

Else, both Thursday and Friday we will have our regular lab sessions. These are also the last lab sessions this semester. 

It has been a great pleasure to get to know you all and we wish you all the very best,

Aram, Daniel, Halvard, Jonas, Morten and Sebastian

Small addendum to weekly FYS3150/4150 plans

Hi all, I forgot to say that for those of you interested in object orientation in C++ (of interest for those of you who plan to do the solar system variant for p5), we will have a small session from 10.15-12 (approx) at the lab tomorrow. Most likely we will move over to room F?465 at the center for computing in science education (same floor as the lab on Fridays).   If you wish to just hear about this as well, please feel free to join.  We meet at the lab and then move eventually on to another room.

Furthermore, for those of you working on other variants of project 5, if it is of interest we can suggroup discussions on Friday and next week as well.

Best wishes to you all,

Morten

Hi all and thx so much for fantastic efforts with project 4. We are getting close to the end of the semester, sadly and 15 weeks have flown away. Our last topic is partial differential equations linked up with two of the final variants of p5. Tomorrow and Friday we will discuss in depth the diffusion equation in one and two dimensions as well as Poisson's and Laplace's equations in 2-dim. Our last lecture on PDEs (wave equation in 1 and 2 dims) is next Thursday, with a summary as well. There is no lecture on Friday the 29th of November, only lab. Our last lab sessions are thus next week, Thursday (28) and Friday (29). Next Friday we will end the semester with pizza at the lab as well!

Best wishes to you all!

The fys3150/4150 gang.

Good morning everybody!

Here's a short update, with plans for this week and the final two ones.

Last week we started with ordinary differential equations and derived the Velocity Verlet algorithm.  We continue this week with the Runge-Kutta family of methods and derive the most famous one, the RK4 method and discuss also adaptive ODE methods. We will thereafter move over to partial differential equations (PDEs) with the diffusion equation first on Friday. We continue with PDEs next week as well, before we wrap up our course during the last week. Project 5 is also out, there are presently 7 variants!  The variant based on the Ising model for modeling of elections is not yet ready. Furthermore, there will be an additional project on PDEs (thanks to Hilmar). These two will be posted over the weekend.

The deadline for project 5 is set to December 15. Unfortunately we cannot extend the deadline beyond that, mainly due to grading deadlines at UiO...

Good morning everybody, just a short message this time about the plans for week 45.

On Thursday we will wrap up our Monte Carlo discussions, pointing at possible other applications and finalize the discussion on the calculation of errors.  The material on error discussions and the central limit theorem can be found in the lecture notes chapter 11.2.  After this we start with our last topic, ordinary and partial differential equations (chapters 8 and 10 of the lecture notes). We start with ordinary differential equations (ODEs) this week (most likely towards the end of the Thursday lecture) and will focus on the standard methods like the family of Runge-Kutta methods and so-called symplectic ODEs, see chapter 8 of the lecture notes. Thereafter, towards the end of next week, we start with partial differential equations and that will keep us busy till the end of the semester. 

At the l...

Dear all,

last week we discussed in detail the technicalities of project 4 and linked the Ising model with the Monte Carlo algorithm and the sampling rule for accepting or rejecting a move. The latter lead us to the Metropolis algorithm.

We will derive it this Thursday (see chapter 12.5 of the lecture notes)  and we will also continue our discussion of project 4 this Thursday and parts of Friday. Most of these technicalities are covered by chapter 13 of the lecture notes. 

Friday we will need to go back to chapter 11 and discuss some statistics which is related to our understanding of errors and how we compute them in a reliable way. This will lead us to the central limit theorem and the calculation of statistical quantities like the covariance (covered by chapter 11.2 of the lecture notes). Hopefully this can lead to a better understanding on how we compute the standard deviation from a Monte Carlo calculation.  

&nb...

Hi all and thx a million for heroic efforts with project 3!!!  We really appreciate what you have done and are very impressed. We look forward to read your projects.

Project 4 is out and we will discuss the mathematics and science behind this project during tomorrow's and Friday's lectures. During these lectures we will also derive the Metropolis algorithm, discuss Markov chains and link this to the statistical physics of project 4. The material of relevance can be found in chapters 12 and 13 of the lecture notes, as well as the slides on Markov chains and brownian motion at http://compphysics.github.io/ComputationalPhysics/doc/pub/rw/html/rw.html

and the slides on statistical physics and project 4 at ...

Hi all and welcome back to FYS3150/4150. We obviously hope that the midterm week fared well with you all. 

Here's a brief update on our plans for this week and the next two ones wrt Monte Carlo methods. 

On Thursday we will wrap up our discussions about project 3, in particular repeat some basic elements if 3a and 3b as well as discussing multidimensional integrals with Monte Carlo methods. The latter material, with a similar example relevant for 3c and 3d is covered by the lecture slides (see the last slides) at for example http://compphysics.github.io/ComputationalPhysics/doc/pub/mcint/html/mcint.html

Thereafter we will also discuss on Thursday  how to compute random numbers with the random class in C...

Dear all, since many of you are busy with various midterm exams this week (best wishes to you from us) there are no lectures this week, neither on Thursday nor on Friday. 

However, we will keep the lab open and this week there is only project work (no dedicated lectures on special topics at the lab).

Concerning project 3, we recommend looking again at the slides on 

1) Gaussian quadrature at for example http://compphysics.github.io/ComputationalPhysics/doc/pub/integrate/html/integrate-reveal.html and chapter 5 of the lectures

2) Monte Carlo integration, in particular the last part o...

Good morning everybody, here follows a short update with plans for this week.

Last week we ended our discussion of eigenvalue problems (chapter 7 of the lecture notes) and started with numerical integration based on Gaussian quadrature (chapter 5 of the lecture notes). We derived the basis for Gaussian quadrature methods with the basic equations for  finding the integration weights and points. On Thursday we will continue this discussion with examples (chapter 5) before we venture into Monte Carlo (MC) integration (chapter 11, in particular 11.4 and 11.5).  MC integration will also serve as our starting point for discussing MC methods in general. This topic will keep us busy for the whole month of October.

 ...

Hi all and please apologize for the delay in sending out the weekly digest. 

This week we will end on Thursday our discussion of eigenvalue methods (first lecture on Thursday)  with a focus on how to find the eigenvalues for a tridiagonal matrix (optional exercise in project 2)  as well as discussing an iterative method for solving eigenvalue problems. We will limit ourselves to symmetric matrices and discuss an algorithm called Lanczos' method. 

Thereafter we start with numerical integration (see slides at http://compphysics.github.io/ComputationalPhysics/doc/pub/integrate/html/integrate.html)

and we will start by reminding ourselves about some simpler integration methods like the trapezoidal rule and Simpson's rule. T...

Hi all and welcome back to fys3150/4150. Last week we started discussing eigenvalue solvers and project 2. We continue this week in particular with a detailed discussion of Jacobi's method and other parts of project 2 on Thursday, see either chapter 7 of the lecture notes or the slides on eigenvalue problems at http://compphysics.github.io/ComputationalPhysics/doc/pub/eigvalues/html/eigvalues-bs.html.

On Friday we discuss Householder's algorithm, the QR algorithm and how to find eigenvalues of a tridiagonal matrix more efficiently. We will also discuss an iterative method like Lanczos' method if we get time. ...

Dear All and welcome to a new week and thx a million for heroic efforts on project 1. Project 2 is now out and a git pull gives you the latest exciting project!  It deals with eigenvalue problems and brings us over to this week's new topic, linear algebra and eigenvalue solvers. The material is covered by chapter 7 of the lecture notes and slides at http://compphysics.github.io/ComputationalPhysics/doc/pub/eigvalues/html/eigvalues.html

We will however start Thursday's lecture with a brief summary from last week, with a short repetition on LU decomposition. We will finalize the discussion on iterative linear equation solvers and move on to interpolation methods before we start with eigenvalue methods. Iterative solvers a...

Hello everybody. We hope this week started the best possible way. 

Here follows a brief summary of last week with plans for this week: 

Lectures: last week   we discussed in detail project 1, with an emphasis on the numerics and mathematics of gaussian elimination for a tridiagonal matrix, both a general one and so-called tridiagonal Toeplitz matrix that results from the discretization of the second derivative. The details are discussed in for example chapter 6 of the lectures and the slides at ...

Dear all, in the UiO schedule there were errors concerning where which room to use for the labs.

The correct version (now also correct under the schedule link) is

Thursdays at room FV203, Department of Physics, from 1015am to 6pm

Fridays at room F?434, Department of Physics, from 1015 to 6pm.

Dear All, the first lecture is Thursday 22, 815am in Store Fysiske Auditorium, Department of Physics.

Lectures are 

Thursdays 815-10am

Fridays 815-10am

Place is Store Fysiske Auditorium

The first lecture is Thursday August 22. The last lecture is Friday November 29. The lab sessions (four different groups are)

Thursday 1015am-2pm (group 1)

Thursday 2-6pm (group 2)

Friday 1015am-2pm (group 3)

Friday 2-6pm (group 4). The first lab session is Thursday August 22.

All course material is available at the github address https://github.com/CompPhysics/ComputationalPhysics. The first lecture  starts with an introduction to the course and preparations for project 1. At the lab sessions we will present git and github for version control during the first hour of every lab session. More details about practicalities will be discus...