org.wallerlab.yoink.adaptive.smooth.smoothfunction.SCMPSmoothFunction.java Source code

Introduction

Here is the source code for org.wallerlab.yoink.adaptive.smooth.smoothfunction.SCMPSmoothFunction.java

Source

/*
 * Copyright 2014-2015 the original author or authors.
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 * http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 */

package org.wallerlab.yoink.adaptive.smooth.smoothfunction;

import org.springframework.stereotype.Service;
import org.wallerlab.yoink.api.service.adaptive.SmoothFunction;

/**
 * this class is for one smooth function used in SCMP for QM lambda.
 * 
 * @author Min Zheng
 *
 */
@Service("scmpSmoothFunction")
public class SCMPSmoothFunction implements SmoothFunction {

    /**
     * this smooth function is used in SCMP method. for details please see:
     * "Size-Consistent Multipartitioning QM/MM: A Stable and Efficient Adaptive
     * QM/MM Method"
     * 
     * @param currentValue
     *            , currentValue(variable) in smooth function
     * @param min
     *            , minimum value in smooth function
     * @param max
     *            , maximum value in smooth function
     * @return smooth factor
     */
    public double evaluate(double currentValue, double min, double max) {
        double smoothFactor;
        if (currentValue > max) {
            smoothFactor = 0;
        } else if (currentValue < min) {
            smoothFactor = 1;
        } else {
            smoothFactor = Math.pow((currentValue - max), 2);
            smoothFactor = smoothFactor * (-3 * min + max + 2 * currentValue);
            smoothFactor = smoothFactor / Math.pow((max - min), 3);
        }
        return smoothFactor;
    }

}