diff --git a/ila-bound-fix.patch b/ila-bound-fix.patch
new file mode 100644
index 0000000..ffd8f48
--- /dev/null
+++ b/ila-bound-fix.patch
@@ -0,0 +1,87 @@
+From SVN:
+  r977 | james | 2011-04-22 00:00:11 +0200 (Fri, 22 Apr 2011) | 2 lines
+  fixed problems with zero matrix and scan
+
+  bug0079 :: ILA[SDCZ]LR out of bound access in the DO LOOP for array A.
+  Change DO LOOP condition and add a IF inside the DO LOOP
+
+Index: SRC/ilazlr.f
+===================================================================
+--- SRC/ilazlr.f	(revision 976)
++++ SRC/ilazlr.f	(revision 977)
+@@ -53,9 +53,12 @@
+          ILAZLR = 0
+          DO J = 1, N
+             I=M
+-            DO WHILE ((A(I,J).NE.ZERO).AND.(I.GE.1))
+-              I=I-1
+-            ENDDO         
++            DO WHILE ((A(I,J).EQ.ZERO).AND.(I.GT.1))
++               I=I-1
++            ENDDO
++            IF( (I.EQ.1).AND.(A(1,J).EQ.ZERO) ) THEN
++               I = 0
++            END IF         
+             ILAZLR = MAX( ILAZLR, I )
+          END DO
+       END IF
+Index: SRC/ilaclr.f
+===================================================================
+--- SRC/ilaclr.f	(revision 976)
++++ SRC/ilaclr.f	(revision 977)
+@@ -54,9 +54,12 @@
+          ILACLR = 0
+          DO J = 1, N
+             I=M
+-            DO WHILE ((A(I,J).NE.ZERO).AND.(I.GE.1))
+-              I=I-1
+-            ENDDO         
++            DO WHILE ((A(I,J).EQ.ZERO).AND.(I.GT.1))
++               I=I-1
++            ENDDO
++            IF( (I.EQ.1).AND.(A(1,J).EQ.ZERO) ) THEN
++               I = 0
++            END IF         
+             ILACLR = MAX( ILACLR, I )
+          END DO
+       END IF
+Index: SRC/ilaslr.f
+===================================================================
+--- SRC/ilaslr.f	(revision 976)
++++ SRC/ilaslr.f	(revision 977)
+@@ -54,9 +54,12 @@
+          ILASLR = 0
+          DO J = 1, N
+             I=M
+-            DO WHILE ((A(I,J).NE.ZERO).AND.(I.GE.1))
+-              I=I-1
++            DO WHILE ((A(I,J).EQ.ZERO).AND.(I.GT.1))
++               I=I-1
+             ENDDO
++            IF( (I.EQ.1).AND.(A(1,J).EQ.ZERO) ) THEN
++               I = 0
++            END IF
+             ILASLR = MAX( ILASLR, I )
+          END DO
+       END IF
+Index: SRC/iladlr.f
+===================================================================
+--- SRC/iladlr.f	(revision 976)
++++ SRC/iladlr.f	(revision 977)
+@@ -53,10 +53,13 @@
+ *     Scan up each column tracking the last zero row seen.
+          ILADLR = 0
+          DO J = 1, N
+-            I=M
+-            DO WHILE ((A(I,J).NE.ZERO).AND.(I.GE.1))
+-              I=I-1
++            I = M
++            DO WHILE ((A(I,J).EQ.ZERO).AND.(I.GT.1))
++               I=I-1
+             ENDDO
++            IF( (I.EQ.1).AND.(A(1,J).EQ.ZERO) ) THEN
++               I = 0
++            END IF
+             ILADLR = MAX( ILADLR, I )
+          END DO
+       END IF
diff --git a/lapack-3.2.2.tar.bz2 b/lapack-3.2.2.tar.bz2
deleted file mode 100644
index 8ffe2ee..0000000
--- a/lapack-3.2.2.tar.bz2
+++ /dev/null
@@ -1,3 +0,0 @@
-version https://git-lfs.github.com/spec/v1
-oid sha256:e5f4ca0fca10a01636c698f6d390db62a826106a384f7ab7ebbdd3437e69a8ae
-size 3007563
diff --git a/lapack-3.3.1.tar.bz2 b/lapack-3.3.1.tar.bz2
new file mode 100644
index 0000000..c8df82e
--- /dev/null
+++ b/lapack-3.3.1.tar.bz2
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:8813ba185727a551f7c2450b96febf0aed7c03d0fc94cd7da5d5b84b180eca0a
+size 3092598
diff --git a/lapack.changes b/lapack.changes
index b69e060..c7d1e56 100644
--- a/lapack.changes
+++ b/lapack.changes
@@ -1,3 +1,15 @@
+-------------------------------------------------------------------
+Thu May 12 07:12:55 UTC 2011 - burnus@net-b.de
+
+- Update to 3.3.1 [bnc#654560]
+  * New procedures to compute the complete CS decomposition of a
+    partitioned unitary matrix
+  * Level-3 BLAS symmetric indefinite solve (xSYTRS) and symmetric
+    indefinite inversion (xSYTRI)
+  * Bug fixes
+  * Changes see at http://www.netlib.org/lapack/lapack-3.3.1.html
+    and http://www.netlib.org/lapack/lapack-3.3.0.html
+
 -------------------------------------------------------------------
 Thu Aug 12 13:52:25 UTC 2010 - rguenther@novell.com
 
diff --git a/lapack.spec b/lapack.spec
index 237296c..4939182 100644
--- a/lapack.spec
+++ b/lapack.spec
@@ -21,14 +21,15 @@
 Name:           lapack
 BuildRequires:  gcc-fortran
 Summary:        Linear Algebra Package
-Version:        3.2.2
+Version:        3.3.1
 Release:        1
 License:        Public Domain, Freeware
 Group:          Development/Libraries/Parallel
 Source0:        lapack-%{version}.tar.bz2
-Source1:        manpages-3.2.0.tar.bz2
+Source1:        manpages.tar.bz2
 Source99:       lapack-rpmlintrc
 Patch1:         lapack-3.2.2.patch
+Patch2:         ila-bound-fix.patch
 Url:            http://www.netlib.org/lapack/
 Requires:       blas  = %{version}
 Requires:       liblapack3 = %{version}
@@ -150,7 +151,8 @@ Authors:
 %prep
 %setup -q
 %patch1
-tar xjf %{S:1} --strip-components=1
+%patch2
+tar xjf %{S:1}
 
 %build
 cp make.inc.example make.inc
diff --git a/manpages-3.2.0.tar.bz2 b/manpages-3.2.0.tar.bz2
deleted file mode 100644
index fdb7e5b..0000000
--- a/manpages-3.2.0.tar.bz2
+++ /dev/null
@@ -1,3 +0,0 @@
-version https://git-lfs.github.com/spec/v1
-oid sha256:dd3438ab0c9d0ff314f3394f47054bab14571f568fa9d33e3b7f410f6c11b352
-size 656985
diff --git a/manpages.tar.bz2 b/manpages.tar.bz2
new file mode 100644
index 0000000..05c4e9b
--- /dev/null
+++ b/manpages.tar.bz2
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:18ea32f5d93732ffe6e033b3f8d10a1c31ec04ecaec1443ab1859bac71acc496
+size 691541