Changeset 7bd389a


Ignore:
Timestamp:
12/06/10 10:52:46 (2 years ago)
Author:
Oleg Batrashev <ogbash@…>
Branches:
master, external, fix-prolong, refactor, refactor-subsolvers
Children:
d6a3aeb
Parents:
52a2f01
git-author:
Oleg Batrashev <ogbash@…> (12/01/10 16:51:51)
git-committer:
Oleg Batrashev <ogbash@…> (12/06/10 10:52:46)
Message:

Remove seemingly unnecessary case from subsolvers.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • src/solvers/subsolvers.F90

    r52a2f01 r7bd389a  
    250250        allocate(A%subd(A%fullaggr%nagr+1)) 
    251251        if (present(A_interf_).and.A_interf_%nnz>0) then 
    252           if (A%arrange_type/=D_SpMtx_ARRNG_ROWS) then 
    253             write (stream,*) 'Arranging A to row storage format!' 
    254             call SpMtx_arrange(A,D_SpMtx_ARRNG_ROWS,sort=.true.) 
    255           endif 
    256           if (numprocs>1.and.A_interf_%arrange_type/=D_SpMtx_ARRNG_ROWS) then 
    257             write (stream,*) 'Arranging A_interf_ to row storage format!' 
    258             call SpMtx_arrange(A_interf_,D_SpMtx_ARRNG_ROWS,sort=.true.) 
    259           endif 
     252          !if (A%arrange_type/=D_SpMtx_ARRNG_ROWS) then 
     253          !  write (stream,*) 'Arranging A to row storage format!' 
     254          !  call SpMtx_arrange(A,D_SpMtx_ARRNG_ROWS,sort=.true.) 
     255          !endif 
     256          !if (numprocs>1.and.A_interf_%arrange_type/=D_SpMtx_ARRNG_ROWS) then 
     257          !  write (stream,*) 'Arranging A_interf_ to row storage format!' 
     258          !  call SpMtx_arrange(A_interf_,D_SpMtx_ARRNG_ROWS,sort=.true.) 
     259          !endif 
    260260          if (numprocs==1) then  
    261261            A_interf_%nnz=0 
     
    277277                 val_interf=A_interf_%val) 
    278278        else 
    279           if (A%arrange_type==D_SpMtx_ARRNG_ROWS) then 
    280 !write(stream,*)'######## calling 1  multi_subsolve',A%mtx_bbe(2,2) 
    281 !write(stream,*)'M%ninner,M%nlf:',M%ninner,M%nlf 
    282 !write(stream,*)'A to be solved with is:' 
    283 !call SpMtx_printRaw(A=A,startnz=1,endnz=A%mtx_bbe(2,2)) 
    284 !call doug_abort('testing...',5342) 
    285             call multi_subsolve(               & 
    286                    nids=A%nsubsolves,          & 
    287                    ids=A%subsolve_ids,         & 
    288                    sol=sol,                    & 
    289                    rhs=rhs,                    & 
    290                    subd=A%subd,                & 
    291                    nfreds=A%nrows,             & 
    292                    nnz=A%mtx_bbe(2,2),         & 
    293                    indi=A%indi,                & 
    294                    indj=A%indj,                & 
    295                    val=A%val) 
    296           else 
    297             call multi_subsolve(               & 
    298                    nids=A%nsubsolves,          & 
    299                    ids=A%subsolve_ids,         & 
    300                    sol=sol,                    & 
    301                    rhs=rhs,                    & 
    302                    subd=A%subd,                & 
    303                    nfreds=A%nrows,             & 
    304                    nnz=A%mtx_bbe(2,2),         & 
    305                    indi=A%indi,                & 
    306                    indj=A%indj,                & 
    307                    val=A%val,                  & 
    308                    nagr1=A%fullaggr%nagr,      & 
    309                    starts1=A%fullaggr%starts,  & 
    310                    nodes1=A%fullaggr%nodes) 
    311           endif 
     279           call multi_subsolve(             & 
     280                nids=A%nsubsolves,          & 
     281                ids=A%subsolve_ids,         & 
     282                sol=sol,                    & 
     283                rhs=rhs,                    & 
     284                subd=A%subd,                & 
     285                nfreds=A%nrows,             & 
     286                nnz=A%mtx_bbe(2,2),         & 
     287                indi=A%indi,                & 
     288                indj=A%indj,                & 
     289                val=A%val) 
    312290        endif 
    313291      endif !} 
Note: See TracChangeset for help on using the changeset viewer.