Ignore:
Timestamp:
Nov 18, 2004, 10:53:59 PM (20 years ago)
Author:
steve
Message:

Working on sparse. Sparse_CSR * vertor works, Sparse_CSR * matrix doesn't.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • inundation/ga/storm_surge/pyvolution/test_cg_solve.py

    r586 r594  
    8181        """Standard 2d laplacian"""
    8282       
    83         n = 50
    84         m = 50
     83        n = 100
     84        m = 100
    8585
    8686        A = Sparse(m*n, m*n)
    … …  
    106106        #print 'finish covert'
    107107        b = A*xe
    108         x = conjugate_gradient(A,b,b,iprint=10)
     108        x = conjugate_gradient(A,b,b,iprint=20)
    109109
    110110        assert allclose(x,xe)
Note: See TracChangeset for help on using the changeset viewer.