*DECKSDSDOTFREALFUNCTIONSDSDOTF(N,SB,SX,INCX,SY,INCY)C***BEGINPROLOGUESDSDOTC***PURPOSEComputetheinnerproductoftwovectorswithextendedCprecisionaccumulation.C***LIBRARYSLATEC(BLAS)C***CATEGORYD1A4C***TYPESINGLEPRECISION(SDSDOT-S,CDCDOT-C)C***KEYWORDSBLAS,DOTPRODUCT,INNERPRODUCT,LINEARALGEBRA,VECTORC***AUTHORLawson,C.L.,(JPL)CHanson,R.J.,(SNLA)CKincaid,D.R.,(U.ofTexas)CKrogh,F.T.,(JPL)C***DESCRIPTIONCCBLASSubprogramCDescriptionofParametersCC--Input--CNnumberofelementsininputvector(s)CSBsingleprecisionscalartobeaddedtoinnerproductCSXsingleprecisionvectorwithNelementsCINCXstoragespacingbetweenelementsofSXCSYsingleprecisionvectorwithNelementsCINCYstoragespacingbetweenelementsofSYCC--Output--CSDSDOTsingleprecisiondotproduct(SBifN.LE.0)CCReturnsS.P.resultwithdotproductaccumulatedinD.P.CSDSDOT=SB+sumforI=0toN-1ofSX(LX+I*INCX)*SY(LY+I*INCY),CwhereLX=1ifINCX.GE.0,elseLX=1+(1-N)*INCX,andLYisCdefinedinasimilarwayusingINCY.CC***REFERENCESC.L.Lawson,R.J.Hanson,D.R.KincaidandF.T.CKrogh,BasiclinearalgebrasubprogramsforFortranCusage,AlgorithmNo.539,TransactionsonMathematicalCSoftware5,3(September1979),pp.308-323.C***ROUTINESCALLED(NONE)C***REVISIONHISTORY(YYMMDD)C791001DATEWRITTENC890531Changedallspecificintrinsicstogeneric.(WRB)C890831Modifiedarraydeclarations.(WRB)C890831REVISIONDATEfromVersion3.2C891214PrologueconvertedtoVersion4.0format.(BAB)C920310CorrecteddefinitionofLXinDESCRIPTION.(WRB)C920501ReformattedtheREFERENCESsection.(WRB)C***ENDPROLOGUESDSDOTREALSX(*),SY(*),SBDOUBLE PRECISIONDSDOTC***FIRSTEXECUTABLESTATEMENTSDSDOTDSDOT=SBIF(N.LE.0)GO TO30IF(INCX.EQ.INCY.AND.INCX.GT.0)GO TO40CCCodeforunequalornonpositiveincrements.CKX=1KY=1IF(INCX.LT.0)KX=1+(1-N)*INCXIF(INCY.LT.0)KY=1+(1-N)*INCYDO10I=1,NDSDOT=DSDOT+DBLE(SX(KX))*DBLE(SY(KY))KX=KX+INCXKY=KY+INCY10CONTINUE30SDSDOTF=DSDOTRETURNCCCodeforequalandpositiveincrements.C40NS=N*INCXDO50I=1,NS,INCXDSDOT=DSDOT+DBLE(SX(I))*DBLE(SY(I))50CONTINUESDSDOTF=DSDOTRETURNEND
此处可能存在不合适展示的内容,页面不予展示。您可通过相关编辑功能自查并修改。
如您确认内容无涉及 不当用语 / 纯广告导流 / 暴力 / 低俗色情 / 侵权 / 盗版 / 虚假 / 无价值内容或违法国家有关法律法规的内容,可点击提交进行申诉,我们将尽快为您处理。