223031012
/
DAMASK_EICMD
mirror of https://github.com/achalhp/DAMASK_EICMD.git
1
0
Fork 0
Code Issues Projects Releases Wiki Activity

generic code for variable order

This commit is contained in:
Martin Diehl 2022-05-07 22:53:31 +02:00
parent 70b9943920
commit 71e4fa222c
1 changed files with 9 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
src/polynomials.f90
View File

@ -70,7 +70,8 @@ function polynomial_from_dict(dict,y,x) result(p)
real(pReal), dimension(:), allocatable :: coef real(pReal), dimension(:), allocatable :: coef
real(pReal) :: x_ref real(pReal) :: x_ref
integer :: i integer :: i, o
character(len=1) :: o_s
allocate(coef(1),source=dict%get_asFloat(y)) allocate(coef(1),source=dict%get_asFloat(y))
@ -79,10 +80,13 @@ function polynomial_from_dict(dict,y,x) result(p)
x_ref = dict%get_asFloat(x//'_ref') x_ref = dict%get_asFloat(x//'_ref')
coef = [coef,dict%get_asFloat(y//','//x)] coef = [coef,dict%get_asFloat(y//','//x)]
end if end if
if (dict%contains(y//','//x//'^2')) then do o = 2,2
x_ref = dict%get_asFloat(x//'_ref') write(o_s,'(I0.0)') o
coef = [coef,[(0.0_pReal,i=size(coef),2-1)],dict%get_asFloat(y//','//x//'^2')] if (dict%contains(y//','//x//'^'//o_s)) then
end if x_ref = dict%get_asFloat(x//'_ref')
coef = [coef,[(0.0_pReal,i=size(coef),o-1)],dict%get_asFloat(y//','//x//'^'//o_s)]
end if
end do
p = Polynomial(coef,x_ref) p = Polynomial(coef,x_ref)