LCOV - code coverage report
Current view: top level - math - difcub.f (source / functions) Hit Total Coverage
Test: combined.info Lines: 0 2 0.0 %
Date: 2019-09-08 04:53:50 Functions: 0 1 0.0 %

          Line data    Source code
       1             : !--------------------------------------------------------------------------------
       2             : ! Copyright (c) 2017 Peter Grünberg Institut, Forschungszentrum Jülich, Germany
       3             : ! This file is part of FLEUR and available as free software under the conditions
       4             : ! of the MIT license as expressed in the LICENSE file in more detail.
       5             : !--------------------------------------------------------------------------------
       6             : 
       7             :       MODULE m_difcub
       8             : 
       9             :       CONTAINS
      10             : 
      11           0 :       REAL FUNCTION difcub(x,f,xi)
      12             : c     **********************************************************
      13             : c     differentiate the function f, given at the
      14             : c     points x0,x1,x2,x3 at the point xi by lagrange
      15             : c     interpolation for polynomial of 3rd order
      16             : c     r.p.
      17             : c     ***********************************************************
      18             :       IMPLICIT NONE
      19             : C     .. Scalar Arguments ..
      20             :       REAL xi
      21             : C     ..
      22             : C     .. Array Arguments ..
      23             :       REAL f(0:3),x(0:3)
      24             : C     ..
      25             :       difcub = ((xi-x(1))* (xi-x(2))+ (xi-x(1))* (xi-x(3))+
      26             :      +         (xi-x(2))* (xi-x(3)))*f(0)/ ((x(0)-x(1))* (x(0)-x(2))*
      27             :      +         (x(0)-x(3))) + ((xi-x(0))* (xi-x(2))+
      28             :      +         (xi-x(0))* (xi-x(3))+ (xi-x(2))* (xi-x(3)))*f(1)/
      29             :      +         ((x(1)-x(0))* (x(1)-x(2))* (x(1)-x(3))) +
      30             :      +         ((xi-x(0))* (xi-x(1))+ (xi-x(0))* (xi-x(3))+
      31             :      +         (xi-x(1))* (xi-x(3)))*f(2)/ ((x(2)-x(0))* (x(2)-x(1))*
      32             :      +         (x(2)-x(3))) + ((xi-x(0))* (xi-x(1))+
      33             :      +         (xi-x(0))* (xi-x(2))+ (xi-x(1))* (xi-x(2)))*f(3)/
      34           0 :      +         ((x(3)-x(0))* (x(3)-x(1))* (x(3)-x(2)))
      35             :       RETURN
      36             :       END FUNCTION difcub
      37             : 
      38             :       END MODULE m_difcub

Generated by: LCOV version 1.13