基因数据处理120之scala调用SSW在linux下运行

更多代码请见:https://github.com/xubo245

基因数据处理系列

1.解释

先有java提供转换,使用jni调用c

然后scala调用java

2.代码:

2.1 java:

package ssw;

/**
 * Created by xubo on 2016/11/25.
 */
public class SSW {
    public static int align(String query, String ref) {
        try {
            System.loadLibrary("sswjni");
        } catch (UnsatisfiedLinkError e) {
            System.out.println(String.format("Cannot find libsswjni.so. Has the library been built and LD_LIBRARY_PATH or -Djava.library.path set appropriately?\n%s", e));
            throw e;
        }
        int[][] score = new int[128][128];
        for (int i = 0; i < 128; i++) {
            for (int j = 0; j < 128; j++) {
                if (i == j) score[i][j] = 2;
                else score[i][j] = -2;
            }
        }
        System.out.println("Aligning nucleotides");
        Alignment aln = Aligner.align(query.getBytes(), ref.getBytes(), score, 12, 2, true);
        if (aln == null) {
            throw new RuntimeException();
        }

        return aln.score1;
    }

    public static int[][] blosum50() {
        int[][] a = {
                {5, -2, -1, -2, -1, -3, 0, -2, -1, -2, -1, -2, -1, -1, -5, -1, -1, -2, 1, 0, -5, 0, -3, -1, -2, -1},
                {-2, 6, -3, 6, 1, -4, -1, 0, -4, -4, 0, -4, -3, 5, -5, -2, 0, -1, 0, 0, -5, -3, -5, -1, -3, 1},
                {-1, -3, 13, -4, -3, -2, -3, -3, -2, -2, -3, -2, -2, -2, -5, -4, -3, -4, -1, -1, -5, -1, -5, -1, -3, -3},
                {-2, 6, -4, 8, 2, -5, -1, -1, -4, -4, -1, -4, -4, 2, -5, -1, 0, -2, 0, -1, -5, -4, -5, -1, -3, 1},
                {-1, 1, -3, 2, 6, -3, -3, 0, -4, -3, 1, -3, -2, 0, -5, -1, 2, 0, -1, -1, -5, -3, -3, -1, -2, 5},
                {-3, -4, -2, -5, -3, 8, -4, -1, 0, 1, -4, 1, 0, -4, -5, -4, -4, -3, -3, -2, -5, -1, 1, -1, 4, -4},
                {0, -1, -3, -1, -3, -4, 8, -2, -4, -4, -2, -4, -3, 0, -5, -2, -2, -3, 0, -2, -5, -4, -3, -1, -3, -2},
                {-2, 0, -3, -1, 0, -1, -2, 10, -4, -3, 0, -3, -1, 1, -5, -2, 1, 0, -1, -2, -5, -4, -3, -1, 2, 0},
                {-1, -4, -2, -4, -4, 0, -4, -4, 5, 4, -3, 2, 2, -3, -5, -3, -3, -4, -3, -1, -5, 4, -3, -1, -1, -3},
                {-2, -4, -2, -4, -3, 1, -4, -3, 4, 4, -3, 4, 2, -4, -5, -3, -3, -3, -3, -1, -5, 2, -2, -1, -1, -3},
                {-1, 0, -3, -1, 1, -4, -2, 0, -3, -3, 6, -3, -2, 0, -5, -1, 2, 3, 0, -1, -5, -3, -3, -1, -2, 1},
                {-2, -4, -2, -4, -3, 1, -4, -3, 2, 4, -3, 5, 3, -4, -5, -4, -2, -3, -3, -1, -5, 1, -2, -1, -1, -3},
                {1, -3, -2, -4, -2, 0, -3, -1, 2, 2, -2, 3, 7, -2, -5, -3, 0, -2, -2, -1, -5, 1, -1, -1, 0, -1},
                {-1, 5, -2, 2, 0, -4, 0, 1, -3, -4, 0, -4, -2, 7, -5, -2, 0, -1, 1, 0, -5, -3, -4, -1, -2, 0},
                {-5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, 1, -5, -5, -5, -5, -5, -5, -5, -5, -5, 0, -5},
                {-1, -2, -4, -1, -1, -4, -2, -2, -3, -3, -1, -4, -3, -2, -5, 10, -1, -3, -1, -1, -5, -3, -4, -1, -3, -1},
                {-1, 0, -3, 0, 2, -4, -2, 1, -3, -3, 2, -2, 0, 0, -5, -1, 7, 1, 0, -1, -5, -3, -1, -1, -1, 4},
                {-2, -1, -4, -2, 0, -3, -3, 0, -4, -3, 3, -3, -2, -1, -5, -3, 1, 7, -1, -1, -5, -3, -3, -1, -1, 0},
                {1, 0, -1, 0, -1, -3, 0, -1, -3, -3, 0, -3, -2, 1, -5, -1, 0, -1, 5, 2, -5, -2, -4, -1, -2, 0},
                {0, 0, -1, -1, -1, -2, -2, -2, -1, -1, -1, -1, -1, 0, -5, -1, -1, -1, 2, 5, 1, 0, -3, -1, -2, -1},
                {-2, -1, -2, -3, -3, -1, -2, -3, 2, -1, -1, -2, 0, 4, 0, -3, -2, -2, 2, 8, 0, -1, -3, -1, -2, -1},
                {0, -3, -1, -4, -3, -1, -4, -4, 4, 2, -3, 1, 1, -3, -5, -3, -3, -3, -2, 0, -5, 5, -3, -1, -1, -3},
                {-3, -5, -5, -5, -3, 1, -3, -3, -3, -2, -3, -2, -1, -4, -5, -4, -1, -3, -4, -3, -5, -3, 15, -1, 2, -2},
                {-1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -5, -1, -1, -1, -1, -1, -5, -1, -1, -1, -1, -1},
                {-5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, -5, 8, -5},
                {-1, 1, -3, 1, 5, -4, -2, 0, -3, -3, 1, -3, -1, 0, -5, -1, 4, 0, 0, -1, -5, -3, -2, -1, -2, 5}
        };
        return a;
    }
}

2.2 scala调用:

package org.dsw.other.ssw

import ssw.SSW

/**
  * Created by xubo on 2016/11/25.
  */
object SSWLearning {
  def main(args: Array[String]) {

    val score = SSW.align("AGCT", "ACT")
    println(score)

  }
}

3.结果:

xubo@xubo:~/xubo/tools/Complete-Striped-Smith-Waterman-Library/src$ scala  -Djava.library.path=. -cp SparkSW.jar org.dsw.other.ssw.SSWLearning
Aligning nucleotides
4

参考

【1】https://github.com/xubo245
【2】http://blog.csdn.net/xubo245/
一、Scala核心编程课程简介 近年来随着大数据的兴起,大数据核心框架Spark和Kafka也受到到广泛关注, Spark底层是Scala实现的, 因此也带火了Scala语言,目前Scala被全球知名公司(如:谷歌、百度、阿里、Twitter、京东等)广泛用于Spark开发。新一代的物联网时代到来,会对大数据应用人才的需求越加紧迫。 尚硅谷网罗和整合了学员很喜爱的师资,打造出专注于Spark开发的Scala课程,我们讲解Scala核心编程技术,同时也讲解编程思想、设计模式和Scala底层实现,让您有豁然开朗的感受。 二、课程内容和目标 本课程重点讲解Scala核心编程,内容包括: Scala语言概述、运算符、程序流程控制、数据结构之集合、Map映射、过滤、化简、折叠、扫描、拉链、视图、并行集合、高阶函数、函数柯里化、偏函数、参数推断、控制抽象、Trait、面向对象编程、异常处理、惰性函数、Akka及Actor模型、Spark Master和Worker通讯、隐式转换、隐式参数、工厂模式、单例模式、观察者模式、装饰者模式、代理模式、泛型、上下界、视图界定、上下文界定、协变逆变不变和源码剖析。通过系统全面的学习,学员能掌握Scala编程思想和Scala底层机制,为进一步学习Spark打下坚实基础。 三、谁适合学 1.希望以较低的投入和风险,来了解自己是否适合转型从事Spark开发的求职人员。 2.有一定的Java基础,或者自学过一些Java书籍与视频资料,想系统掌握Scala语言的程序员。
©️2020 CSDN 皮肤主题: 大白 设计师:CSDN官方博客 返回首页