Polynomial2D example 1

This example shows the usage of [[Polynomial2D_]] class.

Modules and classes

• [[Polynomial2D_]]

Usage

PROGRAM main
use easifemBase
use easifemClasses
implicit none
class(Polynomial2D_), pointer :: f1=>NULL()
real(dfp), allocatable :: coeff( : )
integer(i4b), allocatable :: degree( :, : )

!!! note "Initiate" Initiate the object.

$$y=1+x+y+xy$$

  coeff = [1,1,1,1]
  call reallocate( degree, 4, 2 )
  degree(:,1) = [0,1,0,1]
  degree(:,2) = [0,0,1,1]
  f1=>Polynomial2D_Pointer( coeff=coeff, &
    & degree=degree, name1="x", name2="y" )
  call f1%display( 'f(x,y)=' )
  call f1%deallocate()

!!! example "result" f(x,y)=+1x^0y^1+1x^0y^0+1x^1y^1+1x^1y^0

!!! note "Initiate" Initiate the object.

$$y=1+2x+3y+4xy$$

  coeff = [1,2,3,4]
  call reallocate( degree, 4, 2 )
  degree(:,1) = [0,1,0,1]
  degree(:,2) = [0,0,1,1]
  f1=>Polynomial2D_Pointer( coeff=coeff, degree=degree, &
    & name1="x", name2="y" )
  call f1%display( 'f(x,y)=' )
  call f1%deallocate()

!!! example "result" f(x,y)=+1x^0y^1+1x^0y^0+1x^1y^1+1x^1y^0

!!! note "Initiate" Initiate the object.

$$y=1+2x+3y+4xy$$

  coeff = [1,2,3,2,2]
  call reallocate( degree, 5, 2 )
  degree(:,1) = [0,1,0,1,1]
  degree(:,2) = [0,0,1,1,1]
  f1=>Polynomial2D_Pointer( coeff=coeff, degree=degree, &
    & name1="x", name2="y" )
  call f1%display( 'f(x,y)=' )
  call f1%deallocate()

!!! example "result" f(x,y)=+1x^0y^1+1x^0y^0+1x^1y^1+1x^1y^0

END PROGRAM main