All Downloads are FREE. Search and download functionalities are using the official Maven repository.

org.ejml.equation.ManagerFunctions Maven / Gradle / Ivy

Go to download

A fast and easy to use dense and sparse matrix linear algebra library written in Java.

There is a newer version: 0.43.1
Show newest version
/*
 * Copyright (c) 2022, Peter Abeles. All Rights Reserved.
 *
 * This file is part of Efficient Java Matrix Library (EJML).
 *
 * 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.ejml.equation;

import org.jetbrains.annotations.Nullable;

import java.util.HashMap;
import java.util.List;
import java.util.Map;

/**
 * Centralized place to create new instances of operations and functions. Must call
 * {@link #setManagerTemp} before any other functions.
 *
 * @author Peter Abeles
 */
@SuppressWarnings("NullAway.Init")
public class ManagerFunctions {

    // List of functions which take in N inputs
    Map input1 = new HashMap<>();
    Map inputN = new HashMap<>();

    // Reference to temporary variable manager
    protected ManagerTempVariables managerTemp;

    public ManagerFunctions() {
        addBuiltIn();
    }

    /**
     * Returns true if the string matches the name of a function
     */
    public boolean isFunctionName( String s ) {
        if (input1.containsKey(s))
            return true;
        return inputN.containsKey(s);
    }

    /**
     * Create a new instance of single input functions
     *
     * @param name function name
     * @param var0 Input variable
     * @return Resulting operation
     */
    public @Nullable Operation.Info create( String name, Variable var0 ) {
        Input1 func = input1.get(name);
        if (func == null)
            return null;
        return func.create(var0, managerTemp);
    }

    /**
     * Create a new instance of single input functions
     *
     * @param name function name
     * @param vars Input variables
     * @return Resulting operation
     */
    public @Nullable Operation.Info create( String name, List vars ) {
        InputN func = inputN.get(name);
        if (func == null)
            return null;
        return func.create(vars, managerTemp);
    }

    /**
     * Create a new instance of a single input function from an operator character
     *
     * @param op Which operation
     * @param input Input variable
     * @return Resulting operation
     */
    public Operation.Info create( char op, Variable input ) {
        switch (op) {
            case '\'':
                return Operation.transpose(input, managerTemp);

            default:
                throw new RuntimeException("Unknown operation " + op);
        }
    }

    /**
     * Create a new instance of a two input function from an operator character
     *
     * @param op Which operation
     * @param left Input variable on left
     * @param right Input variable on right
     * @return Resulting operation
     */
    public Operation.Info create( Symbol op, Variable left, Variable right ) {
        switch (op) {
            case PLUS:
                return Operation.add(left, right, managerTemp);

            case MINUS:
                return Operation.subtract(left, right, managerTemp);

            case TIMES:
                return Operation.multiply(left, right, managerTemp);

            case RDIVIDE:
                return Operation.divide(left, right, managerTemp);

            case LDIVIDE:
                return Operation.divide(right, left, managerTemp);

            case POWER:
                return Operation.pow(left, right, managerTemp);

            case ELEMENT_DIVIDE:
                return Operation.elementDivision(left, right, managerTemp);

            case ELEMENT_TIMES:
                return Operation.elementMult(left, right, managerTemp);

            case ELEMENT_POWER:
                return Operation.elementPow(left, right, managerTemp);

            default:
                throw new RuntimeException("Unknown operation " + op);
        }
    }

    public void setManagerTemp( ManagerTempVariables managerTemp ) {
        this.managerTemp = managerTemp;
    }

    /**
     * Adds a function, with a single input, to the list
     *
     * @param name Name of function
     * @param function Function factory
     */
    public void add1( String name, Input1 function ) {
        input1.put(name, function);
    }

    /**
     * Adds a function, with a two inputs, to the list
     *
     * @param name Name of function
     * @param function Function factory
     */
    public void addN( String name, InputN function ) {
        inputN.put(name, function);
    }

    /**
     * Adds built in functions
     */
    private void addBuiltIn() {
        input1.put("inv", Operation::inv);
        input1.put("pinv", Operation::pinv);
        input1.put("rref", Operation::rref);
        input1.put("eye", Operation::eye);
        input1.put("det", Operation::det);
        input1.put("normF", Operation::normF);
        input1.put("sum", Operation::sum_one);
        input1.put("trace", Operation::trace);
        input1.put("diag", Operation::diag);
        input1.put("min", Operation::min);
        input1.put("max", Operation::max);
        input1.put("abs", Operation::abs);
        input1.put("sin", Operation::sin);
        input1.put("cos", Operation::cos);
        input1.put("atan", Operation::atan);
        input1.put("exp", Operation::exp);
        input1.put("log", Operation::log);
        input1.put("sqrt", Operation::sqrt);
        input1.put("rng", Operation::rng);

        inputN.put("normP", ( inputs, manager ) -> {
            if (inputs.size() != 2) throw new RuntimeException("Two inputs expected");
            return Operation.normP(inputs.get(0), inputs.get(1), manager);
        });

        inputN.put("max", ( inputs, manager ) -> {
            if (inputs.size() != 2) throw new RuntimeException("One or two inputs expected");
            return Operation.max_two(inputs.get(0), inputs.get(1), manager);
        });

        inputN.put("min", ( inputs, manager ) -> {
            if (inputs.size() != 2) throw new RuntimeException("One or two inputs expected");
            return Operation.min_two(inputs.get(0), inputs.get(1), manager);
        });

        inputN.put("sum", ( inputs, manager ) -> {
            if (inputs.size() != 2) throw new RuntimeException("One or two inputs expected");
            return Operation.sum_two(inputs.get(0), inputs.get(1), manager);
        });

        inputN.put("zeros", ( inputs, manager ) -> {
            if (inputs.size() != 2) throw new RuntimeException("Two inputs expected");
            return Operation.zeros(inputs.get(0), inputs.get(1), manager);
        });

        inputN.put("ones", ( inputs, manager ) -> {
            if (inputs.size() != 2) throw new RuntimeException("Two inputs expected");
            return Operation.ones(inputs.get(0), inputs.get(1), manager);
        });

        inputN.put("rand", ( inputs, manager ) -> {
            if (inputs.size() != 2) throw new RuntimeException("Two inputs expected");
            return Operation.rand(inputs.get(0), inputs.get(1), manager);
        });

        inputN.put("randn", ( inputs, manager ) -> {
            if (inputs.size() != 2) throw new RuntimeException("Two inputs expected");
            return Operation.randn(inputs.get(0), inputs.get(1), manager);
        });

        inputN.put("kron", ( inputs, manager ) -> {
            if (inputs.size() != 2) throw new RuntimeException("Two inputs expected");
            return Operation.kron(inputs.get(0), inputs.get(1), manager);
        });

        inputN.put("dot", ( inputs, manager ) -> {
            if (inputs.size() != 2) throw new RuntimeException("Two inputs expected");
            return Operation.dot(inputs.get(0), inputs.get(1), manager);
        });

        inputN.put("pow", ( inputs, manager ) -> {
            if (inputs.size() != 2) throw new RuntimeException("Two inputs expected");
            return Operation.pow(inputs.get(0), inputs.get(1), manager);
        });

        inputN.put("atan2", ( inputs, manager ) -> {
            if (inputs.size() != 2) throw new RuntimeException("Two inputs expected");
            return Operation.atan2(inputs.get(0), inputs.get(1), manager);
        });

        inputN.put("solve", ( inputs, manager ) -> {
            if (inputs.size() != 2) throw new RuntimeException("Two inputs expected");
            return Operation.solve(inputs.get(0), inputs.get(1), manager);
        });

        inputN.put("extract", Operation::extract);
        inputN.put("extractScalar", ( inputs, manager ) -> {
            if (inputs.size() != 2 && inputs.size() != 3) throw new RuntimeException("Two or three inputs expected");
            return Operation.extractScalar(inputs, manager);
        });
    }

    public ManagerTempVariables getManagerTemp() {
        return managerTemp;
    }

    /**
     * Creates new instances of functions from a single variable
     */
    public interface Input1 {
        Operation.Info create( Variable A, ManagerTempVariables manager );
    }

    /**
     * Creates a new instance of functions from two variables
     */
    public interface InputN {
        Operation.Info create( List inputs, ManagerTempVariables manager );
    }
}




© 2015 - 2024 Weber Informatics LLC | Privacy Policy