File: example_solve_cg.f90

package info (click to toggle)
fortran-stdlib 0.8.1-1
  • links: PTS, VCS
  • area: main
  • in suites: sid
  • size: 34,008 kB
  • sloc: f90: 24,178; ansic: 1,244; cpp: 623; python: 119; makefile: 13
file content (17 lines) | stat: -rw-r--r-- 448 bytes parent folder | download 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
 
program example_solve_cg
    use stdlib_kinds, only: int8, dp
    use stdlib_linalg_iterative_solvers, only: stdlib_solve_cg

    real(dp) :: matrix(2,2)
    real(dp) :: x(2), rhs(2)

    matrix = reshape( [4, 1,&
                       1, 3] , [2,2])

    x   = dble( [2,1] ) !> initial guess
    rhs = dble( [1,2] ) !> rhs vector

    call stdlib_solve_cg(matrix, rhs, x, restart=.false.)
    print *, x !> solution: [0.0909, 0.6364]

end program