Skip to content
Merged
Show file tree
Hide file tree
Changes from 1 commit
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
Prev Previous commit
update changelog and fix bug in wheels
  • Loading branch information
prehner committed Jul 13, 2022
commit 2ac40f13914127653fa7a9e25baebf8a5afb1a28
3 changes: 3 additions & 0 deletions feos-dft/CHANGELOG.md
Original file line number Diff line number Diff line change
Expand Up @@ -13,6 +13,9 @@ and this project adheres to [Semantic Versioning](https://semver.org/spec/v2.0.0
- Made FMT functional more flexible w.r.t. the shape of the weight functions. [#31](https://github.com/feos-org/feos-dft/pull/31)
- Changed interface of `PairCorrelationFunction` to facilitate the calculation of pair correlation functions in mixtures. [#29](https://github.com/feos-org/feos-dft/pull/29)

### Removed
- Moved the implementation of fundamental measure theory to the new `feos-saft` crate. [#27](https://github.com/feos-org/feos/pull/27)

## [0.2.0] - 2022-04-12
### Added
- Added `grand_potential_density` getter for DFT profiles in Python. [#22](https://github.com/feos-org/feos-dft/pull/22)
Expand Down
2 changes: 1 addition & 1 deletion src/pcsaft/dft/mod.rs
Original file line number Diff line number Diff line change
Expand Up @@ -150,7 +150,7 @@ impl PairPotential for PcSaftFunctional {
fn pair_potential(&self, i: usize, r: &Array1<f64>, _: f64) -> Array2<f64> {
let sigma_ij = &self.parameters.sigma_ij;
let eps_ij_4 = 4.0 * &self.parameters.epsilon_k_ij;
Array::from_shape_fn((self.parameters.m.len(), r.len()), |(j, k)| {
Array2::from_shape_fn((self.parameters.m.len(), r.len()), |(j, k)| {
let att = (sigma_ij[[i, j]] / r[k]).powi(6);
eps_ij_4[[i, j]] * att * (att - 1.0)
})
Expand Down
2 changes: 1 addition & 1 deletion src/pets/dft/mod.rs
Original file line number Diff line number Diff line change
Expand Up @@ -140,7 +140,7 @@ impl PairPotential for PetsFunctional {
let eps_ij_4 = 4.0 * self.parameters.epsilon_k_ij.clone();
let shift_ij = &eps_ij_4 * (2.5.powi(-12) - 2.5.powi(-6));
let rc_ij = 2.5 * &self.parameters.sigma_ij;
Array::from_shape_fn((self.parameters.sigma.len(), r.len()), |(j, k)| {
Array2::from_shape_fn((self.parameters.sigma.len(), r.len()), |(j, k)| {
if r[k] > rc_ij[[i, j]] {
0.0
} else {
Expand Down