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