/* * Licensed to the Apache Software Foundation (ASF) under one * or more contributor license agreements. See the NOTICE file * distributed with this work for additional information * regarding copyright ownership. The ASF licenses this file * to you 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.apache.sysml.runtime.instructions.cp; import org.apache.sysml.parser.Expression.DataType; import org.apache.sysml.parser.Expression.ValueType; import org.apache.sysml.runtime.DMLRuntimeException; import org.apache.sysml.runtime.controlprogram.context.ExecutionContext; import org.apache.sysml.runtime.functionobjects.CM; import org.apache.sysml.runtime.instructions.InstructionUtils; import org.apache.sysml.runtime.matrix.data.MatrixBlock; import org.apache.sysml.runtime.matrix.operators.CMOperator; import org.apache.sysml.runtime.matrix.operators.CMOperator.AggregateOperationTypes; public class CentralMomentCPInstruction extends AggregateUnaryCPInstruction { public CentralMomentCPInstruction(CMOperator cm, CPOperand in1, CPOperand in2, CPOperand in3, CPOperand out, String opcode, String str) { super(cm, in1, in2, in3, out, opcode, str); } public static CentralMomentCPInstruction parseInstruction(String str) throws DMLRuntimeException { CPOperand in1 = new CPOperand("", ValueType.UNKNOWN, DataType.UNKNOWN); CPOperand in2 = null; CPOperand in3 = null; CPOperand out = new CPOperand("", ValueType.UNKNOWN, DataType.UNKNOWN); String[] parts = InstructionUtils.getInstructionPartsWithValueType(str); String opcode = parts[0]; //check supported opcode if( !opcode.equalsIgnoreCase("cm") ) { throw new DMLRuntimeException("Unsupported opcode "+opcode); } if ( parts.length == 4 ) { // Example: CP.cm.mVar0.Var1.mVar2; (without weights) in2 = new CPOperand("", ValueType.UNKNOWN, DataType.UNKNOWN); parseUnaryInstruction(str, in1, in2, out); } else if ( parts.length == 5) { // CP.cm.mVar0.mVar1.Var2.mVar3; (with weights) in2 = new CPOperand("", ValueType.UNKNOWN, DataType.UNKNOWN); in3 = new CPOperand("", ValueType.UNKNOWN, DataType.UNKNOWN); parseUnaryInstruction(str, in1, in2, in3, out); } /* * Exact order of the central moment MAY NOT be known at compilation time. * We first try to parse the second argument as an integer, and if we fail, * we simply pass -1 so that getCMAggOpType() picks up AggregateOperationTypes.INVALID. * It must be updated at run time in processInstruction() method. */ int cmOrder; try { if ( in3 == null ) { cmOrder = Integer.parseInt(in2.getName()); } else { cmOrder = Integer.parseInt(in3.getName()); } } catch(NumberFormatException e) { cmOrder = -1; // unknown at compilation time } AggregateOperationTypes opType = CMOperator.getCMAggOpType(cmOrder); CMOperator cm = new CMOperator(CM.getCMFnObject(opType), opType); return new CentralMomentCPInstruction(cm, in1, in2, in3, out, opcode, str); } @Override public void processInstruction( ExecutionContext ec ) throws DMLRuntimeException { String output_name = output.getName(); /* * The "order" of the central moment in the instruction can * be set to INVALID when the exact value is unknown at * compilation time. We first need to determine the exact * order and update the CMOperator, if needed. */ MatrixBlock matBlock = ec.getMatrixInput(input1.getName()); CPOperand scalarInput = (input3==null ? input2 : input3); ScalarObject order = ec.getScalarInput(scalarInput.getName(), scalarInput.getValueType(), scalarInput.isLiteral()); CMOperator cm_op = ((CMOperator)_optr); if ( cm_op.getAggOpType() == AggregateOperationTypes.INVALID ) { ((CMOperator)_optr).setCMAggOp((int)order.getLongValue()); } CM_COV_Object cmobj = null; if (input3 == null ) { cmobj = matBlock.cmOperations(cm_op); } else { MatrixBlock wtBlock = ec.getMatrixInput(input2.getName()); cmobj = matBlock.cmOperations(cm_op, wtBlock); ec.releaseMatrixInput(input2.getName()); } ec.releaseMatrixInput(input1.getName()); double val = cmobj.getRequiredResult(_optr); DoubleObject ret = new DoubleObject(output_name, val); ec.setScalarOutput(output_name, ret); } }