17 type (hecmwST_matrix) :: hecMAT
18 type (fstrST_matrix_contact_lagrange) :: fstrMAT
28 type (hecmwST_local_mesh) :: hecMESH
29 type (hecmwST_matrix) :: hecMAT
30 type (fstrST_matrix_contact_lagrange) :: fstrMAT
31 integer (kind=4) :: ntdf, ilag_sta
32 integer (kind=4) :: numNon0
33 integer (kind=4) :: ierr, nprocs, myrank
35 real(kind=8), allocatable :: b(:)
37 call hecmw_mat_dump(hecmat, hecmesh)
41 ntdf = hecmat%NP*hecmat%NDOF + fstrmat%num_lagrange
42 ilag_sta = hecmat%NP*hecmat%NDOF + 1
43 numnon0 = hecmat%NPU*hecmat%NDOF**2+hecmat%NP*hecmat%NDOF*(ntdf+1)/2 &
44 + (fstrmat%numU_lagrange)*hecmat%NDOF+fstrmat%num_lagrange
46 allocate(b(
size(hecmat%B)))
55 call hecmw_mat_dump_solution(hecmat)
subroutine, public hecmw_solve_direct_serial_lag(nrows, ilag_sta, nttbr, pointers, indices, values, b)
subroutine solve_lineq_serial_lag_hecmw(hecMESH, hecMAT, fstrMAT)
subroutine solve_lineq_serial_lag_hecmw_init(hecMAT, fstrMAT, is_sym)