← Back to branch summary

~siesta-ts/siesta/trunk_ts_soc

« back to all changes in this revision

Viewing changes to Src/SiestaXC/mesh1d.f90

  • Committer: Nils Wittemeier
  • Date: 2019-02-14 07:45:07 UTC
  • mfrom: (746.1.15 trunk)
  • Revision ID: nils@4wittemeier.de-20190214074507-1mvzbmj9kw19gllr
Merged trunk 761

Show diffs side-by-side

added added

removed removed

Lines of Context:
Src/SiestaXC/mesh1d.f90
543
543
 
544
544
! Set the mesh points and the derivatives of x(i)
545
545
  if (mesh_type=='numerical') then
546
 
    xi = x
 
546
    xi(:) = x
547
547
 
548
548
!   Centered 5-point derivatives for all but first/last two points
549
549
    do i = 3,n-2
597
597
  if (present(xmax)) xi(n) = xmax
598
598
 
599
599
! Find auxiliary functions associated to the mesh
600
 
  sqrxp = abs(xp1)**0.5_dp
601
 
  s0 = abs(xp1)**1.5_dp
602
 
  s1 = xp1**2
603
 
  s2 = (3._dp*xp2**2 - 2._dp*xp1*xp3) / 4._dp / xp1**2
 
600
  sqrxp(:) = abs(xp1)**0.5_dp
 
601
  s0(:) = abs(xp1)**1.5_dp
 
602
  s1(:) = xp1**2
 
603
  s2(:) = (3._dp*xp2**2 - 2._dp*xp1*xp3) / 4._dp / xp1**2
604
604
 
605
605
  defined_mesh = .true.
606
606