1*1e99f302SBen Matthews! This software is Copyright (c) 2012-2015 The Regents of the 2*1e99f302SBen Matthews! University of California. All Rights Reserved. 3*1e99f302SBen Matthews! 4*1e99f302SBen Matthews! Permission to copy and modify this software and its documentation 5*1e99f302SBen Matthews! for educational, research and non-profit purposes, without fee, 6*1e99f302SBen Matthews! and without a written agreement is hereby granted, provided that 7*1e99f302SBen Matthews! the above copyright notice, this paragraph and the following three 8*1e99f302SBen Matthews! paragraphs appear in all copies. 9*1e99f302SBen Matthews! 10*1e99f302SBen Matthews! Permission to make commercial use of this software may be obtained 11*1e99f302SBen Matthews! by contacting: 12*1e99f302SBen Matthews! 13*1e99f302SBen Matthews! Technology Transfer Office 14*1e99f302SBen Matthews! 9500 Gilman Drive, Mail Code 0910 15*1e99f302SBen Matthews! University of California 16*1e99f302SBen Matthews! La Jolla, CA 92093-0910 17*1e99f302SBen Matthews! (858) 534-5815 18*1e99f302SBen Matthews! invent@ucsd.edu 19*1e99f302SBen Matthews! 20*1e99f302SBen Matthews! This software program and documentation are copyrighted by The 21*1e99f302SBen Matthews! Regents of the University of California. The software program and 22*1e99f302SBen Matthews! documentation are supplied "as is", without any accompanying 23*1e99f302SBen Matthews! services from The Regents. The Regents does not warrant that the 24*1e99f302SBen Matthews! operation of the program will be uninterrupted or error-free. The 25*1e99f302SBen Matthews! end-user understands that the program was developed for research 26*1e99f302SBen Matthews! purposes and is advised not to rely exclusively on the program for 27*1e99f302SBen Matthews! any reason. 28*1e99f302SBen Matthews! 29*1e99f302SBen Matthews! IN NO EVENT SHALL THE UNIVERSITY OF CALIFORNIA BE LIABLE TO ANY 30*1e99f302SBen Matthews! PARTY FOR DIRECT, INDIRECT, SPECIAL, INCIDENTAL, OR CONSEQUENTIAL 31*1e99f302SBen Matthews! DAMAGES, INCLUDING LOST PROFITS, ARISING OUT OF THE USE OF THIS 32*1e99f302SBen Matthews! SOFTWARE AND ITS DOCUMENTATION, EVEN IF THE UNIVERSITY OF 33*1e99f302SBen Matthews! CALIFORNIA HAS BEEN ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. 34*1e99f302SBen Matthews! THE UNIVERSITY OF CALIFORNIA SPECIFICALLY DISCLAIMS ANY 35*1e99f302SBen Matthews! WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES 36*1e99f302SBen Matthews! OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE. THE 37*1e99f302SBen Matthews! SOFTWARE PROVIDED HEREUNDER IS ON AN "AS IS" BASIS, AND THE 38*1e99f302SBen Matthews! UNIVERSITY OF CALIFORNIA HAS NO OBLIGATIONS TO PROVIDE 39*1e99f302SBen Matthews! MAINTENANCE, SUPPORT, UPDATES, ENHANCEMENTS, OR MODIFICATIONS. 40*1e99f302SBen Matthews 41*1e99f302SBen Matthews SUBROUTINE svLS_COMMU_CREATE(commu, commi) 42*1e99f302SBen Matthews INCLUDE "svLS_STD.h" 43*1e99f302SBen Matthews TYPE(svLS_commuType), INTENT(INOUT) :: commu 44*1e99f302SBen Matthews INTEGER(KIND=MPI_ADDRESS_KIND), INTENT(IN) :: commi 45*1e99f302SBen Matthews END SUBROUTINE svLS_COMMU_CREATE 46*1e99f302SBen Matthews 47*1e99f302SBen Matthews SUBROUTINE svLS_COMMU_FREE(commu) 48*1e99f302SBen Matthews INCLUDE "svLS_STD.h" 49*1e99f302SBen Matthews TYPE(svLS_commuType), INTENT(INOUT) :: commu 50*1e99f302SBen Matthews END SUBROUTINE svLS_COMMU_FREE 51*1e99f302SBen Matthews 52*1e99f302SBen Matthews SUBROUTINE svLS_LHS_CREATE(lhs, commu, gnNo, nNo, nnz, 53*1e99f302SBen Matthews 2 gNodes, rowPtr, colPtr, nFaces) 54*1e99f302SBen Matthews INCLUDE "svLS_STD.h" 55*1e99f302SBen Matthews TYPE(svLS_lhsType), INTENT(INOUT) :: lhs 56*1e99f302SBen Matthews TYPE(svLS_commuType), INTENT(IN) :: commu 57*1e99f302SBen Matthews INTEGER, INTENT(IN) :: gnNo, nNo, nnz 58*1e99f302SBen Matthews INTEGER, INTENT(IN) :: gNodes(nNo), rowPtr(nNo+1), 59*1e99f302SBen Matthews 2 colPtr(nnz) 60*1e99f302SBen Matthews INTEGER, INTENT(IN) :: nFaces 61*1e99f302SBen Matthews END SUBROUTINE svLS_LHS_CREATE 62*1e99f302SBen Matthews 63*1e99f302SBen Matthews SUBROUTINE svLS_LHS_FREE(lhs) 64*1e99f302SBen Matthews INCLUDE "svLS_STD.h" 65*1e99f302SBen Matthews TYPE(svLS_lhsType), INTENT(INOUT) :: lhs 66*1e99f302SBen Matthews END SUBROUTINE svLS_LHS_FREE 67*1e99f302SBen Matthews 68*1e99f302SBen Matthews SUBROUTINE svLS_LS_CREATE(ls, LS_type, relTol, absTol, maxItr, 69*1e99f302SBen Matthews 2 dimKry, relTolIn, absTolIn, maxItrIn) 70*1e99f302SBen Matthews INCLUDE "svLS_STD.h" 71*1e99f302SBen Matthews TYPE(svLS_lsType), INTENT(INOUT) :: ls 72*1e99f302SBen Matthews INTEGER, INTENT(IN) :: LS_type 73*1e99f302SBen Matthews REAL*8, INTENT(IN), OPTIONAL :: relTol, absTol, relTolIn(2), 74*1e99f302SBen Matthews 2 absTolIn(2) 75*1e99f302SBen Matthews INTEGER, INTENT(IN), OPTIONAL :: maxItr, dimKry, maxItrIn(2) 76*1e99f302SBen Matthews END SUBROUTINE svLS_LS_CREATE 77*1e99f302SBen Matthews 78*1e99f302SBen Matthews SUBROUTINE svLS_LS_FREE(ls) 79*1e99f302SBen Matthews INCLUDE "svLS_STD.h" 80*1e99f302SBen Matthews TYPE(svLS_lsType), INTENT(INOUT) :: ls 81*1e99f302SBen Matthews END SUBROUTINE svLS_LS_FREE 82*1e99f302SBen Matthews 83*1e99f302SBen Matthews SUBROUTINE svLS_BC_CREATE(lhs, faIn, nNo, dof, BC_type, 84*1e99f302SBen Matthews 2 gNodes, Val) 85*1e99f302SBen Matthews INCLUDE "svLS_STD.h" 86*1e99f302SBen Matthews TYPE(svLS_lhsType), INTENT(INOUT) :: lhs 87*1e99f302SBen Matthews INTEGER, INTENT(IN) :: faIn, nNo, dof 88*1e99f302SBen Matthews INTEGER, INTENT(IN) :: BC_type 89*1e99f302SBen Matthews INTEGER, INTENT(IN) :: gNodes(nNo) 90*1e99f302SBen Matthews REAL*8, INTENT(IN), OPTIONAL :: Val(dof,nNo) 91*1e99f302SBen Matthews END SUBROUTINE svLS_BC_CREATE 92*1e99f302SBen Matthews 93*1e99f302SBen Matthews SUBROUTINE svLS_BC_FREE(lhs, faIn) 94*1e99f302SBen Matthews INCLUDE "svLS_STD.h" 95*1e99f302SBen Matthews TYPE(svLS_lhsType), INTENT(INOUT) :: lhs 96*1e99f302SBen Matthews INTEGER, INTENT(IN) :: faIn 97*1e99f302SBen Matthews END SUBROUTINE svLS_BC_FREE 98*1e99f302SBen Matthews 99*1e99f302SBen Matthews SUBROUTINE svLS_SOLVE (lhs, ls, dof, Ri, Val, res) 100*1e99f302SBen Matthews INCLUDE "svLS_STD.h" 101*1e99f302SBen Matthews TYPE(svLS_lhsType), INTENT(INOUT) :: lhs 102*1e99f302SBen Matthews TYPE(svLS_lsType), INTENT(INOUT) :: ls 103*1e99f302SBen Matthews INTEGER, INTENT(IN) :: dof 104*1e99f302SBen Matthews REAL*8, INTENT(INOUT) :: Ri(dof,lhs%nNo) 105*1e99f302SBen Matthews REAL*8, INTENT(IN) :: Val(dof*dof,lhs%nnz) 106*1e99f302SBen Matthews REAL*8, INTENT(IN), OPTIONAL :: res(lhs%nFaces) 107*1e99f302SBen Matthews END SUBROUTINE svLS_SOLVE 108*1e99f302SBen Matthews 109*1e99f302SBen Matthews 110*1e99f302SBen Matthews 111*1e99f302SBen Matthews 112*1e99f302SBen Matthews 113