Révision 9 src/Egrad.f90
Egrad.f90 (revision 9) | ||
---|---|---|
454 | 454 |
END DO |
455 | 455 |
DEALLOCATE(GradTmp) |
456 | 456 |
CASE ("CART","HYBRID") |
457 |
! PFL 2013 Feb |
|
458 |
! We could think of just doing: |
|
459 |
! Grad=GradCart |
|
460 |
! however, this does not work as implicitly we have : |
|
461 |
! Grad(Nat,3) BUT GradCar(3,Nat) |
|
462 |
! So we have to use reshape :-( |
|
463 |
Grad=reshape(reshape(GradCart,(/3,Nat/),ORDER=(/2,1/)),(/3*Nat/)) |
|
457 |
Grad=GradCart |
|
464 | 458 |
CASE DEFAULT |
465 | 459 |
WRITE(*,*) "Coord=",AdjustL(Trim(COORD))," not yet implemented in Egradients. STOP" |
466 | 460 |
STOP |
Formats disponibles : Unified diff