ForwardDifference.h

Go to the documentation of this file.

/* Copyright 2018 Kristofer Björnson
 *
 * 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.
 */
 
#ifndef COM_DAFER45_TBTK_FORWARD_DIFFERENCE
#define COM_DAFER45_TBTK_FORWARD_DIFFERENCE
 
#include "TBTK/HoppingAmplitudeList.h"
 
namespace TBTK{
 
class ForwardDifference : public HoppingAmplitudeList{
public:
    ForwardDifference(
        unsigned int subindex,
        const Index &index,
        double dx = 1
    );
};
 
inline ForwardDifference::ForwardDifference(
    unsigned int subindex,
    const Index &index,
    double dx
){
    TBTKAssert(
        subindex < index.getSize(),
        "ForwardDifference::ForwardDifference()",
        "Invalid subindex. The subindex '" << subindex << "' is larger"
        << " than the size of the Index '" << index.getSize() << "'.",
        ""
    );
    Index forward = index;
    forward[subindex]++;
    add(HoppingAmplitude(1/dx, index, index));
    add(HoppingAmplitude(-1/dx, index, forward));
}
 
};  //End of namespace TBTK
 
#endif